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I ABSTRACT 

This  report  describes  various  techniques  for  performing 
linear,  space-variant  operations  with  coherent  optical  pro- 
cessors. General  linear  systems  theory  is  developed  for  appli- 
cation to  processor  design.  Three  analytical  characterizations 
t are  separately  examined  based  on  (a)  a continuum  orthonormal 

basis  set  response;  (b)  a piecewise  isoplanatic  approximation 
and  (c)  modifications  of  the  Whittaker-Shannon  sampling  theorem 
for  space-variant  systems.  Advantages  and  disadvantages  of 
each  approach  are  considered.  Generalized  one-dimensional 
(1-D)  space-variant  coherent  optical  processors  capable  of  per- 
t forming  a broad  class  of  1-D  superposition  integral  operations 

are  presented.  Applications  discussed  include  magnification, 
coordinate  distortion,  convolution,  correlation,  cross-power 
spectrum  display,  cimbiguity  function  display,  and  evaluation 
of  Laplace  and  inverse  Abel  transforms.  A technique  for  lens- 
less 1-D  space-variant  processing  is  also  presented  wherein  a 

• one-dimensional  linear  operation  can  be  performed  with  an 
input,  a processor  mask,  and  a few  centimeters  of  free  space. 
Two  dimensional  space-variant  processors  are  discussed  from 
the  frameworks  of  (a)  a composite  hologram  scheme  resulting 
from  a direct  application  of  a space-variant  system  sampling 
theorem  and  (b)  angle-multiplexing  of  a number  of  hologreims 

• in  a thick  recording  medium.  An  improvement  on  the  volume 
hologram  approach  based  on  the  use  of  phase-coding  diffusers 
in  the  multiple  reference  beams  is  also  discussed. 
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Chapter  I 


1 . INTRODUCTION 

Initial  recognition  of  the  ability  to  perform  general 
mathematical  operations  with  coherent  light  can  be  attribu- 
ted to  the  members  of  the  research  team  involved  in  the  so 

1-2 

called  Project  Michigan.  In  initially  classified  re- 
search, they  demonstrated  the  feasibility  of  utilizing 
coherent  optical  processors  to  handle  vast  quantities  of 
data  at  the  speed  of  light.  Coherent  processors  can  thus 
perform  operations  at  a rate  far  in  excess  of  that  re- 
quired by  conventional  digital  computers.  Being  an  analog 
computer,  however,  the  coherent  processor  suffers  from 
inflexibility  and  inexactitude  when  compared  with  the 
digital  computer.  The  coherent  processor,  though,  has  a 
far  greater  capacity  for  information  throughput. 

There  are  two  basic  building  blocks  for  any  coherent 
processor.  The  first  is  the  thin  lens  which,  besides 
familiar  imaging  operations,  is  capable  of  performing  a 
Fourier  transformation  on  incident  light  field  eunplitudes. 
The  second  is  the  hologram  on  which  complex  field  eunplitude 
distributions  can  be  represented.  The  most  basic  processor 
utilizing  these  components  is  shown  in  Figure  1-1.  A two- 
dimensional  transmittance,  f,  is  illuminated  from  the  left 
by  a coherent  plane  wave.  The  Fourier  transform  of  another 
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A basic  coherent  optical  processor  for  per- 
forming space-invariant  (isoplanatic,  con- 
volution) operations. 
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function,  say  h,  is  placed  in  plane  P2.  The  convolution 
of  h and  f then  appears  on  the  processor  output  plane > 

The  convolution  operation  performed  here  is  also  referred 
to  as  a linear  space- invariant  (or  isoplanatic)  operation. 
The  space- invariant  operation  is  a subclass  of  the  more 
general  class  of  linear  operations.  If  a linear  operation 
is  not  space-invariant,  then  it  is  said  to  be  space- variant. 
A more  rigorous  treatment  of  these  concepts  is  presented 
in  Chapter  II. 

The  use  of  the  term  "space-variant"  in  the  report 
title  is  somewhat  inexact  and  is  used  not  because  of  the 
author's  preference  but  rather  due  to  precedence.  Due  to 
the  wide  spread  use  of  the  coherent  processor  in  Figure 
1-1  and  the  predominance  of  time-invariant  systems  in 
the  sister  electrical  sciences,  the  term  "space- (time-) 
variance"  has  been  used  to  stress  that  a linear  operation 
is  "not  space-invariant."  There  are  also  those  to  whom 
the  term  "linear"  implies  "invariant"  although  such  is 
not  rigorously  the  case.  Following  this  unfortunate 
historical  development,  we  herein  also  use  the  term 
"space-variant"  although  the  term  "linear"  (space- variant 
and  space- invariant)  is  rigorously  more  appropriate  in 
most  cases. 

The  class  of  space- invariant  operations  that  can  be 
performed  by  the  coherent  processor  in  Figure  1-1  is  fun- 
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daunentally  limited  only  by  the  implementability  of  the 
so  called  transfer  function  mask  on  plane  P2.  More  ex- 
citedly, any  coherent  processor  which  perfoinns  a two-di- 
mensional space- invariant  operation,  irregardless  of  the 
number  of  components  or  complexity,  can,  in  principle,  be 
reduced  to  the  elementary  configuration  in  Figure  1-1. 

Unfortunately,  most  coherent  processors  are,  in  fact, 
not  space- invariant.  Many,  however,  are  linear  in  nature. 

As  such,  it  would  be  a useful  undertaking  to  develop  a 
space-variant  coherent  processor  displaying  the  simplicity 
of  Figure  1-1  to  which  all  linear  coherent  processors 
could  be  reduced.  Investigation  of  this  possibility  is 
the  main  topic  of  this  report. 

There  are  of  course,  many  specific  coherent  processors 
which  perform  space-variant  operations.  The  thin  lens 
Fourier  transformer  and  non-unity  magnification  imaging 
systems  are  prime  examples.  In  our  investigations,  how- 
ever, we  are  more  concerned  with  developing  a general 
processor  capable  of  performing  a wide  class  of  opera- 
tions. 

We  begin  treatment  of  processor  design  in  Chapter  II 
where  several  aspects  of  linear  system  theory  are  developed. 
System  characterization  by  determination  of  the  system  re- 
sponse to  all  elements  in  a continuum  orthonormal  basis 
set  is  presented  and  leads  to  the  well  known  line  spread 
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function  and  frequency  response  characterizations  as  special 
cases.  The  "ontinuum  orthonormal  basis  set  response  char- 
acterization brings  to  light  the  possibility  of  alternate 
system  descriptions  which  allow  for  the  enrichment  of  lin- 
ear systems  theory  beyond  conventional  impulse  and  fre- 

3-4 

quency  response  characterization.  Furthermore,  it  is 
shown  that  due  to  the  choice  of  our  system  response  no- 
tation,^ Fourier  analysis,  conventionally  utilized  only 
for  space- invariant  system  description,  can  straight- 
forwardly be  generalized  for  description  of  the  more  gen- 
eral class  of  linear  systems. 

The  remainder  of  Chapter  II  is  concerned  with  system 
characterization  wherein  certain  restrictive  assumptions 
are  placed  on  the  system  and/or  input.  The  result  of 
these  constraints  in  all  cases  considered  is  the  reduction 
of  the  required  system  responses  needed  for  system  chara- 
cterization from  an  uncountable  to  countable  number. 

The  first  scheme,  titled  the  piecewise  isoplanatic  approxi- 
mation (PIA)®”^,  requires  that  the  shape  of  the  system 
impulse  response  change  slowly  with  respect  to  a small 
interval  in  the  input  space.  Under  this  constraint,  the 
system  can  be  approximated  by  the  superposition  of  a number 
of  space-invariant  systems. 

The  second  constrained  system  characterization  in 
Chapter  II  is  termed  the  discrete  orthonormal  basis  set 
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response  characterization.  Here,  the  constraint  is  placed 

on  the  input  which  we  restrict  to  be  in  that  signal  class 

spanned  by  a given  orthonormal  basis  set.  In  this  case, 

the  system  is  completely  characterized  by  knowledge  of 

the  system  response  to  each  element  in  the  basis  set. 

Finally,  in  Chapter  II,  a number  of  seunpling  theorem 

characterizations  are  presented  that  correspond  to  various 

compact  support  and  bandlimited  constraints  on  the  system 

9-10 

input  and/or  line  spread  function.  For  the  most  part, 
these  characterizations  depend  on  both  the  input  and 
system. 

Chapter  III  presents  various  schemes  for  performing 
one-dimensional  space-variant  operations  utilizing  astig- 
matic processors. The  processors  here  are  as  general 
for  performing  space-variant  operations  as  that  in  Figure 
1-1  is  for  performing  space- invariant  operations.  We 
can,  however,  perform  only  one-dimensional  operations 
here  as  opposed  to  the  two-dimensional  capability  of  the 
processor  in  Figure  1-1.  The  basic  idea  for  the  one- 

12 

dimensional  processor  was  initially  proposed  by  Cutrona. 

Generalizations  and  implementations  of  1-D  processors  were 

13 

simultaneously  reported  by  Goodman  et.  al.  and  Marks 
14 

et.  al.  An  interesting  serendipity  is  a processor  for 
simultaneous  display  of  the  complete  2unbiguity  function 
display  first  suggested  by  Said  and  Cooper^^  and  later 
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reported  by  Marks  et.  al. 

Chapter  IV  discusses  the  present  state-of-the-art  in 
the  effort  to  formulate  a generalized  coherent  optical  pro- 
cessor for  performing  two-dimensional  space-variant  opera- 
tions. The  first  scheme  (the  composite  hologram  approach) 
is  a straightforward  implementation  of  a sampling  theorem 
developed  in  Chapter  II.  It  suffers,  however,  from  a 
strict  limitation  on  the  class  of  allowable  systems  and 
inputs . 

The  second  two-dimensional  space-variant  processor 

discussed  utilizes  the  selective  diffraction  efficiency 

effects  inherent  in  volvime  holography.  This  scheme  was 

19 

initially  proposed  by  Burton  et.  al.  and  was  developed 

via  the  PIA  by  Marks®  and  via  a saunpling  theorem  by  Been 
20-21 

et.  al.  The  volume  hologram  approach,  however,  suf- 

fers from  the  one-dimensional  nature  of  the  useful  angu- 
lar bandpass  which  prohibits  a straightforward  generali- 
zation of  the  scheme  to  two  dimensions. 

The  final  two-dimensional  space-variant  processing 

scheme  utilizes  phase  coded  reference  beams.  This  ap- 

22 

proach,  proposed  by  Krile  et.  al.,  smears  unwanted 
crosstalk  terms  into  diffuse  background  noise. 

Chapter  V contains  some  concluding  remarks  and  sug- 
gests some  possible  areas  of  future  investigation. 


CHAPTER  II 


^ 2.  GENERAL  LINEAR  SYSTEMS  CHARACTERIZATION  THEORY 

[ A system  may  be  abstracted  into  three  basic  components: 

E 

\ em  input,  a black  box,  and  an  output.  For  a given  set  of 

t . 

[ conditions  on  the  black  box,  the  output  is  simply  the  re- 

L sponse  of  the  black  box  to  a given  input.  A system  can  ^ 

[ be  multidimensional  and  have  associated  with  its  input- 

I I 

; output  domain  any  number  of  physical  quantities.  An 

i • 

[ V electronic  circuit-type  system,  for  example,  has  associ- 

ated with  its  input-output  domain  temporal  voltage  and 

r 

current  waveforms.  A coherent  processor,  on  the  other 
hand,  has  associated  with  its  input-output  domains  two- 
dimensional  spatial  variations  in  electric  field  intensity 
(light) . Input-output  domains  for  a given  system  need  not 
be  the  same.  For  example,  television  has  both  acoustical 
and  optical  output  domains  corresponding  to  a single  tem- 
poral voltage  input. 

In  a more  mathematical  sense,  we  can  think  of  a system  ^ 

in  terms  of  a mapping.  An  input  is  mapped  to  an  output, 
according  to  a deterministic  mapping  rule  dictated  by 

I . 

[ 

I the  black  box.  If  we  denote  our  input  by  f(y)  and  our 

f 

' output  by  g(x),  then^^ 

I 

t 

g(x)  - S[f(y)l.  (2-1) 


J 
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Here , S ( • ) corresponds  to  the  mapping  rule  or  system 
operator.  That  is,  the  output,  g(x),  corresponds  to  the 
input,  f(y),  being  operated  on  by  the  system  operator 
S(»).  Again,  the  input  and  output  can  correspond  to 
any  physical  quantity  such  as  voltage  or  temperature, 
and  the  variables  x and  y can  correspond  to  any  variable 
domain  (such  as  space  or  time)  in  any  desired  dimension. 
When  applied  to  coherent  optical  processors,  both  g(*) 
and  f(»)  denote  electric  field  amplitude,  and  both  x and 
y correspond  to  two-dimensional  spatial  output  and  input 
planes  respectively. 

The  general  system  characterization  in  Eq.  2-1  leaves 
little  room  for  further  development.  In  order  to  charac- 
terize this  system  completely,  we  would  need  to  catalog 
the  system  output  corresponding  to  every  possible  input 
under  all  possible  system  conditions.  To  avoid  performing 
of  this  fomidable  task,  further  ass\imptions  need  to  be 
made  concerning  out  black  box. 

2.1.  Classical  Linear  System  Characterization 

The  most  popular  system  assumption  is  that  of  lin- 
earity where  two  basic  suppositions  are  made:  1)  The 
response  of  the  system  to  the  sum  of  a nviraber  of  inputs 
is  equivalent  to  the  sum  of  the  system's  response  to  each 
of  the  inputs  individually  and  2)  amplification  of  the 


r 
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input  will  result  only  in  the  equivalent  amplification  of 
the  output.  These  system  properties  are  known  respectively 
as  superposition  and  homogeneity  and  together  constitute 
the  definition  of  linearity.  We  write  these  conditions 
respectively  as 

Sffj^Cy)  + f2(y)]  = S[fj^(y)]  + S[f2(y)]  (2-2) 

and 

S[af(y)]  = aS[f(y)]  (2-3) 

where  a is  a constant  independent  of  the  y variable. 

These  two  conditions  may  be  combined  into  a single  neces- 
sary and  sufficient  condition  for  linearity: 

S[afj^(y)  + bf2(y)]  = aS[fj^(y)]  + bS[f2(y)l. 

(2-4) 

While  no  physical  system  is  exactly  linear,  the  concept 
of  linearity  allows  powerful  mathematical  treatment  and 
in  many  instances  excellent  description  of  a system  pro- 
cess. 

Linear  systems  can  be  described  by  the  superposition 

integral  which  we  shall  now  derive.  We  begin  with  the 

23 

sifting  property  of  the  Dirac  delta  function: 

00 

f(y)  - / f(5)6(y-C)dC. 


(2-5) 
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This  relationship  is  substituted  into  the  system  charac- 
terization in  Eq.  2-1.  Since  the  system  operator,  S(»)» 
is  linear,  we  can  make  the  following  two  observations: 

1)  Thinking  of  the  integral  as  a continuum  sum,  we  can 
take  the  operator  inside  the  integral  sign  due  to  the 
superposition  property  of  Eq.  2-2.  2)  Since  the  system 
operator  operates  only  on  functions  of  y,  the  terms  f(C) 
and  dC  may  be  factored  out  of  the  system  operation  due  to 
the  homogeneity  property  of  Eq.  2-3.  This  leaves 

00 

g(x)  = / f(C)S[6(y-C)]dC.  (2-6) 

^00 

The  quantity  S[5(y-C)]  is  interpreted  as  the  system  re- 
sponse to  an  input  Dirac  delta  function  located  at  the 
input  coordinate  y = ?.  This  response  has  various  names, 
including  the  impulse  response  and  the  system  Greens 
function.  In  optics,  it  is  referred  to  as  either  the 
system  line  spread  function  or  the  system  point  spread 
function  depending  on  whether  analysis  is  being  performed 
in  one  or  two  dimensions.  Since  our  treatment  to  date 
has  a one- dimensional  flavor,  we  will  refer  to  this  re- 
sponse as  the  system  line  spread  function  and  write  it 
as 

h(x  - C;  OASIS  (y-  0 3.  (2-7) 

Then,  Eq.  2-6  becomes  the  superpositional  integral: 
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00 

g(x)  = / f(C)h(x  - C;C)dC,  (2-8) 

*00 

It  should  be  noted  that  the  use  of  the  line  spread 

function  notation  here  is  not  that  which  is  convention- 
23-24 

ally  used.  There  are,  however,  certain  computational 

e 25w27 

advantages  of  this  notation  ' which  will  soon  become 

clear. 

Let  us  now  physically  interpret  the  superposition 
integral  description  of  the  linear  system.  In  order  to 
completely  characterize  a system  in  this  manner,  we  must 
catalog  the  responses  of  our  system  to  Dirac  delta  inputs 
at  every  point  of  our  input  space.  The  number  of  re- 
quired responses  is  thus  uncountably  infinite. 

For  coherent  optical  systems,  we  make  the  following 
physical  interpretation:  A line  source  placed  on  our 
input  plane  at  y = possibly  by  the  focusing  of  an 
incident  plane  wave  with  a lens,  gives  rise  to  a field 
eunplitude  distribution  on  our  output  plane  which  is  the 
system  line  spread  function  corresponding  to  the  input 
coordinate  y = C.  As  the  input  line  source  explores  new 
positions  on  the  input  plane,  the  line  spread  function 
shifts  and/or  changes  shape. 

Suppose  that  the  line  spread  function  shifts  directly 
with  the  input  Dirac  delta  without  changing  shape.  That 
is,  if  we  move  our  input  Dirac  delta  a distance  d,  then 
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our  output  also  shifts  a distance  d and  otherwise  remains 
unchanged.  The  system  is  then  said  to  be  space-invariant 
and  the  line  spread  function  takes  on  the  form 

h(x  - C;  5)  h(x  - 5).  (2-9) 

In  temporal  systems,  this  property  is  referred  to  as  time- 
invariance  and  in  optics  as  spatial  invariance  or  isopla- 
nicity. 

Actually,  the  properties  of  linearity  and  invariance 
are  independent.  That  is,  there  are  systems  which  are 
invariant  and  not  linear.  An  invariant  system  is  simply 
one  in  which  the  output  shifts  directly  with  the  input  for 
all  inputs. 

For  linear  invariant  systems,  the  superposition  in- 
tegral becomes  the  convolution  integral: 

00 

g(x)  = / f(C)h(x  - C)dC.  (2-10) 

••00 

One  of  the  computational  advantages  of  linear  invariant 
systems  is  the  applicability  of  Fourier  analysis  in  their 
treatment.  We  define  the  Fourier  transform  of  a signal 
s(x)  as 

S(fx)  = F3j[s(x)] 

00 

= / s(x)exp(-j2iTf  x)dx.  (2-11) 

*00  ^ 

Here  Pjj(*)  refers  to  the  Fourier  transform  operator  with 


4, 
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respect  to  x,  and  f is  the  frequency  variable  correspond- 
ing  to  X.  Here,  and  henceforth,  a capital  letters  refers 
to  the  Fourier  transform  of  the  corresponding  lower  case 
letter,  (e.g.,  S(f^)  is  the  Fourier  transform  of  s(x)). 

Suppose,  then,  we  apply  the  Fourier  transform  oper- 
ator to  both  side  of  the  convolution  integral  in  Eq.  2-10. 
Since  the  Fourier  operator  is  linear,  we  may  bring  it 
through  the  integral  sign  and  past  the  input  to  operate 

only  on  the  system  line  spread  function.  Utilizing  the 
28 

shift  theorem  we  write 

G(f^)  = F^[g(x)] 

= r f (C)F^Ih(x)]exp(-j2Trf  C)dC 

»00 

This  familiar  relationship  is  the  statement  of  the  Four- 
ier domain  input-output  statement  of  a linear  invariant 
system.  The  Fourier  transform  of  the  line  spread  func- 
tion, referred  to  as  the  system  transfer  func- 

tion, or  simply  the  system  function.  The  relationship 
in  Eq.  2-12  reduces  the  convolution  integral  to  the  more 
computationally  attractive  operation  of  multiplication. 

The  system  discriptions  presented  to  this  point 

23—24  28—29 

are  relatively  standard.  ' In  the  remainder  of 

this  chapter,  various  other  methods  of  linear  system  char- 
acterization are  presented.  These  include  generalization  of 


the  superposition  integral  concept  to  other  response  map- 
pings of  continuum  orthonormal  basis  sets  of  which  the 
Dirac  delta  set  is  a special  case.  Also  presented  are 
three  schemes  by  which  the  required  number  of  cataloged 
responses  is  reduced  from  a continuum  to  a countably  in- 
finite set.  Each  of  these  three  methods  requires  assump- 
tions on  the  system  input  and/or  line  spread  function. 

As  will  be  seen,  the  three  schemes  are  highly  complemen- 
tary to  each  other.  The  overall  purpose  of  this  chapter 
is  to  establish  a wealth  of  linear  system  theory  to  be 
later  utilized  in  the  development  of  coherent  processors 
capable  of  performing  a wide  class  of  linear  operations. 
The  theory,  however,  is  applicable  to  all  linear  systems. 

2.2.  Continuum  Orthonormal  Response  Characterization 
The  conventional  superposition  integral  can  be 
viewed  as  a special  case  of  system  characterization 
through  cataloging  the  system  response  to  each  element 
in  a continuum  basis  set.  A continuum  set,  is 

said  to  be  orthonormal  if^® 

00 

/ = 5(x-y),  (2-13) 

*00  ^ ^ 

where  denotes  complex  conjugation.  In  our  linear 
system  treatment  in  the  previous  section,  {6(y  - C) } 
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obviously  conforms  to  this  definition. 

A system  input,  f (y) , can  be  expressed  in  terms  of 
a continuum  basis  set  by 

00 

f(y)  = J (2-14) 

where 

00 

a(C)  = / f (y)i)**(C)dy.  (2-15) 

«,00  ^ 

We  can  prove  the  validity  of  this  relationship  as  follows: 
Substitute  Eq.  2-15  into  Eq.  2-14  and  interchange  the 
order  of  integration.  Utilizing  the  orthonormality  con- 
dition of  Eq.  2-13  and  the  sifting  property  of  the  Dirac 
delta  function  in  Eq.  2-7  will  then  reduce  the  result  to 
the  identity  f (y)  = f (y) . 

To  develop  the  input-output  relationship  of  a linear 
system  corresponding  to  an  input  f(y),  we  substitute  Eq. 
2-14  into  our  system  operator  in  Eq.  2-1.  As  in  the 
development  of  the  superposition  integral,  we  make  the 
following  two  observations  concerning  our  linear  operator 
S(«);  1)  Viewing  the  integral  as  a continuum  sum,  we 

can  bring  the  system  operator  inside  the  integral  due 
to  the  superposition  property  of  Eq.  2-2.  2)  Since  the 

operator  only  operates  on  functions  of  y,  we  can  factor 
out  the  input  f(C)  and  the  differential  d?  due  to  the 


[, 


1 


17 


homogeneity  property  of  Eq.  2-3.  The  system  operator  thus 
only  operates  on  the  orthonormal  basis  set  element: 


g(x)  = / a(5)S[ijj  (C)]dC.  (2-16) 

*00  * 


Together  with  Eq.  2-15,  this  is  the  general  statement  of 
the  continuum  orthonormal  basis  set  response  characteri- 
zation of  a linear  system. 

Let's  interpret  it.  In  our  input  space,  we  place 
our  continuum  basis  set  element  4>y(5).  As  we  vary  C in 
'l'y(C)  the  corresponding  system  output,  S[if/y(C)]»  is  cat- 
aloged. With  knowledge  of  the  system  response  for  all 
the  system  is  completely  defined.  For  an  arbitrary  input, 
f (y) , the  functional  coefficient  a(C)  is  computed  from 
Eq.  2-15  and  substituted  into  Eq.  2-16  to  determine  the 
corresponding  system  output.  Note  that,  in  general,  the 
set  is  not  orthogonal. 

2.2.1.  Spatial  Frequency  Response 

As  was  mentioned,  the  superposition  integral  in  Eq. 
2-8  is  a special  case  of  this  relation.  Here,  'l'y(C)  = 

5 (y  - C) • By  the  sifting  property  of  the  Dirac  delta  in 
Eq.  2-5,  it  follows  that  a(C)  = f(C)  and,  from  Eq.  2-7, 
S[i|^y(C)l  = h(x  - C;  C)  . 

Another  continuum  basis  set  that  is  meaningful  in 
C engineering  applications  is  the  complex  sinusoid: 


1 


J. 


f 
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<l'y(v)  = exp(j2Tryv).  (2-17) 

Here,  v is  used  instead  of  C due  to  the  obvious  frequency 
personality  of  the  variable.  From  Eq.  2-15,  our  functional 
coefficient  is  seen  to  be  the  Fourier  transform  of  the 
input ; 

00 

a(v)  = / f (y)exp(-j2Trvy)dy  (2-18) 

••00 

= F(v). 

Due  to  its  extensive  use  later,  we  will  distinguish  the 
response  to  an  input  complex  sinusoid  with  the  name  "fre- 
quency response"  and  denote  it  by 


lc(x  - cv;  v)  A S [exp ( j2irvy)  ].  (2-19) 


This  particular  notation  is  used  due  to  certain  computa- 
tional advantages.^  The  constant  c merely  retains  dimens- 
ional consistancy  between  x and  v which  have  reciprocal 
units . 

Substituting  Eqs.  2-18  and  2-19  into  the  general 
continuum  orthonormal  response  characterization  of  Eq. 

2-16  gives 

00 

g(x)  = / F(v)k(x  - cv;  v)dv.  (2-20) 

*00 

Note  the  computational  similarities  between  this  frequency 

response  characterization  and  the  superposition  integral  ; 

J 

repeated  here: 


19 


00 

g(x)  = / f(C)h(x  - C;  Od?.  (2-21) 

••00 

From  our  definition  of  the  frequency  response  in  Eq.  2-19, 
we  can  use  the  superposition  integral  to  write 

00 

k(x  - cv;  u)  * /_^  h(x  -C;  5)exp(j2TiCv)dC. 

(2-22) 

As  such,  the  frequency  response  is  simply  the  inverse 

Fourier  transform  of  the  line  spread  function.  Inspection 

of  Eqs.  2-20  and  2-21  also  dictates  this  must  be  true  if 

28 

Parseval's  theorem  is  to  be  satisfied.  That  is,  if  f(C) 
and  F(v)  are  Fourier  transform  pairs,  then  h(x-C;C)  and 
k(x-cv;v)  must  also  be  Fourier  transform  pairs. 

2.2.2.  Output  Frequency  Mapping 

For  purposes  of  discussion,  assume  that  both  C and 
X are  spatial  variables.  It  follows  then  that  v and  f^ 
respectively  are  the  corresponding  frequency  variables. 
Thus,  the  superposition  integral  maps  an  input  spatial 
domain  to  an  output  spatial  domain  and  the  frequency  re- 
sponse characterization  maps  the  input  frequency  domain 
to  an  output  spatial  domain.  It  is  also  possible  to 
directly  map  to  the  output  frequency  domain.  To  do  this, 
in  accordance  to  our  observations  concerning  the  relation 
between  spatial  and  frequency  variables,  the  following 
two  Fourier  transform  operators  are  defined  on  function 
s(x,C)  ; 


J 
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F^[s(x;5)  ] A/_^  s(x;C)exp(-j2Trf^x)dx  (2-23) 

and 

00 

F^ts  (x;  C)  ] s(x;?)exp(-j2TTv5)dC.  (2-24) 

Application  of  the  linear  operator  Fjj(*)  to  both  sides  of 
the  superposition  integral  of  Eq.  (2-8)  gives 


G(fj^)  = / f(C)Fx[h(x  - C;  OldC 
»00 

= /Zoo  f(C)F^[h(x;  ?)]exp(-j2Trf^C)dC 


= F^  F^(f(5)h(x;  C)] 


(2-25) 


v=f 


where,  in  the  second  step  we  have  utilized  the  shift 

28 

theorem  of  Fourier  transform  theory.  Here,  we  are 
directly  mapping  the  input  spatial  domain  to  the  output 
frequency  domain.  Before  applying  a similar  procedure  to 
the  frequency  response  characterization,  it  is  instruc- 
tive to  discuss  the  interpretation  of  Eq.  2-25. 

First  of  all,  note  that  we  are  still  using  the  line 
spread  function  notation  given  in  Eq.  2-7.  The  computa- 
tional advantage  of  this  particular  notation  choice  is 
made  clear  in  the  development  of  Eq.  2-25.  Mathematic- 
ally, one  makes  the  transition  from  h(x  - C;  C)  to  h(x;  C) 
by  simple  variable  manipulation.  A straightforward  pro- 
cedure for  doing  this  is  to  make  the  substitution  x = (x-C) 
+ C for  all  x's  appearing  in  h(x  - C;  C) • To  form  h(x;  C) , 
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one  then  simply  substitutes  x's  for  all  (x  - C) 's. 

In  a more  physical  sense,  the  transition  from  h(x-C;C) 
to  h(x;  C)  can  probably  be  best  interpreted  by  visualizing 
a temporal  linear  system.  With  reference  to  Fig.  2.1a, 
we  first  excite  the  system  with  a Dirac  delta  at  input  time 
C ■ 0.  The  system  responds  in  this  example  with  the  tri- 
angular response  shown.  We  wait  a bit,  and  again  excite 
the  system  at  a later  time  C > 0.  The  system  response  is 
now  a triangle  of  shorter  duration.  Note  that  the  system 
shown  here  is  causal  in  that  we  can  have  no  system  re- 
sponse prior  to  application  of  an  input.  As  such,  the 
system  impulse  response  (line  spread  function)  must  be 
identically  zero  for  x - C < 0.  Thus  the  line  x - C = 0 
is  a reference  line  of  sorts  for  our  response.  In  the 
more  general  (non-causal)  case,  the  line  x •=  C is  a ref- 
erence or  a common  system  axis  for  comparison  of  the 
relation  of  the  location  of  the  input  Dirac  delta  to  the 
corresponding  system  line  spread  function.  As  is  shown 
in  Fig.  2.1b,  the  transition  from  h(x  - C;  C)  to  h(x; 
is  thus  a shift  of  this  reference  axis  to  the  origin  for 
all  possible  C* 

In  a cruder  sense,  we  can  envision  the  transition 
from  h(x  - C)  to  h(x;  as  follows.  For  the  invariant 
case,  h(x  - varies  directly  as  6 (y  - C) • That  is,  as 

6{y  - C)  moves  up  the  input  plane,  h(x  - C)  moves  up  the 
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output  plane  without  changing  shape.  In  the  linear  vari- 
ant case  we  can  also  envision  h(x  - C;  C)  moving  up  the 
plane/  but  changing  shape  as  it  does  so.  Returning  to 
the  invariant  case,  make  the  transition  from  h(x  - C)  to 
h(x)  as  in  Fig.  2.2.  As  the  input  Dirac  delta  moves  up 
the  plane,  h(x)  remains  unchanged  in  position  and  shape. 
Similarly,  h(x;  C)  may  be  thought  of  as  remaining  unchang- 
ing in  its  position.  That  is,  it  does  not  shift  as  before, 
but,  as  the  input  Dirac  delta  moves,  h(x;  C)  changes  shape. 
Again,  this  view  of  h(x  - 5)  vs  h(x;  C)  is  crude,  but 

does  provide  for  a somewhat  intuitive  feeling  of  physical 
interpretation . 

Returning  now  to  the  output  frequency  mapping  rela- 
tion in  Eq.  2-25,  we  can  generalize  the  idea  of  the  system 
transfer  function  to  a general  linear  system.  We  define 
the  system  transfer  function  as 

Hx(fx;  OAFjjthCx;  C)]-  (2-26) 

Our  output  frequency  mapping  in  Eq.  2-25  can  then  be 
written 

G(f^)  =F.[f (?)H^(f  ; 5)]  (2-27) 

X XX  ^ 

Note  that  the  generalized  transfer  function  is  sensitive 
to  input  location,  C*  Note  also  that,  for  the  invariant 
case,  takes  on  the  famil- 
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iar  product  form  given  in  Eq.  2-12. 

Equation  2-25  gives  an  input  spatial  to  output  fre- 
quency mapping.  An  input  frequency  to  output  frequency 
mapping  can  be  developed  by  applying  similar  steps  to 
the  frequency  response  characterization  in  Eq.  2-20: 


G(f„)  = / F(v)  F tk(x  - cv;  v) ]dv 


= /“  F(v)  F^[k(x;  v)  lexp(-j2Trcvf^)dv 


5 = -cf. 


(2-28) 

where  the  inverse  Fourier  transform  operator  with  respect 
to  5 is  given  by 

F"^[S(x,v)]  = S(x;  v) exp( j2nvC)dv.  (2-29) 


From  Eqs.  2-20,  2-21,  2-25,  and  2-28,  we  have  expres- 
sions covering  all  possible  frequency/spatial- input/out- 
put mappings.  A sxammary  is  given  in  Table  2.1.  These 
relationships,  which  have  roots  in  the  generalized  con- 
tinuum orthonormal  basis  set  response  characterization 
of  linear  systems,  lend  a flexibility  to  linear  system 
representation  which  will  later  be  utilized  in  linear 
coherent  processor  implementation  schemes. 

2.3.  The  Piecewise  Isoplanatic  Approximation  (PIA) 

In  the  previous  section,  we  require  a continuum 
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DOMAIN 
input  — » output 
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spatial 

spatial 

0(x)»r  f(() 
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frequency 

spatial 

/»• 

g(x)*/  F(i»)  k(x“C»i»')dy 

2-20 

spatial 

frequency 

2-25 

frequency 

frequency 

« • L J '{.-cf. 

2-28 
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Teible  2.1  Summary  of  all  possible  input  to  output 
mappings  from  space  and/or  frequency 
domain  to  space  and/or  frequency  domain 
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cataloging  of  system  responses  in  order  to  completely  char- 
acterize our  linear  system.  By  making  certain  assumptions 
or  approximations  concerning  the  system  and/or  input  class, 
it  is  possible  to  reduce  this  to  a countably  infinite  num- 
ber of  required  input-output  relationships.  We  will  dis- 
cuss three  such  linear  system  characterizations  the  first 
of  which  is  the  piecewise  isoplanatic  approximation  (PIA) J 
In  optics  nomenclature,  the  terms  space-invariant  and 
isoplanatic  are  equivalent.  Thus,  an  isoplanatic  system's 
input-output  relationship  is  given  by  the  convolution  in- 
tegral. A class  of  systems  which  lie  between  the  isopla- 
natic and  space-variant  classification  consists  of  those 
systems  which  are  piecewise  isoplanatic.  The  input 

spaces  of  such  systems  consist  of  a number  of  adjoined 

31-32 

isoplanatic  patches.  For  this  class  of  linear  systems, 

as  an  input  Dirac  delta  explores  a given  isoplanatic  patch, 
the  corresponding  line  spread  function  shifts  directly 
without  changing  shape.  That  is,  within  a single  isopla- 
natic patch,  the  system  is  space-invariant.  When  the  in- 
put Dirac  delta  moves  from  one  ^atch  to  the  next,  the  line 
spread  function  immediately  changes  shape.  While  the  in- 
put Dirac  delta  explores  this  second  isoplanatic  patch, 
however,  the  corresponding  output  response  will  again  shift 
accordingly.  Thus,  we  may  interpret  the  piecewise  iso- 
planatic system  as  a number  of  space-invariant  systems. 
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A given  input,  f(5),  is  divided  into  isoplanatic  regions. 
Each  isoplanatic  input  region  serves  as  an  input  to  a space- 
invariant  system  defined  by  the  impulse  response  assigned 
to  the  isoplanatic  patch.  The  outputs  of  each  of  these 
component  systems  are  summed  to  give  the  output  of  the 
overall  piecewise  isoplanatic  system. 

The  power  of  the  piecewise  isoplanatic  system  con- 
cept lies  in  its  use  for  approximating  space-variant 
systems.  If  the  line  spread  function  of  a space-variant 
system  does  not  change  shape  drastically  corresponding 
to  small  input  region,  then  we  may  approximate  the  system 
as  isoplanatic  over  this  region.  The  result  of  making 
this  approximation  over  a large  number  of  disjoint  regions 
is  appropriately  termed  a piecewise  isoplanatic  approxi- 

7 

mation  or  PIA.  An  illustration  of  a possible  PIA  as 
applied  to  an  example  line  spread  function  is  pictured 
in  Fig.  2.3. 

We  now  will  proceed  with  a brief  overview  of  the 
mathematics  of  the  PIA.  Let  the  n^^  input  isoplanatic 
patch  in  our  one-dimensional  treatment  extend  from 
to  The  midpoint  of  this  patch  is 

(2-30) 

and  the  patch's  width  is 


fe 


Figure  2.3  A possible  piecewise  isoplanatic  approximation 
(PIA)  of  the  space-variant  system  line  spread 
f\inction  in  Figure  2.1. 


I 


. 

n n n 


(2-31) 


We  assume  adjacent  patches  are  adjoined  so  that  the  upper 
point  on  the  n^^  patch  is  equivalent  to  the  lower  point 
of  the  (n  + 1)®^  patch.  That  is 


'^n  *'n  + 1, 


(2-32) 


For  each  isoplanatically  modeled  patch,  we  must  some- 
how choose  a single  space-invariant  line  spread  function. 

If  our  parent  space-variant  system  has  line  spread  function 
h(x  - C) , it  follows  that  a reasonable  choice  of  the 
isoplanatic  patch  line  spread  function  would  be  h(x-Cj^;Cjj) 
where 


S,  < C < V 
n — n — n* 


(2-33) 


This  is,  our  descriptive  space- invariant  line  spread  func- 
tion should  correspond  to  the  system  response  to  an  input 
Dirac  delta  function  placed  somewhere  within  the  isopla- 
natic patch. 

Suppose  then,  we  have  knowledge  of 


h(x  - 5^;  5^)  = S[6(y  - Cjj)]  (2-34) 


for  each  isoplanatic  patch.  We  must  now  divide  our  in- 
put, f(C)f  into  isoplanatic  patches.  Thus,  let 


f(5)  = I f„(E  - 5„) 


(2-35) 


;H|  .Bt] 


nnT 
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where 

C-m 

frj(C-Cn)4f  (C)rect(-^).  (2-36) 

n 

The  output  of  our  piecewise  isoplanatically  modeled 
system  can  now  be  determined  by  operating  on  the  input 
with  our  linear  system  operator  in  Eg.  2-1.  We  can  ex- 
tract the  summation  sign  due  to  the  homogeneity  property 
of  Eg.  2.2  to  give 

g(x)  = I S[f^(y  - ?„)].  (2-37) 

n 

VJe  place  the  over  g to  denote  that  our  computed  out- 
put is  a PIA  of  the  true  output.  Over  each  segment  of 
the  input  lying  within  an  isoplanatic  patch,  the  system 
ope'rator  is  space-invariant  and  can  thus  be  expressed 
via  the  convolution  integral  in  Eg.  2-10.  Our  desired 
line  spread  function  is  given  by  Eg.  2-34.  Thus; 

S[f  (5  - 5„)I  = / f„(5  - 5n)h(x  - ?!  5^)d{ 

= fjj(x  - c^)*h(x;  ?^^)  (2-38) 

where  denotes  convolution.  Substituting  this  result 
into  Eg.  2-37  gives 

g(x)  = I f^(x  - ?^)*h(x;  K^),  (2-39) 

n 

Using  the  shift  theorem  and  the  linearity  of  the  Fourier 


i 

t 


i 
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operator,^®  we  can  write  the  Fourier  transform  of  the  PIA 
as 

°<'x’  * I 

(2-40) 

We  now  inverse  transform,  this  time  grouping  the  complex 
exponential  with  the  system  transfer  function.  Doing 
this  changes  the  shift  from  the  input  to  the  line  spread 
function: 


g(x)  = I f^(x)*h^(x  - (2-41) 

n 

This  is  the  desired  final  form  of  the  PIA.  Note  that  we 
have  here  approximated  our  superposition  integral  as  a 
summation  of  convolutions.  Equivalently,  we  can  interpret 
the  PIA  as  a number  of  outputs  of  isoplanatic  systems 
which,  when  superimposed,  give  an  approximate  output  to 
the  parent  space-variant  system. 

To  further  compare  the  PIA  output  to  the  true  out- 
put, we  substitute  Eq.  2-36  into  Eq.  2-41  and  express 
our  convolutional  operation  in  integral  form.  The  re- 
sult is 


g(x)  = I f(C)h(x  - C;  5-)dC. 

^ "n 


(2-42) 


J 


\ 


Consider  then,  the  system  output  given  by  the  superposi- 
tion integral  in  Eq.  2-21.  We  can  break  up  our  range  of 


33 


integration  from  -»  to  “ into  those  regions  we  have  modeled 
as  isoplanatic  and  write: 

u 

g(x)  = I J ^ f(5)h(x  - Od?.  (2-43) 

” ^n 

One  sees  that  as  each  isoplanatic  patch  width  narrows 
around  (matched  by  a corresponding  increase  in  patch 
density) , the  PIA  in  Eg.  2-42  approaches  the  true  output 
given  by  the  superposition  integral  in  Eq.  2-43. 

Let  us  now  discuss  some  of  the  properties  of  the  PIA. 
First  of  all,  as  promised,  the  number  of  required  system 
responses  has  been  reduced  from  the  continuum  number  re- 
quired by  the  superposition  integral  to  a countable  number. 
That  is,  for  the  PIA,  we  need  only  one  defining  system  re- 
lationship per  isoplanatic  patch.  For  inputs  with  finite 
support  (space  limited  inputs) , the  PIA  will  require  only 
a finite  nvimber  of  defining  relationships.  The  price  we 
pay  lies  in  the  work  "approximation."  How  good  an  approxi- 
mation is  the  PIA?  For  the  general  case,  this  is  a dif- 
ficult question  to  answer.  The  development  and  use  of 
the  PIA  is  simply  intuitively  pleasing.  For  some  specific 
example  systems,  the  reader  is  referred  to  the  work  of 
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Marks  and  Krile  who  have  shown  that  the  PIA  gives  "good" 
results  when  applied  to  the  magnifier  and  Fourier  trans- 
former. It  should  be  noted  also  that  when  applied  to 
invariant  systems,  the  PIA  and  true  outputs  are  identical 


since  invariant  systems  can  be  thought  of  as  a special 
class  of  piecewise  isoplanatic  systems. 

An  important  aspect  of  the  PIA  scheme  is  its  charac- 
teristic of  being  essentially  input  independent.  That 
is,  the  manner  in  which  our  line  spread  function  changes 
shape  as  an  input  Dirac  delta  shifts  is  determined  solely 
by  the  system.  If  the  line  spread  function  changes  shape 
only  slightly  corresponding  to  an  input  interval,  then 
we  can  model  that  interval  as  an  isoplanatic  patch.  The 
act  of  input  plane  patch  calibration  is  therefore  inde- 
pendent of  the  input  we  choose. 

As  a final  note  in  this  section,  we  should  note  that 
the  PIA  can  also  be  applied  to  the  input  frequency  domain 
using  the  system  frequency  response  instead  of  the  line 
spread  function.  Although  the  concept  of  a "piecewise 
frequency-invariant*”  system  requires  a different  physical 
interpretation,  the  mathematical  development  is  identical 
to  that  of  the  PIA.  One  merely  divides  up  the  input 
spectrum  into  frequency-invariant  (instead  of  isoplanatic) 
patches  and  utilizes  a sample  frequency  response  (instead 
of  line  spread  function)  from  somewhere  within  each  patch. 


♦Frequency-invariant  systems  are  defined  rigourously  in 
Appendix  A. 
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2.4.  Discrete  Orthonormal  Basis  Set  Response  (DBR)  Charac- 
terization 

2.4.1.  General 

The  piecewise  isoplanatic  approximation  is  a method 
by  which  the  required  number  of  defining  input-output  re- 
lationships of  a system  can  be  reduced  to  a countable  set. 
Another  method  of  performing  such  a reduction  uses  a method 
somewhat  similar  to  the  previously  discussed  continuum 
orthonormal  basis  set  response  characterization.  The  only 
difference  is  that  the  orthonormal  basis  set  is  now  dis- 
crete. Whereas  the  PIA  is  primarily  system  dependent,  the 

O 

discrete  orthonormal  basis  set  response  characterization 
(DBR  for  short)  will  be  seen  to  be  primarily  dependent 
on  the  system  input. 

In  order  to  reduce  our  required  number  of  defining 
input-output  relations  to  a countable  set,  it  is  always 
necessary  to  make  a restrictive  assumption.  In  the  PIA, 
we  assumed  the  space-variant  system  was  piecewise  isopla- 
natic. For  the  DBR  characterization,  our  restrictive 
assiimption  is  that  the  class  of  allowable  inputs  can  be 
expressed  in  terms  of  a specified  orthonormal  basis  set. 

As  we  will  see,  the  system  can  then  be  completely  charac- 
terized by  knowledge  of  the  system  response  to  each  ele- 
ment in  the  given  discrete  orthonormal  basis  set. 


y 
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Before  developing  the  DBR  characterization  of  a linear 
system,  it  is  instructive  to  review  the  concepts  of  the 
discrete  orthonormal  basis  set.  A set  of  functions, 

n 

is  said  to  be  orthonormal  in  the  discrete  sense,  if 


/ ♦„(y)%(y)dy  = 

wOO 


(2-44) 


where  6 is  the  Kronecker  delta: 
nm 


1 ; n = m 

0 ; n m . 


(2-45) 


A discrete  orthonormal  basis  set  is  said  to  be  complete' 


I 'l'jj(C)'l'j^(ri)  =6(5-  n). 
n 


(2-46) 


A signal,  f(y),  can  be  expressed  in  terms  of  a complete 
discrete  orthonormal  basis  set  by  the  expansion 


*1  “n’^'n^^^ 
n 


(2-47) 


where  the  coefficients,  a^,  are  found  from  the  inner  pro- 
duct 


cijj  * / f (y)'P„(y)dy. 


(2-48) 


We  can  prove  the  validity  of  this  expansion  as  fol- 
lows. First,  substitute  Eq.  2-48  into  Eq.  2-47  and  ex- 
change orders  of  summation  and  integration.  Utilizing 
the  completeness  property  in  Eq.  2-46  and  the  sifting 


property  of  the  Dirac  delta  in  Eq.  2-5  then  reduces  the 
relationship  to  the  identity  f (y)  = f (y) . 
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The  DBR  characterization  of  a linear  system  is  found 
by  substituting  the  orthonormal  expansion  in  Eq.  2-47 
into  the  linear  system  operator  of  Eq.  2-1.  We  can  extract 
the  summation  sign  from  the  system  operator  due  to  the 
superposition  property  in  Eq.  2-2.  The  coefficients  may 
be  factored  out  by  using  the  homogeneity  property  of  Eq. 
2-3.  The  result  is 

9(x)  » ^ * (2-48) 

n 

Thus,  we  can  express  the  output  to  the  linear  system  with 
knowledge  of  the  system's  response  to  each  element  in  the 
discrete  orthonormal  basis  set.  Let  us  denote  this  re- 
sponse by 

« S[i))^(y)].  (2-49) 

Thus,  our  expansion  for  the  system  output  in  Eq.  2-48 
becomes 

g(x)  » (2-50) 

n 

We  can  relate  the  response  set  elements  in  terms  of  the 
system  line  spread  function  by 

00 

0n(x)  * / 

••00 
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Note  that  the  DBR  set,  {i}>^(x)},  is  in  general  not  ortho- 
gonal . 

2.4.2.  The  Sine  and  Rect  Response  Characterizations 

With  the  general  concept  of  the  DBR  characterization 

of  linear  systems  established,  we  now  specifically  address 

two  classes  of  inputs  commonly  encountered  in  signal  pro- 
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cessing.  The  first  is  the  bandlimited  signal  whose 
expansion  coefficients  are  simply  the  signal's  sampled 
values  for  an  appropriately  chosen  basis  set.  The  second 
class  are  those  inputs  with  finite  support  which  can  be 
expressed  in  the  familiar  Fourier  series  expansion. 

A bandlimited  signal  is,  by  definition,  a signal 
whose  Fourier  transform  has  finite  support.  We  will  here 
specifically  consider  that  class  of  low-pass  signals, 
f(C),  whose  Fourier  transform  is  identically  zero  outside 
the  interval  |v|  £ W.  As  such,  we  can  write 

F(v)  = F(v)rect(^).  (2-52) 

34 

Such  signals  can  be  expressed  via  the  Whittaker  - Shan- 
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non  sampling  theorem  : 

f(C)  * I f(5j^)sinc  2W(C  - (2-53) 

n 

where 


(2-54) 
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and 

sine  X A . (2-55) 

The  uniformly  converging^®  seunpling  theorem  expansion  in 
Eq.  2-54  can  be  considered  as  a special  case  of  a dis- 
crete orthonormal  expansion.  Here/  our  basis  set  ele- 
ments are 

- /2W  sine  2W(y  - (2-56) 

The  expansion  coefficients  are  simply  the  signal's  sample 
values : 

00 

On  - /2W  / f(C)sinc  2W(5  - 5^^)  d? . (2-57) 

wOO 

Utilizing  Parseval's  theorem  gives 

-1  • 

otn  * /2w  / F(v)exp(j2TrvC^)rect(^)dv 
^00 

^ CO 

- / F(v)exp(j2TivC  )dv  (2-58) 

«00 

_1 

= /3w  f(C„). 

To  characterize  a linear  system  for  the  class  of 
bandlimited  inputs/  we  need  knowledge  of  the  system  re- 
sponse to  all  elements  in  the  sine  basis  set  of  Eq.  2-56. 
In  accordance  with  Eqs.  2-50  and  2-58/  our  system  output 
can  be  written  as 
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g(x)  = I f(C^)  S[sinc  2W(y  - 5^)].  (2-59) 
n 

This  "sine-response"  characterization  has  quite  nice  im- 
plementation properties  as  far  as  coherent  processing  is 
concerned.  First  of  all,  the  sine  function  is  one  that 
can  be  generated  by  quite  straightforward  optical  means. 

We  will  dwell  on  this  facet  further  in  Chapter  IV.  Secondly, 
the  expansion  coefficients  are  simply  the  sample  values 
of  the  signal.  As  such,  no  inner  product  computation, 
such  as  in  Eq.  2-48,  is  required.  Again  the  sine  response 
linear  system  characterization  is  valid  only  for  the  class 
of  bandlimited  inputs. 

A second  input  class  consists  of  those  inputs  of  fin- 
ite support.  This  class  may  be  considered  the  Fourier  dual 
of  the  bandlimited  signal  whose  Fourier  transform  has  fin- 
ite support.  In  the  bandlimited  case,  sampling  is  per- 
formed in  the  spatial  domain.  For  the  class  to  now  be 
considered,  seunpling  will  be  performed  in  the  frequency 
domain. 

We  here  consider  only  those  inputs  that  are  identically 
zero  outside  the  interval  |c|  IT.  That  is 

f(C)  = f(C)rect(|^).  (2-60) 

O A 

Such  signals  can  be  expressed  in  a Fourier  series^  expan- 
sion: 
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fiv  ^ I P(f^)exp(j2iTf^C)rect(|^)  (2-61) 

n 

where  = n/2T.  The  rect  term  merely  retains  the  zeroth 

order  term  of  the  periodic  expansion  which  is  our  input. 

The  Fourier  series  is  also  a special  case  of  a discrete 

orthonormal  expansion.  The  n^^  orthonormal  basis  function 

for  the  Fourier  series  is 

_1 

4»„(y)  » /5t  e^^’^^n^  rect(^).  (2-62) 

The  expansion  coefficients  are  thus  seunple  values  of  the 
input's  Fourier  transform; 

J.  » 

Un  = /5t  / f (C)exp(-j2Trf^C)dC 
_1 

= /3f  F(f„) . (2-63) 

n 

In  this  computation,  we  have  included  the  finite  support 
of  the  input  by  virtue  of  Eq.  2-60. 

To  characterize  a linear  system  for  the  specified 
class  of  finite  support  (or  space-limited)  inputs,  we 
require  knowledge  of  the  system  response  to  all  elements 
of  the  form  of  Eq.  2-62.  With  reference  to  Eqs.  2-50  j 

i 

and  2-63,  our  system  output  is 

9(x)  * 5^  E ^^^n^  S[exp( j2iTf^y)rect(^)  ] . 
n 

(2-64) 
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As  with  the  sine-response  characterization,  this  "rect- 
response"  characterization  has  quite  'nice  implementation 
properties  when  applied  to  coherent  processing.  The  re- 
quired orthonormal  basis  set  elements  can  be  generated 
straightforwardly  and  the  expansion  coefficients  are  simply 
the  sample  values  of  the  input's  Fourier  transform. 

The  sine  and  rect  responses  are  only  two  examples  of 
orthonormal  basis  sets  that  can  be  applied  to  the  charac- 
terization of  a linear  system.  They  do,  however,  utilize 
elements  that,  as  we  shall  see,  are  easily  generated  by 
coherent  optical  techniques.  Note  also  that,  in  general, 

DBR  characterizations  of  linear  systems  places  restriction 
only  on  the  permissable  class  of  inputs  and  is  thus  system 
independent . 

2.5.  Sampling  Theorem  Representations 

In  the  previous  section,  we  have  seen  that  the  class 
of  bandlimited  signals  can  be  expressed  via  the  Whittaker- 
Shannon  sampling  theorem.  For  certain  systems,  the  sampl- 
ing theorem  can  also  be  applied  to  the  line  spread  func- 
tion. We  will  here  look  at  two  cases.  The  first  is  when 
the  line  spread  function,  h(x;  ^) , is  bandlimited  in  x. 

In  this  case,  we  can  express  the  system  output  an  expansion 
utilizing  scunples  of  the  line  spread  function  of  the  form 
h(x^;  C) • This  scheme  is  essentially  input  independent 
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and  thus  can  be  grouped  »/ith  the  PIA  in  this  regard.  The 
second  case  where  the  sunpling  theorem  is  applicable  to 
system  characterizatioi.  is  when  the  line  spread  function, 
h(x;  C) / is  bandlimited  in  In  order  to  apply  the 
sampling  theorem  to  this  class  of  systems,  the  system  in- 
put must  also  be  bandlimited.  Thus,  unlike  any  other 
method  thus  far  discussed,  this  scheme  is  both  input  and 
system  dependent. 

The  class  of  linear  systems  to  which  the  two  scimpling 
theorems  are  applicable  can  be  enlarged  by  considering  the 
system  frequency  response  instead  of  the  line  spread  func- 
tion. For  exeunple,  if  h(x;  ?)  is  not  bandlimited  in  x or 
C,  it  might  be  possible  that  the  corresponding  k(x  - cv;  v) 
might  be  bandlimited  with  respect  to  one  of  its  variable. 

For  both  of  the  sampling  theorems  considered,  the 
system  output  is  shown  to  be  bandlimited  and,  as  such, 
can  be  expressed  by  a conventional  Whittaker-Shannon 
Seunpling  theorem  expansion.  The  required  sample  values 
of  the  output  are  shown  to  be  the  product  of  infinite 
matrices  the  components  of  which  are  the  sample  values 
of  the  system  line  spread  function  and  input.  A great 
simplification  occurs  for  those  systems  falling  into  both 
classes  for  which  the  sampling  theorems  are  applicable, 
that  is,  for  line  spread  functions,  h(x;  C) / that  are 
bandlimited  in  both  x and 
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2.5.1.  Line  Spread  Function  Bandlimited  in  x 

When  we  say  a line  spread  function,  h(x;  C),  is  band- 
limited  in  X,  we  simply  mean  that  the  generalized  transfer 
fxinction,  H^(f^;  5)  as  defined  in  Eq.  2-26,  has  finite 
support  in  f^  for  all  Let  us  consider  only  the  low 
pass  case  where  H^(f^;  C)  is  identically  zero  for  all 
! f 1 > W . That  is 

X X 

«x^^x'  = «x<^x'  1.  (2-65) 

This  relation  is  satisfied  if  our  output  line  spread 
function  is  bandlimited  irrespective  of  the  location  of 
our  input  Dirac  delta.  V7e  would  thus  expect  that,  regard- 
less of  our  input,  our  output  would  be  bandlimited.  We 
can  show  this  by  substituting  Eq.  2-65  into  the  output 
frequency  relationship  in  Eq.  2-27: 


rect(^).  (2-66) 

V = f^  X 

Thus,  irrespective  of  our  input,  the  corresponding  spectrum 
(Fourier  transform)  of  the  output  will  be  identically  zero 
for  all  jf^l  > As  such,  we  can  express  the  system  out- 

put by  the  following  sampling  theorem  expansion: 

g(x)  = I g(Xj^)sinc  2Wj^  (x  - x^)  . (2-67) 

n 

The  output  sample  coefficients  of  this  expansion  are  given 
by 


G(f^)  =F^[f(?)H^(f^;  C)] 


g(Xj^)  = / f(C)h(Xj^  - C;  C)dC.(2-68) 
«00 

We  will  later  give  alternate  more  computationally  prone 

expressions  of  the  sample  output.  The  system  line  spread 

2 

function  can  also  be  expressed  in  sampling  theorem  form: 


h(x;  0=1  “ ^n^  • 

" (2-69) 


Substituting  this  result  into  the  output  frequency  rela- 
tionship in  Eq.  2-25  gives 


G(fx)  = F^F^[f(^)  I h(x^;  ?)sinc  2W^(x  - x^) ] 


V » f 


= 2W^  I F^[f(C)h(x^;  Ol 


n 


exp  ( - j 2Trxf  ) rect 

V. 


(2-70) 


where,  in  the  second  step,  we  have  recognized  that 

1 ^x 

F^fsinc  2W^(x-x^)]  = exp  {-j2Trf^Xj^)  rect  (^^). 

X X 

(2-71) 

Utilizing  the  similarity  in  definition  between  F^(*)  and 
F^(»)  as  given  in  Eqs.  2-23  and  2-24,  we  can  rewrite  Eq. 
2-70  as 

G(fx)  = ^ I Fj^tf  (x)h(x^;  x)  ] exp  (-j2irXj^f  j^). 

X n 

(2-72) 


Inverse  Fourier  transformation  then  gives 
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g(x)  = I [f(x)h(Xj^;  x)]*sinc  2W^(x  - x^) 

" (2-73) 

where,  again,  denotes  convolution.  Note  that  in  this 

expression  we  have  expressed  the  system  output  as  a sum- 
mation of  convolutions.  Thus,  in  effect,  the  space-variant 
system  is  represented  by  a number  of  isoplanatic  systems 
each  with  an  input-output  relationship  dictated  by  a con- 
volution operation. 

In  our  seunpling  theorem  expression  in  Eq.  2-73  we 
are  using  seunple  line  spread  functions  of  the  form  h(x^;C) 
as  exemplified  in  Figure  2-4.  We  physically  measure  these 
Seunple  functions  by  placing  on  our  output  plane  an  array 
of  detectors  spaced  as  x^  = n/2W^.  These  detectors  meas- 
ure the  system  output  for  an  input  6 (y) . As  the  input 
Dirac  delta  moves  to  <5(y  - C),  our  corresponding  output 
is  h(x  - C;  C) . In  order  for  the  n^^  detector  to  meas- 
ure we  must  also  shift  our  entire  output  array 

to  be  centered  at  x = C-  The  output  measured  by  the  n 
detector  as  a function  of  ^ is  the  desired  K)  • 

It  is  also  possible,  however,  to  measure  our  seunple 
line  spread  function  using  an  array  that  does  not  shift 
in  accordance  with  the  location  of  our  input  Dirac  delta. 
This  follows  from  the  observation  that  if  h(x;  5)  is 

handlimited  in  x,  then  so  is  h(x  - C;  C) • That  is  from 

28 

the  shift  theorem 


i 
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Fj^[h(x  - C;  C)]  » F^[h(x;  C)  ] exp (- j2TrCf . (2-74) 

Our  measurement  scheme  here  utilizes  the  same  detector 
array  as  before  but  with  no  shifting  of  the  detectors. 

As  the  input  Dirac  delta  is  positioned  at  different  values 
of  the  n^^  output  detector  will  measure  h(x^  - C). 

Note  that  through  use  of  this  sampling  scheme,  we  cannot 
regain  the  functional  form  h(x^;  5). 

2.5.2.  Line  Spread  Function  Bandlimited  in  C 

In  the  previous  section,  we  described  a sampling 
theorem  for  linear  space-variant  systems  with  line  spread 
functions  that  are  bandlimited  in  x.  That  is,  the  system 
transfer  function,  C)»  has  finite  support  in  f^  for 

all  C*  An  alternate  sampling  theorem  exists  when  the 
line  spread  function  h(x;  5) » is  bandlimited  in  C for  all 

9 

X.  To  facilitate  definition  of  the  corresponding  band- 
width, we  define  the  system  "variation  spectrum"  as 

H^(x;  V)  = F^[h(x;  C)].  (2-74) 

Note  the  similarity  between  the  definitions  of  the  system 
variation  spectrum  and  the  system  transfer  function  repeated 
here  from  Eq.  2-26: 

Hx(fx;  C)  = F^th(x;  C) 1 . (2-75) 


J 
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The  function,  is  called  the  variation  spec- 

trum because  its  support  is  a measure  of  how  the  system 
line  spread  function  changes  shape  when  C is  varied.  The 
system  is  said  to  be  variation  limited  if  H^(x;  v)  has 
finite  support  in  v for  all  x.  We  will  here  consider  the 
low  pass  case  where 

H^(x;  v)  = H^(x;  v)rect(2;^).  (2-76) 

The  quantity  2W^  is  appropriately  called  the  variation 
bandwidth. 

To  illustrate  the  application  of  the  variation  band- 
width to  measure  the  spatial  variance  of  a system,  con- 

• i 

sider  the  case  where  the  system  is  isoplanatic.  That  is  1 

h(x;  C)  -►  h(x).  In  this  case,  the  corresponding  varia-  j 

tion  spectrum  is  H^(x;  v)  = h(x)6(v).  The  corresponding 
variation  bandwidth  is  zero.  Thus,  we  can  conclude,  as 
expected,  that  an  isoplanatic  system  is  indeed  space- in- 
variant. 

Returning  now  to  the  sampling  theorem  derivation, 
we  inspect  the  computational  form  of  the  output  frequency 

i 

relationship  in  Eq.  2-25  which  is  repeated  here:  j 
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bandlimited  in  C is  that  our  input,  f(C)/  be  bandlimited 
and  that  our  system  be  variation  limited  with  a variation 


bandwidth  of,  say,  2W  . Let  the  bandwidth  of  our  input  be 


2W^  where 


F(v)  = F^[f(C)] 


= F(v)  rect  (2^) . 


(2-78) 


Multiplication  in  the  spatial  (^)  domain  corresponds  to 
convolution  in  the  frequency  (v)  domain.  If  two  functions 
have  finite  support,  then  the  support  of  the  convolution 
of  the  two  functions  is  upper  bounded  by  the  sum  of  the 
supports  of  the  component  functions.  Thus,  we  conclude 
that  the  upper  bound  for  the  bandwidth  of  the  product 
f(C)h(x;C)  in  C is  given  by  the  sum  of  the  component  band- 
widths.  That  is 


F^[f(C)h(x;  C)]  = F^[f(C)h(x;  C)]rect(^)  (2-79) 

s 


where  the  bandwidth  sum,  2W  , is  given  by 

S 


2W„  « 2W-  + 2W,. 
s r V 


(2-80) 


We  can  immediately  express  the  product  f(C)h(x;C)  in  a 
S2unpling  theorem  expansion: 


f(5)h(x;  O • I Cn)sinc  2Wg(5  - 


n 


(2-81) 


J 
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g(x)  ^ I f(5n)h{x  - Cj,)*sinc  2WgX . (2-84) 

n 

As  with  the  case  where  h(x;  C)  was  bandlimited  in  x, 
we  have  here  reduced  our  system  characterization  to  a 
summation  of  convolutions  or,  equivalently,  the  output 
of  our  space-variant  system  is  given  by  the  superposition 
of  the  outputs  of  a number  of  isoplanatic  systems.  In 
Eq.  2-84,  however,  the  input  as  well  as  the  line  spread 
function  is  sampled  in  contrast  to  the  previous  scheme 
of  Eq.  2-72  where  only  the  line  spread  function  is  sampled. 
Thus,  this  scheme  is  both  input  and  system  dependent. 

In  a sense  we  have  oversaunpled  in  the  expression  in  Eq. 
2-84.  That  is,  the  minimum  allowable  sampling  rate  for 


our  input  is  2W^  and  the  line  spread  function  requires 
a saunpling  rate  of  only  2W^.  But  we  have  sampled  both  the 
input  and  line  spread  function  at  a rate  of  = 2Wj  + 

2W^.  As  is  shown  in  Appendix  B,  however,  sampling  at  the 
minimum  allowable  rates  results  in  an  expression  that  is 
highly  unattractive  in  a computational  sense  when  compared 
to  the  intuitively  pleasing  form  of  Eq.  2-84. 

The  physical  interpretation  of  sampling  h(x;  5)  in  ^ 
is  much  different  than  saunpling  it  in  x.  In  the  latter 
case,  we  have  seen  that  an  array  of  output  detectors  is 
required  which  records  the  system  line  spread  function  as 
the  input  6 (y-C)  explores  the  continuum  values  of  C-  When 
sampling  h(x;  C)  in  however,  it  is  necessary  to  record 
the  entire  output  field.  An  input  S (y-^^)  appears  at  the 
output  as  h(x  - which  is  a sample  line  spread 

function.  The  Dirac  delta  is  then  translated  by  a samp- 
ling interval  corresponding  to  the  next  value  of  n and 
the  process  is  repeated.  An  example  of  a line  spread 
function  saunpled  in  this  manner  is  in  Fig.  2.5.  It  is 
important  to  note  that  the  condition  that  h(x;  C)  is  band- 
limited  in  C does  not  necessarily  imply  that  h(x-C;C)  is 
bandlimited  in  This  point  will  be  discussed  shortly. 

Consider,  next,  the  implementation  of  Eq.  2-84.  Again, 
our  space-variant  system  is  represented  as  the  sum  of  a 
number  of  isoplanatic  systems.  The  n^^  component  isopla- 


natic  system  has  a line  spread  function  h(x  - and 

is  fed  a sample  value  of  the  input  of  the  form  f <S  (y-C^^)  . 
The  outputs  from  the  component  systems  are  added  and  then 
passed  through  a low  pass  filter  of  bandwidth  2W  which 

S 

supplies  the  required  convolution  sine  in  Eq.  2-84. 

2.5.3.  Discrete  (Matrix)  Representations 

We  have  discussed  two  space-variant  sampling  theorems; 
one  when  h(x;  C)  is  bandlimited  in  x and  the  other  when 
both  h(x;  C)  and  f(C)  are  bandlimited  in  C*  In  both  cases, 
the  corresponding  system  outputs  are  bandlimited  and  thus 
can  also  be  expressed  in  a sampling  theorem  expansion.  In 
this  section,  we  will  be  concerned  with  directly  determin- 
ing these  scunple  values  of  the  system  output  under  various 
conditions . 

Consider  first  the  case  where  the  input  is  bandlimited 

and  h(x;  C)  is  bandlimited  in  C-  From  the  rect  term  in 

Eq.  2-83,  we  conclude  the  system  output  is  bandlimited 

with  bandwidth  2W  . As  such,  we  can  express  the  output 

s 

in  the  following  seunpling  theorem  expansion: 

g(x)  = Jg(Xjj^)sinc  2Wg(x  - x^^^)  (2-85) 

m 

where  here 

X = 

*m  2W  ' 
s 

We  now  define  the  low-pass  line  spread  function 
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h(x  - h(x  - )*sinc  2W  x . 

n n s XI  ii  0 


(2-86) 


The  coefficient  2W  is  used  to  assure  that  both  h(*,*)  and 

s 

h(»,*)  have  identical  units.  Using  the  low-pass  line 
spread  function,  the  sampling  theorem  expression  in  Eq. 
2-84  can  be  written 


g(x)  = ^ I f(Cj,)h(x  - . (2-87) 


s n 

The  expression  for  our  output  sample  values  follows  immed- 
iately as 

1 


9(V  -2S-  I V «n>- 

s n 


(2-88) 


This  relationship  can  be  viewed  as  a discrete  version  of 
the  superposition  integral  in  Eq.  2-8.  The  computational 
form  is  that  of  an  infinite  matrix.  That  is,  an  infinite 
vector  composed  of  the  sample  input  values  multiplies  an 
infinite  matrix  of  the  sample  line  spread  function  values 
to  give  an  infinite  vector  composed  of  the  corresponding 
output  sample  values. 

Let's  now  investigate  the  corresponding  case  when 
the  system  line  spread  function  is  bandlimited  in  x.  From 
the  rect  term  in  Eq.  2-71,  the  system  output  is  bandlimited 
in  X with  bandwidth  2W^  and  can  thus  be  expressed  in  the 
sampling  theorem  expansion  in  Eq.  2-85.  But  it  is  neces- 
sary to  redefine  x^  in  this  relationship  as 
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X = -5!- 
2W 


(2-89) 


In  order  to  determine  the  output  sample  values  in  this 
case,  we  rewrite  Eq.  2-72  in  integral  form: 


g(x)  =*  I / f(?)h(x  ; ?)sinc  2W  (x  -C  - x )dC. 

li  A 14 

n — ® 


(2-90) 


It  follows  that 


g(Xm)  = I / f(?)h(x^;  Osinc  2W^(Xjjj  “ x^  - C)dC 
n -® 


= ^ (f(x)h(x  ; x)]*sinc  2W  x (2-91) 

n I X = X - X . 

m n 


As  can  be  seen,  the  discrete  treatment  for  this  case  does 
not  in  general  result  in  the  computational  attractive  ma- 
trix form  encountered  in  Eq.  2-88. 

Computational  simplifications  occur  in  the  discrete 
treatment  for  systems  in  which  1)  the  input  is  bandlim- 
ited  and  2)  the  line  spread  function  is  bandlimited  in 
both  X and  Consider  first  the  spectriim  of  the  convo- 
lution term  in  Eq.  2-91: 


Fjj(f  (x)h(x^;x)*sinc  2W^x] 


2^  F^[f  (x)h(x^;x)  ]rect  (^~) . 


Consider  then  the  case  where,  from  Eq.  2-79,  h(x;C)  and 
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f(5)  are  also  bandlimited  in  C; 

fx 

F^[f  (x)h(x^;x)  ] = F^[f  (x)h(x^;x)  Irect  (2~) . 

*(2-92) 

If  W < W , then  the  rect(f  /2W  ) term  in  Eq.  2-92  serves 

S X XX 

no  purpose.  That  is,  the  term  does  not  chop  off  higher 
frequencies  simply  because  there  are  no  frequencies  great- 
er than  W < W . For  this  condition,  we  conclude  that 
s X 

f (x)h(x^;x)  *sinc  2W^x  = 2^  f (x)h(Xj^;x)  ; • 

(2-93) 

Thus,  our  sampling  theorem  expression  in  Eq.  2-91  becomes 


9<v  = 5^  j ''=‘m  - % -='n' ’ "s  ' "x  • 

X n 

(2-94) 

Malting  the  index  change  ft  = m - n gives  the  final  desired 
relationship  for  this  case: 


9<V  * 2^  I ' "s  ^ ”x- 

* fi 

We  have  again  arrived  at  an  infinite  matrix  relationship. 
Mote  that  we  are  again  oversampling  but  in  a different 
manner.  Before  sampling  was  performed  at  a rate  of  2W 

5 

in  the  C direction.  In  Eq.  2-95,  sampling  is  performed 
at  a rate  W > W in  the  ? direction  so  that  sampling 

X s 

indexes  will  match. 

Let  us  now  apply  a similar  approach  to  the  expansion 
in  Eq.  2-88  in  which  the  input  is  bandlimited  and  h(x;C) 
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is  bandlimited  in  C-  Consider  the  low-pass  line  spread 
function  as  defined  in  Eq.  2-86  for  the  case  where  W <W„. 
Employing  arguments  identical  to  those  just  used,  we  con- 
clude that 

h(x  - X ; X ) = h(x  - x ; x ) ; W > W (2-96) 
n n n n s X 

in  which  case,  Eq.  2-88  becomes 

^ J "n*  ’ "s  ^ ”x- 

s n 

When  f(C)  and  h(x;C)  are  bandlimited  in  both  x and 
C,  either  the  expression  in  Eq.  2-95  or  in  Eq.  2-97  applies 
depending  on  which  bandwidth,  W or  W , is  larger.  We  can 

S X 

combine  these  two  relations  into  a single  expression  by 
writing 

9(x„)  = 2^  I 5„)  (2-98) 

n 

where 

^m  = ^m  = "‘/2W  (2-99) 

and 

W = max[W^,W^].  (2-100) 

X s 

In  all  cases  of  Eq.  2-98,  we  are  oversampling.  Sup- 
pose that  we  wished  to  Scunple  at  the  minimum  allowable 
rate.  That  is,  sample  f(?)  at  a rate  of  2Wj  and  sample 
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h(x;  C)  at  a rate  of  2W^  in  x and  2W^  in  C*  As  is  shown 
in  Appendix  B,  the  sample  value  of  the  output  is  given 
in  this  case  by 

“III  ^^2W”'2W“^^pqr  ^2W^  (2-101) 

rpg  r x v 

where 

Ipqj,(x)  = [sinc(2W^x  - p)  sinc(2W^x  - q)]*sinc  (2W^x  - r) . 

(2-102) 

Clearly,  the  discrete  superposition  in  Eq.  2-98  is  compu- 
tationally much  simpler  than  the  above  triple  sxim  of  con- 
tinuous convolutions.  The  price  we  pay  for  this  reduction 
in  complexity,  however,  is  the  increased  Scunpling  rate. 

2.5.4.  Alternate  Saunpling  Theorems 

We  have  developed  two  basic  sampling  theorems  appli- 
cable to  two  classes  of  linear  systems.  The  first  is 
applicable  to  systems  in  which  h(x;  C)  (and  thus  h(x-C;C)) 
is  bandlimited  in  x.  The  second  is  for  variation  limited 
systems  with  a bandlimited  input.  Clearly,  these  two 
linear  system  classes  form  only  two  small  subsets  of  the 
set  of  all  linear  systems. 

It  is  possible,  however,  to  widen  this  class  by 
utilization  of  various  forms  of  the  system  line  spread 
function  which  might  be  bandlimited  in  one  or  more  vari- 
ables. For  exeimple,  if  h(x  - 5;  K)  has  finite  support 


in  C for  all  x,  it  follows  that  its  Fourier  transform 
with  respect  to  C will  be  bandlimited  in  v for  all  x. 

The  Fourier  transform  of  the  line  spread  function,  as  is 
given  in  Eq.  2-22,  is  simply  the  system  frequency  re- 
sponse. Thus  we  can  utilize  the  frequency  response  char- 
acterization in  Eq.  2-20  and  develop  a corresponding 
sampling  theorem.  For  another  system,  h(x  - 5)  might 

be  bandlimited  in  C.  Note  that,  from  Fig.  2.1,  the  con- 
dition of  h(x;  C)  being  bandlimited  in  ^ does  not  neces- 
sarily assure  that  h(x  - C;  C)  is  bandlimited  in  At 
any  rate  this  is  a condition  not  yet  discussed  for  which 
a Scunpling  theorem  can  be  easily  derived. 

In  this  section,  we  will  present  a brief  rundown  of 

these  and  other  possible  input/system  conditions  where 

the  sampling  theorem  is  applicable.  The  list  is  in  no 

way  meant  to  be  exhaustive.  There  are  forms  of  the  line 

spread  function  that  we  here  have  not  even  considered 
25 

yet.  The  expressions  to  be  listed  are  rather  meant  to 
illustrate  the  relative  ease  in  application  of  sampling 
theorem  concepts  to  the  system  process.  For  the  purpose 
of  completeness,  the  two  S2unpling  theorems  previously 
discussed  are  included.  In  stating  these  seunpling  theorems, 
we  use  various  forms  of  the  system  line  spread  function 
such  as  the  frequency  response,  transfer  function,  varia- 
tion spectrum,  etc.  Throughout,  we  use  the  notation 


F(v)  = F^[f(C)].  All  of  these  sampling  theorems  are 
special  cases  of  the  superposition  integral: 


g(x)  = / f(C)h(x  - C;  5)dC.  (2-103) 

«ao 

1)  For  the  case  where  h(x;  C)  is  bandlimited  in  x,  we 
rewrite  from  Eq.  2-65: 

f 

Hx(fx;  o = H^(f^;  ?)rect(^).  (2-104) 

This  relation  is  a sufficient  condition  for  the  sampling 
theorem  expansion  in  Eq.  2-73  repeated  here: 

g(x)  = I [f(x)h(x  ; x)]*sinc  2W  (x  - x ) (2-105) 

^ II  X li 

n 

where  x^  = n/2W^. 

2)  The  Fourier  dual  of  the  above  occurs  when 

Kx(fx»*  v)AF^[k(x;  v)  ] 

f 

= K^(£x'  '>)rect(^).  (2-106) 

Note  that  this  condition  is  equivalent  to  that  in  Eq.  2-104. 
That  is,  if  h(x;  C)  has  bandwidth  2W^  in  x,  then  so  does 
k(x;  C) • Under  the  condition  of  Eq.  2-106,  our  sampling 
theorem  expansion  is  simply  the  Fourier  dual  of  Eq.  2-105: 

g(x)  = I ^ [F(|)k(x^;|)]*sinc  2W^(x  - x^^) 
n 

where,  again,  x^  = n/2W^. 


(2-107) 


3)  For  the  case  where  both  £(C)  and  h(x;  are  band- 
limited  in  our  conditions  are 

H^(x;  V)  *H^(x;  v)rect(-2^)  (2-108) 


and 


F(v) 


F (v) rect ( 


(2-109) 


From  Eq.  2-84,  our  corresponding  sampling  theorem  expres- 
sion is 


g(x)  = J f(C^)h(x  - 5^;  C„)*sinc  2WsX  (2-110) 

n 

where  Wg  » and  5^^  * n/2Wg. 

4)  The  conditions  for  the  Fourier  dual  of  the  above 
S€unpling  theorem  expansion  are 
00 

K^(x;  C)A  / k(x;  v) e'^^’^'^^dv 
*00 

= K^(x;  C)rect(^)  (2-111) 

and 

f(C)  = f(C)rect(^).  (2-112) 

It  follows  that  the  Fourier  dual  of  Eq.  2-110  is  then 

9(x)  “ § 1 F(v^)k(x  - cVj^;  v^)*sinc  (2TgX/c) 
n 

where  T„  = T,,  + T-  and  = n/2T„. 
s V t n s 


(2-113) 


5)  We  here  consider  the  case  where  h(x  - C;  C)  is  band- 
limited  in  5.  Since  the  frequency  response  is  the  Fourier 
transform  of  h(x  - C;  C) » it  follows  that  this  condition 
is  satisfied  when 


k(x  - cv;  v)  = k(x  - cv;  v)rect(=^).  (2-114) 

We  must  also  restrict  the  input  to  be  bandlimited: 


F(v)  = F (V)  rect  (v/2Wj^) . 


(2-115) 


It  then  follows  from  previous  arguments  that  the  product 
f(C)h(x  - C?  C)  has  bandwidth  2W^  = 2Wj^  + 2Wj.  Thus,  we 
can  write; 

f(C)h(x-C;C)  = I f (Cn)h(x-Cj^;Cjj)sinc  2W^(C-?j^)  (2-116) 

n 

where  = n/2W^.  Integrating  over  all  C yields 

00 

g(x)  = /_„f (C)h(x-C;C)dC 

00 

= I f(Cn)h(x-Cj^;Cn)/_„sinc  2W^(C-C^) 


(2-117) 


This  is  the  desired  sampling  theorem  expansion  for  the 
conditions  in  Eqs.  2-114  and  2-115. 

6)  The  conditions  for  the  Fourier  dual  of  the  above  saunp- 
ling  theorem  expansion  are 
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h(x-C;?)  = h(x-C;C) 

and 

f(C)  = f(C)rect(i).  (2-119) 

The  Fourier  dual  of  Eq.  2-117  follows  immediately  as 

g{x)  * I F(v^)k(x  - cv^;  v^)  (2-120) 

t n 

where  T.  = T.  + T,  and  » n/2T.  . 
t h r n t 

2.6.  Characterization  Comparisons 

In  this  chapter,  we  have  explored  various  methods  by 
which  a linear  system  can  be  characterized.  In  all  cases, 
this  characterization  is  performed  by  probing  the  system 
with  various  inputs  and  observing  the  corresponding  system 
response.  Under  certain  system  and/or  input  class  con- 
straints, the  set  of  sufficient  conditions  required  to 
completely  characterize  the  system  performance  varies. 

For  the  general  linear  system  with  no  restriction 
on  the  input  class,  cataloging  the  system  response  to  all 
elements  in  a continuum  orthonormal  basis  set  sufficies  for 
complete  system  characterization.  Special  cases  of  this 
scheme  include  use  of  the  Dirac  delta  orthonormal  basis 
set  in  which  case  the  system  is  characterized  by  the  super- 
position integral.  The  system  response  to  each  element 
in  the  Dirac  delta  class  is  termed  the  system  line  spread 


function.  Another  orthonormal  basis  set  consists  of  the 
complex  exponential  whose  input  into  the  system  results 
in  the  system  frequency  response.  By  appropriate  manipu- 
lations of  the  frequency  response  and  line  spread  function 
characterizations  with  Fourier  transform  operations , it 
was  shown  that  the  output  spectrum  in  each  case  could  be 
directly  expressed.  The  ease  in  performing  these  manipu- 
lations was  due  primarily  to  our  particular  choice  of 
line  spread  function  and  frequency  response  notation. 

Three  characterization  schemes  in  which  restrictions 
were  placed  on  the  system  and/or  input  were  also  discus- 
sed in  this  chapter.  These  consisted  of  the  discrete 
orthonormal  basis  set  response  (DBR)  characterization, 
the  piecewise  isoplanatic  approximation  (PIA) , and  the 
seunpling  theorem  characterizations.  In  each  case,  the 
continuum  number  of  system  responses  required  in  the  gen- 
eral case  was  reduced  to  a countable  number.  In  certain 
instances,  the  system  characterization  reduced  to  the 
superposition  of  a number  of  invariant  systems.  Each  of 
the  three  schemes  is  characterized  by  distinct  properties 
that  we  shall  now  review  and  contrast. 

The  DBR  characterization  scheme  places  a restriction 
on  the  class  of  allowable  system  inputs.  That  is,  system 
inputs  are  restricted  to  those  signals  which  can  be  ex- 
pressed in  an  orthonormal  series  expansion  for  a given 
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orthonormal  basis  set.  The  scheme  is  thus  only  input  de- 
pendent and  can  be  utilized  for  any  linear  system.  The 
system  is  completely  characterized  with  knowledge  of  the 
system  response  due  to  each  element  in  the  discrete  ortho- 
normal basis  set. 

The  PIA  is  a scheme  that  is  input  independent.  That 
is,  the  restrictive  assumption  concerns  only  the  system. 
The  system  line  spread  function  is  assumed  to  have  suffi- 
ciently negligible  variation  with  respect  to  its  second 
variable  over  an  input  patch.  T1  system  can  be  then 
approximated  as  piecewise  isoplan  ic.  In  order  to  com- 
pletely characterize  a piecewise  i planatic  system,  we 
need  knowledge  of  the  system  line  spread  functions  cor- 
responding to  an  input  Dirac  delta  placed  within  each  of 
the  input  isoplanatic  patches.  The  resulting  system  char- 
acterization consists  of  a countable  number  of  invariant 
systems  whose  outputs  are  summed  to  give  the  approximated 
output  of  the  parent  linear  variant  system. 

The  final  method  discussed  for  characterizing  linear 
variant  systems  utilized  application  of  the  Scunpling 
theorem.  One  seunpling  theorem  requires  a bandlimiting 
assumption  on  the  system  line  spread  function  and  is 
thus  input  independent.  The  second  sampling  theorem  re- 
quires that  the  system  is  variation  limited  and  that  the 
system  input  be  bandlimited.  This  latter  scheme  is  thus 


both  input  and  system  dependent. 

All  of  the  system  characterizations  discussed  will 
be  applied  in  one  degree  or  another,  to  the  design  of 
coherent  optical  processors  for  performing  general  space- 
variant  operations  although  they  can  be  applied  to  any 
linear  system.  The  continuum  orthonormal  basis  set  re- 
sponse characterization  will  be  utilized  in  the  next 
chapter  for  design  of  one-dimensional  space-variant  pro- 
cessors. The  PIA,  DBR,  and  Scimpling  theorem  character- 
ization are  used  in  Chapter  IV  to  extend  generalized 
space-variant  processing  to  two  dimensions. 

We  note  finally  that  the  derivations  of  the  system 

characterizations  here  were  formal.  That  is,  they  lack 

mathematical  rigor.  Such  treatment  was  used  for  reasons 

of  continuity  and  clarity  of  presentation.  The  manipu- 

2 23  24  28 

lations  used  however,  are  relatively  standard^'  ' ' ' 

40,64,66  shown  to  be  valid  under  certain 

physical  realizability  criteria  which,  as  witnessed  by 
the  experimental  results  in  Chapter  III,  are  character- 
istic of  coherent  optical  processors. 


CHAPTER  III 


3.  ONE-DIMENSIONAL  SPACE-VARIANT  PROCESSORS 

In  this  chapter  we  shall  utilize  the  continuum  ortho- 
normal basis  set  response  characterizations  of  linear 
systems  to  develop  coherent  optical  processors  capedsle  of 
performing  a wide  class  of  one-dimensional  (1-D)  space- 
variant  operations.  Recognition  of  the  capability  of 
coherent  processors  to  perform  general  linear  integral 

type  of  operations  can  probably  be  credited  to  Cutrona 
1 12 

et.al.  ' More  recent  schemes  in  one- dimensional  process- 

37—38 

ing  include  those  of  Rhodes  et.al.,  whose  work  pri- 

marily was  concerned  with  frequency-variant  processing. 

The  bulk  of  the  one-dimensional  processing  schemes 

14 

in  this  chapter  are  due  to  the  work  of  Marks  et.al.  and 

. 13 

also  follow  closely  the  independent  work  of  Goodman  et.al. 

We  begin  this  chapter  with  a presentation  of  two  basic 

coherent  processing  schemes  capable  of  performing  a wide 

class  of  1-D  space-variant  operations.  The  first,  called 

the  direct  output  display  method,  performs  a one>dimensional 

space-variant  operation  and  displays  the  system  output 

along  a line  in  the  output  plane  of  the  coherent  processor. 

The  second  scheme,  termed  the  output  spectrum  display 

method  does  essentially  the  seune  thing,  except  the  Fourier 

transform  (or  spectrum)  of  the  system  output  is  displayed. 
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Both  schemes  utilize  a line  spread  function  mask  on  which 
the  system  line  spread  function  is  recorded.  Thin  lenses 
are  utilized  to  perform  appropriate  Fourier  transforma- 
tion and  imaging  operations.  Numerous  examples  and  experi- 
mental results  are  presented  for  various  forms  of  these 
processors. 

Due  to  the  limited  dynamic  range  of  photographic 
transparencies/  the  line  spread  function  masks  correspond- 
ing to  certain  linear  operations  might  possibly  not  be 
implementable.  It  might  be  possible,  however,  to  record 
the  corresponding  frequency  response  of  the  operation 
may  be  2d3le  to  be  recorded  on  a mask.  Thus,  it  is  in- 
structive to  consider  the  Fourier  duals  of  the  previous 
two  processors  thereby  increasing  the  class  of  space-var- 
iant systems  to  which  the  processor  techniques  are  appli- 
cable. 


3.1.  Direct  Output  Display  (POD)  Method 

The  direct  output  display  (DOD)  method  for  space- 
variant  processing  directly  evaluates  the  superposition 
integral  which  we  repeat  here  from  Eq.  2-8: 

OD 

* L„f<C)h(x  - C;  C)dC.  (3-1) 

A coherent  processor  capable  of  evaluating  the  superposi- 
tion integral  is  shown  in  Figure  3-1.  The  one- dimensional 


P2 


Figure  3.1  A coherent  optical  processor  for  perforaing 
1-D  space-variant  operations.  The  desired 
processor  output  lies  along  the  x axis  on 
plane  P2* 


r 
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input  f(C)  is  placed  in  plane  directly  adjacent  to  a 
mask  on  which  the  system  line  spread  function  h(x  - C) 
is  recorded.  Note  that  f(5)  covers  the  entire  (x,  C) 
plane  with  no  variation  in  the  x direction.  Cylindrical 
lenses  L2  and  have  respective  focal  lengths  of 

Thus,  when  plane  is  illuminated  from  the  left  with  a 
coherent  plane  wave,  Fourier  transformation  is  performed 
in  the  ^ direction,  and  imaging  is  performed  in  the  x 

direction. The  field  amplitude  ^^(x,  v)  on  the  output 

* 

plane  P2  is  then  given  by 

00 

gQ(x;  V)  = /_^  f(C)h(x  - C;  C)  exp (-j2TTV?) d? 

(3-3) 

where  the  spatial  frequency  v is  related  to  the  actual 
horizontal  distance  X2  on  plane  P2  by 

V = (3-4) 

Here,  X is  the  wavelength  of  the  spatially  coherent  illum- 
ination . 

*Here  and  henceforth,  constant  multiplicative  factors 
associated  with  field  cunplitude  relations  will  be  omitted 
since  it  is  the  functional  interelations  we  are  actually 
interested  in. 
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Comparing  the  superposition  integral  in  Eq.  3-1  with 
the  processor  output  [Eq.  3-3],  we  find  that 

g(x)  * gQ(x,0) . (3-5) 

That  is,  the  1-D  output,  corresponding  to  the  input  and 
line  spread  function  mask,  appears  in  the  processor  out- 
put plane  along  the  x axis.  The  desired  space-variant 
operation  has  thus  been  performed.  For  later  reference, 
this  processing  scheme,  or  modifications  thereof,  will 
be  referred  to  as  the  direct  output  display  or  DOD  method. 

An  alternate  and  somewhat  simpler  scheme  for  perform- 
ing 1-D  Fourier  transformation  is  pictured  in  Fig.  3-2. 

As  shown,  spherical  and  cylindrical  lenses,  placed  back 
to  back,  replace  the  three  cylindrical  lenses  in  Fig.  3-1. 
If  we  assign  a focal  length  of  f to  both  the  spherical 
and  cylindrical  lens,  this  processor's  output  intensity 
distribution  is  the  scune  as  in  Fig.  3-1  with  v=X2/Xf.^^ 

The  primary  limitation  of  this  processor  is  the 
allowable  dynamic  range  of  our  line  spread  function  mask. 
If  a holographic  mask  is  used  then  we  must  require  that 

the  values  of  h(x  - C;?)  lie  within  the  unit  circle  on 

23 

the  complex  plane.  As  we  shall  see  in  our  first  example, 
however,  this  limitation  can  sometimes  be  overcome  by 
clever  design. 


Figure  3.2  Another  coherent  processor  for  performing 

1-D  space-variant  operations.  The  intensity 
distribution  on  P-  is  identical  to  that  in 
Figiire  1 . ^ 
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3.1.1.  One- Dimensional  Magnification 

As  our  first  exeunple  for  specific  application  of 
the  DOD  processor,  we  consider  the  ideal  magnifier  which 
is  characterized  by  the  input-output  relationship 

g(x)  » ^f(x/M)  (3-6) 

where  M is  the  magnification.  From  the  sifting  property 
of  the  Dirac  delta  in  Eq.  2-5,  we  write 

00 

00 

= - MC)d5.  (3-7) 

The  corresponding  line  spread  function  is  thus 


J 


1 
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h(x  - C;  ?)  = 6(x  - M5).  (3-8) 

On  the  (x,  5)  plane,  we  interpret  this  relation  as  a Dirac 
delta  sheet  along  the  line  x = M?.  The  magnification  is 
simply  the  line's  slope. 

We  may  implement  the  DOD  magnifier  as  shown  in  Fig. 
3-3.  The  Dirac  delta  sheet  is  formed  by  focusing  a plane 
wave  with  a cylindrical  lens  mounted  on  a rotatable 
assembly.  We  can  change  the  slope  of  this  sheet,  and 
thus  the  resulting  magnification,  by  simply  rotating  the 
lens.  Note  here,  that  by  appropriate  design,  we  have  gen- 
erated our  line  spread  function  without  the  use  of  a mask 
and  thus  have  avoided  the  fundounental  limitation  of  a 


mask's  dynamic  range  in  which  a Dirac  delta  obviously  does 
not  fall. 

In  Fig.  3-3,  at  the  lens'  back  focal  plane,  we  place 
the  input  f(5).  As  before,  Fourier  transformation  is  per- 
formed in  the  C direction  by  cylindrical  lenses  L2» 
'and  Lj.  The  desired  output,  (1/M)  f (x/M) , appears  on 
plane  P2  along  the  x axis.  For  the  case  of  a double 
pulse  input,  experimental  results  using  this  setup  are 

shown  in  Fig.  3-4  for  various  M.  Note  that,  although 

23 

2-D  magnification  is  optically  trivial  no  conventional 
scheme  can  control  the  magnification  parameter  from  a 
single  plane  as  it  is  one  here. 

3.1.2.  Coordinate  Distortion 

The  magnifier  scheme  may  be  generalized  to  linear 
systems  with  input-output  relationships  of  the  form 

g(x)  = f [D(x)  ] (3-9) 

where  we  shall  refer  to  D(x)  as  the  distortion  function. 
From  the  sifting  property  of  the  Dirac  delta,  we  can 
write 

00 

g(x)  = /_^  f(C)6[5  - D(x)]dC.  (3-10) 


The  line  spread  function  associated  with  Eq.  3-9  is  thus 
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h(x  - 5;  5)  = 6[C  - D(x)l.  (3-11) 

Our  Dirac  delta  sheet  is  now  bent  along  the  locus  C = D(x) 
on  the  (x,  C)  plane.  Such  a bent  Dirac  delta  may  be  crudely 
generated  by  an  appropriately  bent  glass  rod  which  acts  as 
a curved  cylindrical  lens.  Experimental  results,  using  a 
glass  rod  with  a single  bend,  are  shown  in  Fig.  3-5.  The 
resulting  lens  slopes,  as  shown  in  Fig.  3-6  (a),  constitute 
different  magnifications  for  C > 0 and  ? < 0.  The  DOD 
processor  output  for  the  input  pulse  configuration  in  Fig. 
3-6 (b)  is  shown  in  Fig.  3-5. 

37 

An  alternate  form  of  distortion,  used  by  Rhodes, 
employs  the  line  spread  function 

h(x  - C;  5)  = exp[-j2TrCD(x)  ].  (3-12) 

Here,  the  system  output  is  a distorted  version  of  the  in- 
put's spectrum 

00 

“ /_«  f(C)exp[-j27rCD(x)]dC 

= F[D(x)]  (3-13) 

where  F(v)  is  the  Fourier  transform  of  the  input 

00 

“ /_«of  (C)exP(“j2TrvC)dC  . 


(3-14) 
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Figure  3.5  Output  of  a piecewise  magnifier  distortion 
processor  corresponding  to  the  inputs 
pictured  in  Fig.  3.6.  The  larger  pulse  is 
roughly  three  times  the  length  of  the 
smaller. 


s 
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Rhodes uses  such  a system  for  generating  a log-frequency 
display  of  the  input  by  utilizing  a distortion  function 
proportional  to  In (x) . 

3.1.3.  Convolution  and  Correlation 

The  convolution  operation  is  a special  case  of  the 
superposition  integral  that  arises  when  the  linear  system 
is  space- invariant  (or  isoplanatlc) . From  Eq.  2-10,  the 
one-dimensional  convolution  integral  is 

9(x)  - f(5)h(x  - C)d5 

- f(x)*h(x).  (3-15) 

The  capability  of  conventional  coherent  processors  to  per- 

23 

form  t%io» dimensional  convolution  is  well  kno«m.  All  of 
these  schemes,  however,  either  necessitate  that  the 
Fourier  transform  of  the  line  spread  function  (system 
transfer  function)  be  recorded  on  an  amplitude  transmit- 
tance or  require  motion.  These  requirements  inhibit  real 
time  application.  As  we  will  show,  1-D  convolution  can 
be  performed  with  a variation  of  the  DOD  processor  with 
no  requirement  of  motion  or  Fourier  encoding. 

Before  presenting  the  DOD  processors  capable  of  per- 
forming convolution  and  correlation,  we  digress  briefly 
for  a discussion  of  the  geometrical  interpretations  of 
the  spatial  shift  encountered  in  Eq.  3-15.  Consider  the 
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1-D  function  r(C)  in  the  (x,  C)  plane  as  pictured  in  Fig. 

3-7 (a) . If  we  rotate  this  transparency  about  the  origin 
through  an  angle  of  6 in  the  clockwise  sense  [Fig.  3-7 (b) ] , 
the  resulting  2-D  function  is  described  by 

r[Ccos6  - X sinei.  (3-16) 

For  the  case  of  a 45®  rotation,  we  obtain  r[(C  - x)//T~  ]* 

To  eliminate  the  factor,  scaling  can  be  trivially  accom- 
plished by  a conventional  (2-D)  imaging  system  with  mag- 
nification M = l//3~.  Consider  next  physically  rotating 
the  transparency  in  Fig.  3-7 (b)  180®  about  both  ibs  x and 
C axes.  This  constitutes  coordinate  reversal.  For  the 
case  of  6 > 45®,  the  result  would  be  the  transmittance 
r[(x  - which  is  the  scaled  shift  required  by 

the  convolution  integral.  Various  other  rotations  could 
of  course  be  employed  to  obtain  a number  of  such  shifts 
and  scalings. 

Consider,  then,  the  DOD  convolution  processor  as  pic- 
tured in  Fig.  3-8.  In  plane  we  place  the  transparency 
representing  h(^)  with  the  45®  rotation  shown.  This  func- 
tion, hI(C  - x)//5“]  is  scaled  and  inverted  by  lenses  L 
and  which  have  respective  focal  lengths  of 


(3-17) 


Figure  3.7 


(a)  A 1-D  function  r(C)  on  the  (x,C) 

plane.  (b)  The  transmittance  in  (a) 
rotated  clockwise  about  the  origin 
and  auigle  of  9.  (c)  The  transmittance 

in  (b)  for  6 >■  45*  with  a coordinate 
reversal  formed  by  physically  rotating 

(b)  about  its  C and  x 2oces. 


scaling  on  the  rotated  transmi 


Incident  on  plane  P2  is  the  desired  h(x  - C)  which  multi- 
plies the  input  transmittance  £(C)>  The  product  is  then 
processed  as  before  by  the  three  cylindrical  lenses.  The 
output  (£  * h)  appears  c\long  the  x axis  in  plane  P^. 

A primary  drawback  to  this  processing  scheme  is  find 
ing  two  lenses  which  have  a focal  length  ratio  of  /?.  An 
alternate  DOD  processor  exists,  however,  that  does  not 
require  the  use  of  these  scaling  lenses.  We  will  illus- 
trate such  a processor  that  performs  the  correlation 
operation* : 

g(x)  = /“„f(5)h(C  - x)dC 

» f(x)*h(x)  (3-18) 

where  denotes  correlation.  Comparison  with  Eq.  3-15 
reveals  similar  computational  forms  for  convolution  and 
correlation. 

Consider,  then,  placing  the  1-D  transmittances , f(C) 
and  h(C),  in  plane  Pj^  of  the  DOD  processor  in  Fig.  3-1 
each  rotated  45®  in  such  a manner  as  to  form  the  product 

f [ iii  j h [ ^ ] . (3-19) 

/?  /7 

The  field  amplitude  in  plane  P2»  after  the  1-D  Fourier 
transformation  [Eq.  3-3] , is 

*For  complex  signals,  h(*)  is  complex  conjugated  in  this 
integral  relation. 
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g_(x;  V)  = /”  f (■^^)h(-^^)exp(-j2TrvC) 

° /I  /I 

» /?  exp(j27rvx)  /"^  f(C') 

h[C'  - /Ixlexp[-j2Tr(/5v)C']<i?'  (3-20) 

where  we  have  made  the  variable  substitution  5'  = C^+x)//?. 
Along  the  x axis,  we  obtain  the  output 

g(x)  * 0) 

= /5  /“„  f(C')h(C  - /Jx)dC'.  (3-21) 

This  relationship  is  recognized  as  a scaled  version  of 
the  correlation  integral  and  is  thus  the  desired  result. 
Convolution  can  of  course  be  performed  in  a similar  man- 
ner by  choosing  appropriate  orientations  of  the  1-D  input 
transparencies . 

3.2.  Entire  Output  Plane  Utilization 

The  purpose  of  the  DOD  processor  is  to  compute  one- 
dimensional linear  integral  operations.  The  solution 
appears  along  a line  in  the  processor  output  plane.  There 
do  exist  some  computational  manipulations,  however,  in  which 
the  entire  output  plane  can  be  utilized.  These  cases  are 
serendipities  of  the  DOD  processor  scheme.  We  briefly 
consider  here  two  cases;  Laplace  transform  and  ambiguity 
function  display. 
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3.2.1.  Laplace  Transform  Display  ^ 

i' 

For  the  Laplace  transform  scheme,  our  line  spread  | 

function  mask  is 

h(a  - 5)  = exp(-aC)y(C)  (3-22) 

where  u(C)  is  the  unit  step  function:  j 


Equation  3-23  is  immediately  recognized  as  the  one-sided 
Laplace  transform  of  f(C)  with  frequency  variable 

s = a'+  jB.  (3-26) 

For  this  case,  the  (a,  6)  output  plane  of  the  DOD  processor 
may  be  interpreted  as  the  complex  s plane  encountered  in 
Laplace  transform  theory.  A discussion  of  this  specific 

*5  0 

DOD  processor  is  given  by  Mueller  and  Carlson.*’^ 

3.2.2.  Ambiguity  Function  Display 

40 

The  aunbiguity  function,  first  introduced  by  Woodward, 
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has  been  applied  in  radar  in  predicting  the  capability  of 
a given  signal  to  determine  simultaneously  the  range  and 
velocity  of  a target.  The  range  is  determined  by  the  time 
delay  x and  the  velocity  by  the  Doppler  shift  v.  The  eunbi- 
guity  fvinction  for  a given  real-valued  signal  f (x)  is 

X(v;  X)  - /•„  f(5)f(5  - x)exp(-j2iTv5)dC.(3-27) 

In  optics,  Papoulis  has  employed  the  eunbiguity  function 

in  analyzing  diffraction  phenomena. Cutrona  et.al.^'^^ 

2 

and  Preston  have  proposed  a coherent  ambiguity  function 

processor  that  utilizes  multiple  channels  to  display  the 

cunbiguity  function  for  discrete  values  of  x.  The  scheme 

42 

of  Casasent  et.  al.  generates  1-D  slices  of  the  eunbiguity 

function  in  the  (v,x)  plane.  Similar  1-D  displays  have 

43 

also  been  electronically  produced.  Use  of  the  DOD  pro- 
cessing scheme,  initially  reported  by  Said  and  Cooper^^ 
and  also  independently  discovered  by  Marks,  Walkup,  and 
Krile,  1)  displays  lx(v,x)l  in  a continuous  (rather 
than  quantized)  form  over  the  entire  (v,x)  plane,  2)  has 
the  capacity  for  extension  to  real  time  processing,  and 
3)  is  easily  implemented. 

Inspection  of  Eq.  3-27  reveals  that  it  is  the  compu- 
tational. form  of  the  output  of  the  DOD  processor  in  Eq. 

3-3  for  an  input  of  f(5)f(5  - x) . Furthermore,  the  input 
is  recognized  as  the  integrand  of  the  correlation  integral 


1. 


in  Eq.  3-18  for  the  case  of  autocorrelation . Thus,  when 
we  compute  an  autocorrelation  of  a signal  f(C)  as  pre- 
viously described,  we  are  also  computing  the  signal's 
eunbiguity  function.  Recall,  however,  when  the  two  trans- 
mit tances  were  both  appropriately  rotated  45®  in  the  input 
plane,  the  resulting  processor  output  [Eq.  3-20]  has  a 
spatially  varying  phase. For  identical  inputs  the  out- 
put field  amplitude  is  then  given  by 

go(x;  v)  = /7exp( j2iTvx) 

00 

/-oo  - .^x)exp[-j2Tr(/Iv)C]dC. 

(3-27) 

In  dealing  with  the  ambiguity  function,  however,  one  is 
primarily  concerned  with  the  modulus  squared  of  x ( * » * ) • 

The  spatially  varying  phase  term  in  Eq.  3-27  is  thus  of  no 
concern  since 

lgQ(x;  v)|^  = 2l/“^  f(C)f(C  - /Ix)exp[-j2Tr(/Iv)C]d5l2 
= 2|x(/2’v;  /Ix)  1^,  (3-28) 

The  intensity  distribution  on  the  output  plane  is  thus  a 
scaled  version  of  the  desired  square  modulus  of  the  2unbi- 
guity  function. 

We  will  now  evaluate  the  performance  of  the  ambiguity 
fvinction  processor  by  comparing  the  processor  output  with 
the  theoretical  results  for  both  single  and  double  pulse 
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inputs.  For  a single  pulse  [Fig.  3-9(a)],  we  may  write 


f(C)  » rect(:*~) 


(3-29) 


where  2T  is  the  pulse  duration.  The  geometric  interpre- 
tation of  f(5)»  f((C  + x)//7j,  and  f ( (C+x)//J]f  [ (5-x)//J] 
are  shown  in  Figs.  3-9 (b),  9(c),  and  9(d)  respectively. 

Substituting  Eq.  3-29  into  Eq.  3-27  followed  by  eval- 
uation yields  the  ambiguity  function 
► 

X (v;x)  » •(  (2T-|xl ) sincv(2T-|xl  )exp(-jirvx) ; |xl<2T 
0 ; |x|>2T. 

» 

(3-30) 

The  corresponding  output  intensity  is 

lx(v;x)l^  * i (2T-lx| ) ^sinc^v (2T-lxl ) ; lxl<2T 

“ (3-31) 

0 ; lxl<2T. 

For  purposes  of  identification,  it  is  instructive 
to  examine  the  locus  of  points  where  the  ambiguity  func- 
tion is  identically  zero.  From  Eq.  3-31,  this  zero  locus 
may  easily  be  shown  to  be 

V = n/(2T  - |x|)  ; jxj  < 2T  (3-32) 


where  n is  any  nonzero  integer.  The  piecewise  hyperbolic 
nature  of  these  curves  is  shown  in  Fig.  3-10. 


K 2T-H 


Figure  3.9  A single  pulse  (a)  in  time,  (b)  in  the  (x,C) 
plane,  (c)  rotated  45*  on  the  (x,C)  plane  to 
form  £C(x  + C)//2I,  (d)  the  product  of  two 
pulses  rotated  45*  and  -45*  on  the  (x;C)  plane 
to  form  f t(5+x)//2]f  [(5-x)//2] . 
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The  ambiguity  fund  ion  for  a single  pulse  is  gener- 
% 

ated  by  appropriately  rotating  two  identical  thin  slits 

in  plane  of  the  coherent  optical  processor  of  Fig.  3-1. 

The  resulting  output  is  shown  in  Fig.  3-11.  As  can  be 

seen,  the  coherent  processor  output  compares  quite  nicely 

with  the  theoretical  result  in  Fig.  3-10.  A 3-D  computer 

graph  of  the  corresponding  eunbiguity  function  modulus  can 

44 

be  found  in  Fig.  6-6  of  Rihaczek. 

For  a double  pulse  input  [Fig.  3-12 (a)],  we  write 

f(5)  = rect[^^]  + rect[^^]  (3-33) 


where,  for  convenience,  the  pulse  separation  2T  has  been 
chosen  to  be  equal  to  each  pulse  width.  The  geometrical 
interpretation  of  f ( (C+x)/y^] f [ (5“X)//5]  is  shown  in  Fig. 
3-12 (b).  The  aunbiguity  function  associated  with  the  double 
pulse  is 


X (v,x) 


, 2 (2T- 1 xj ) sincv  (2T- 1 x| ) cos  (4i7Tv)  1 x ] <2T 

- (2T- 1 x| ) sincv  (2T- 1 x]  )e'’^^^'’;  2T£|x  |£4T 
(6T-|x| ) sincv  (6T- 1x1  4T<  lxl£6T 


The  corresponding  output  intensity  is 


; 1x1 >6T. 

(3-34) 


X (v,x) 


2 


4 (2T-lxl ) ^sinc^v  (2T- 1 X 1 )cos^  (4itTv)  ; lxl£2T 
(2T-|x| )^sinc^v(2T-|x| ) ;2T£|x|£4T 
(6T-lxl ) ^sinc^v(6T-lxl ) ; 4T£lxl<^6T 
0 ; lxl>6T.  (3-35) 


Figure  3.11  The  cunbiguity  function  (modulus  squared) 
display  for  a single  pulse. 


The  equations  describing  the  zero-value  loci  are  easily 
shown  to  be 
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(2m  + 1)/8T 
n/(|xl  - 2T) 
n/(6T  - 1x1) 


1 X 1 <2T 
I X I £4T 
4T£lxl£6T. 


(3-36) 


where  m is  any  integer,  and  as  before,  n -is  any  nonzero 
integer.  An  illustration  of  these  zero-value  loci  is 
offered  in  Fig.  3-13. 

By  appropriately  placing  two  identical  double  slits 
in  plane  of  the  coherent  processor,  the  ambiguity  func- 
tion for  the  double  pulse  is  generated.  The  result,  shown 
in  Fig.  3-14  again  compares  quite  favorably  with  the  theory. 

One  final  note  concerning  the  ambiguity  function. 

If  it  is  required  to  compute  it  without  the  spatially  vary- 
ing phase,  we  can  straightforwardly  do  so  by  using  the 
spherical  scaling  lenses  with  focal  length  ratio  of  /2 
and  three  cylindrical  lenses  as  discussed  previously. 


3.3.  Output  Spectrum  Display 

We  now  return  to  the  topic  of  generalized  one* dimen- 
sional space-variant  processors.  The  DOD  processor  has 
been  shown  to  compute  directly  the  output  of  a 1-D  linear 
space-variant  system.  An  alternate  scheme,  which  we  shall 
call  the  output  spectrum  display  or  DOD  method,  directly 
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We  have  shovm  that  this  relation  can  be  expressed  as  [from 
Eq.  2-25]  : 

G(f  ) = F-F^[f (C)h(x;  C)]  (3-38) 

^ V = fx 

where,  for  a function  p(x;  ?) , we  have  defined  the  Fourier 
operation: 

F^(p(x;  Ol  = p(x;C)exp(-j2TrvC)d5.  (3-39) 


Equation  3-38  states  that  we  may  generate  the  output 
spectrum  of  a space-variant  system  by  successive  Fourier 
transformation  of  the  product  f(C)h(x;C)  with  respect  to 
C and  X followed  by  evaluation  along  the  line 

The  coherent  processor  capable  of  generating  the  out- 
put spectrum  in  Eq.  3-38  is  pictured  in  Fig.  3-15.  The 
input  f(C)  and  the  line  spread  function  mask  h(x;  C)  are 
both  placed  in  the  front  focal  plane  of  spherical  lens 

In  the  back  focal  plane  we  obtain  the  familiar  Four- 

23 

ler  transform  relationship 


101 


exp[-j2Tr(fj^x  + vOJdxdC  (3-40) 

where  the  spatial  frequencies  are  related  to  the  horizontal 
(x^)  and  vertical  (y^)  distances  in  the  output  plane  by 

^x  “ v = xyAf.  (3-41) 

Note  that  we  may  write  Eq.  3-40  as 

Go(fx;  ^^)  = F^F^tf  (C)h(x;  C)].  (3-42) 

Comparing  with  Eq.  3-38,  we  conclude  that 

That  is,  the  desired  1-D  output  of  the  OSD  processor  in 
Fig.  3-15  lies  on  the  45*  line  v * f^  in  plane  P2.  Re- 
markably, this  familiar  Fourier  transform  configuration 
thus  has  the  capability  of  potentially  performing  a wide 
number  of  space-variant  operations  of  the  form  of  Eq.  3-38. 

When  one  is  interested  only  in  the  magnitude  of  the 
output  of  the  OSD  processor,  the  Fourier  transforming  lens 

may  be  placed  in  the  S2une  plane  as  the  input  and  line 

23 

spread  function  mask.  One  may  furthermore  completely 
discard  the  lens  by  simply  including  the  lens'  phase  trans- 
mittance in  the  line  spread  function  mask.  In  this  case. 


we  may  perform  1-D  space-variant  processing  with  an  input, 
a mask,  and  a few  centimeters  of  free  space.  Also,  note 
that  vignetting  would  be  eliminated  by  such  a scheme. 

We  now  present  some  specific  applications  of  the  OSD 
processor. 

3.3.1.  Magnifier  Spectrum  Display 

A straightforward  application  of  the  OSD  processor  is 
in  displaying  the  spectrum  of  an  ideal  magnifier.  We  be- 
gin by  rewriting  the  magnifier's  line  spread  function 
from  Eq.  3-8  as 

h(x;  5)  * 6tx  - (M-l)Cl  . (3-44) 

As  before,  we  may  generate  the  required  Dirac  delta  sheet 

with  a cylindrical  lens.  The  slope,  however,  is  now  equal 

to  M - 1.  The  OSD  processor  for  spectr\im  magnification 

is  pictured  in  Fig.  3-16. 

Experimental  results  displaying  the  spectrum  of  a 

magnified  rectangular  pulse  are  in  Fig.  3-17.  As  shown, 

the  resulting  sine  functions  are  inversely  magnified  due 

28 

to  the  scaling  theorem  of  Fourier  transform  theory. 

Note  that,  for  M > 0 (corresponding  to  a -45°  delta  sheet 
input) , the  magnified  pulse  is  a Dirac  delta,  the  spectrum 
of  which  is  uniform. 
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3.3.2.  Cross  Power  Spectral  Density  Display 

The  cross  power  spectral  density  function  of  two 

signals,  f (x)  and  h(x),  is  defined  as  the  Fourier  trans- 

29 

form  of  their  cross-correlation: 

* /-CO  If (x)*h(x)lexp(-j2irf^x)dx.  (3-45) 

Consider,  then,  placing  two  transparencies  of  f(x)  and 
s(x)  in  plane  of  Fig.  3-15  in  such  a manner  as  to  form 
the  product  f(x)h(C).  From  Eq.  3-40,  the  processor  out- 
put is 

®o^^x'  “ /-CO  + vOldxdC. 

(3-46) 

Along  the  45®  line  v * f^,  this  becomes 

°o^^x'^x^  “ /-oo  /ro»^I*)^I^^®*Pl"32iT(x+C)fjj]dxd5 

= J“„  /“^f  (x)h(C-x)exp(-j2TrC'f^)d5'dx 

(3-47) 

where  we  have  made  the  variable  substitution  C - ^ + x. 
We  may  rewrite  Eq.  2-47  as 

Go(fxrfx)  ff(C)*h(C)l 

- S(f^)  (3-48) 

which  is  the  desired  result. 


J 
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3.4.  Fourier  Duals  of  the  POD  and  OSD  Processors 

The  two  processors  thus  far  presented  utilize  line 
spread  function  masks.  The  DOD  processor  uses  a mask  of 
h(x  - C;  K)  to  compute  g(x)  while  the  OSD  processor -uses 
an  h(x;  C)  mask  to  compute  G(f^).  It  is  also  possible  to 
use  frequency  response  masks  in  one-dimensional  processor 
design.  The  processors  utilizing  such  masks  can  be  con- 
sidered Fourier  duals  of  their  spatial  domain  counterparts. 

We  investigate  use  of  the  frequency  response  mask  for 
two  basic  reasons.  First,  it  widens  the  ways  in  which  a 
given  space-variant  operation  can  be  performed  and  thus 
allows  us  to  choose  from  a number  of  possible  schemes  that 
one  which  is  most  easily  implemented.  Secondly,  due  to 
the  fundamental  dynamic  range  limitation  of  photosensitive 
media,  it  might  be  possible  that  a given  space-variant 
system  line  spread  function  mask  might  not  be  implementable 
whereas  the  corresponding  frequency  response  mask  might. 

Thus,  the  number  of  possible  space-variant  operations  which 
can  be  performed  by  the  one-dimensional  processors  is  in- 
creased. 

3.4.1.  The  FDOD  Processor 

The  Fourier  dual  of  the  direct  output  display  (FDOD) 
processor  is  pictured  in  Fig.  3-18.  This  processor  com- 
putes the  frequency  response  characterization  repeated 

< 

. ^ .j 
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here  from  Eq.  2-20; 

g(x)  = F(v)k(x  - cv;  v)dv,  (3-49) 

On  plane  of  the  processor , we  place  our  one-dimensional 
input,  f(C}>  Cylindrical  lens  L performs  a one-dimens- 

C 

* 

ional  Fourier  transform  so  that  we  have  F(v)  incident 

on  plane  P2.  Also  placed  in  this  plane  is  the  frequency 

response  mask  so  that  the  field  amplitude  immediately  to  > 

the  right  of  P2  is  I 

I: 

F(v)k(x  - cv;  v)  , (3-50)  | 

As  in  the  DOD  processor,  cylindrical  lenses  L2  and 

perform  a one-dimensional  Fourier  transform.  The  field 
amplitude  on  plane  P^  is  thus 

g_(x;  C)  * /*oo  F(v)k(x  - cv;  v)  exp  (-j2iTv5)  dv . 

(3-51) 

Comparison  with  Eq.  3-49  reveals  that 

g(x)  = go(x;  0).  (3-52) 

That  is,  our  desired  output  lies  along  the  x axis  in  the 
processor  output  plane. 

3.4.2.  The  FOSD  Processor 

The  Fourier  dual  to  the  output  spectrxim  display  (FOSD) 

*A  single  cylindrical  lens  can^be  used  due  to  the  fact 
that  f(C)  is  one -dimens ional.'^'*’ 

J 


processor  is  shown  in  Fig.  3-19.  This  processor  directly 
displays  the  Fourier  transform  of  the  space-variant  system 
output  according  to  Eq.  2-28  which  we  repeat  here; 
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G(f^)  =F"^F^[F(v)k(x;  v)  ] 


-cf 


= /!«  /!!coF(v)k(x;  V) 


exp  [-j2iTf^  (x  + cv)]dxdv. 


(3-53) 


On  plane  P^,  we  place  our  input  f(C)  which  is  again  Four- 
ier transformed  by  the  single  cylindrical  lens  L^.  The 
input  spectrviro,  F(v),  multiplies  the  frequency  response 
mask  placed  in  plane  P2  so  that  the  field  amplitude  immed- 
iately to  the  right  of  P2  is 


F(v)k(x;  v). 


(3-54) 


Spherical  lens  then  performs  a two-dimensional  Fourier 
transform  on  this  product.  The  field  amplitude  incident 
on  plane  P^  is  thus 

expt-j2Tr(f^x  + vOldxdv.  (3-55) 

Comparison  with  Eq.  3-53  reveals  that 
G(f^)  = G^(f^;  cf^). 


(3-56) 


Ill 


Thus,  our  desired  output  appears  in  plane  along  the  line 

In  summary,  we  see  that  the  Fourier  dual  FOOD  and 
FOSD  processors  differ  from  the  corresponding  Fourier 
dual  DOD  and  OSD  processors  in  only  two  respects.  First, 
a frequency  response  mask  is  used  instead  of  a line  spread 
function  mask.  Secondly,  the  system  input  is  Fourier 
transformed  before  being  processed.  In  all  other  respects, 
the  processors  are  implemented  identically. 

3.5.  Alternate  Interpretations  of  the  1-D  Processors 

In  presenting  the  OSD  and  FOSD  processors,  it  was 

stated  that  they  directly  computed  the  Fourier  transform 

of  the  system  output.  We  can  however,  ascribe  a quite 

45 

different  interpretation  to  their  operation.  Consider 
first  Eq.  3-38  which  we  rewrite  here  in  integral  form: 

G(fx)  = /!«  5)exp{-j2TrCfj^)dC.  (3-57) 

The  system  transfer  function,  H^(*,»),  was  defined  in  Eq. 
2-26  as  the  Fourier  transform  of  the  line  spread  function 
in  X.  Let  us  now  define 

- bC;  C)AH^(f^;  C)  exp  (-j2iTCf  (3-58) 

where  the  constant  b retains  dimensional  consistancy 
between  f^  and  C-  Substituting  into  Eq.  3-57  gives  a 


linear  integral  relation  akin  to  the  superposition  inte- 
gral: 
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= /-CO  - bC;  C)dC.  (3-59) 

Thus,  we  can  interpret  the  OSD  processor  operation  as  one 
that  1)  directly  computes  the  output  spectrum  of  the  lin- 
ear operation  corresponding  to  h(x  - C;  C)  or  2)  computes 
the  output  of  the  linear  operation  corresponding  to 
h^°Mfx  ” bC;  C)  . In  the  latter  case,  we  can  compute  the 
required  mask  transmittance  from  Eq.  3-58  for  a desired 
linear  operation  corresponding  to  h^°^  (•,♦); 

h(X!  5)  = 

= - bC,  ?) 

exp[j2iT(x  + C)f^]df^.  (3-60) 

The  mask  h(x;  C) / is  the  one  which  is  actually  utilized 
in  the  OSD  processor. 

As  with  the  OSD  processor,  an  alternate  interpre- 
tation can  be  ascribed  to  the  FOSD  processor.  Here,  we 
let 

k^°^  (fx  - v;  v)  = Kj^(fjj;  v)exp(-j2ircvf^) 

(3-61) 

where,  as  defined  in  Eq.  2-106,  K^(*r*)  is  the  Fourier 
transform  of  the  frequency  response  with  respect  to  x. 
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We  now  can  rewrite  Eq.  3-53  as 

“ /-»  F(v)k^°^(f^  - v;  v)dv.  (3-62) 

This  integral  relationship  is  akin  to  the  frequency  re- 
sponse characterization  in  Eq.  3-49.  Thus,  we  can  alter- 
nately interpret  the  output  of  the  FOSD  processor  as  cor- 
responding to  a frequency  response  characterization  of  a 
linear  system  with  frequency  response  k^°N*»*).  Note 
that  k^°^(*,»)  and  h^°^(*,*)  are  Fourier  transform  pairs: 

- v;  V)  = /“„  h^°Nfjj  - bC;  C)  exp  ( j2iTvC)  d? . 

(3-63) 

This  relationship  follows  from  Eqs.  3-59  and  3-62  and  Par- 
se val  ' s theorem . ^ 

An  obvious  question  now  arises.  Can  the  DOD  and  FDOD 
processors  be  alternately  interpreted  as  displaying  the 
output  spectrum  corresponding  to  some  linear  operation? 

The  answer  is  yes.  The  derivation  of  the  alternate  system 
line  spread  function  in  both  cases  is  straightforwardly 
derived  in  the  above  indicated  spirit.  Due  to  the  fact 
that  we  are  primarily  concerned  with  the  output  rather 
than  the  output  spectrum,  we  will  not  further  dwell  on 
these  alternate  interpretations. 
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3.6.  Lensless  Processing 

When  initially  presented,  it  was  noted  that  if  one 
is  interested  only  in  the  OSD  processor  output  intensity, 
it  was  possible  to  place  the  1-D  input  and  line  spread 
function  mask  directly  adjacent  the  spherical  lens.  Be- 
sides illiminating  vingnetting,  this  scheme  also  allows 
the  lens  transmittance  to  also  be  included  in  the  mask 
transmittance.  With  the  alternate  interpretation  of  the 
OSD  processor  presented  in  the  previous  section,  we  arrive 
at  the  rather  amazing  conclusion  that  a 1-D  space-variant 
system's  output  magnitude  can  be  displayed  utilizing  only 
an  input,  a mask,  and  a few  centimeters  of  free  space. 

Such  physical  characteristics  are  desireable,  of  course, 
when  one  is  primarily  concerned  with  processor  weight, 
physical  stability,  and  real  space  compactness.  It  is 
difficult,  in  fact,  to  imagine  a processor  in  a more  com- 
pact form.  In  this  section,  we  present  some  example  appli- 
cations and  experimental  results  for  two  such  1-D  lens- 
less  processors. 

3.6.1.  Unit  Magnif ier 

Even  though  it  is  doubtful  that  the  lensless  unit 
magnifier  will  find  wide  spread  use,  its  analysis  reveals 
some  rather  fascinating  aspects  of  lensless  1-D  process- 
ing. The  input-output  relationship  for  our  system  is  de- 
fined as 
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G(f^)  = (Xf)f[Xff^]  (3-64) 

where  Xf  corresponds  to  the  familiar  wave  length-focal 
length  product.  The  line  spread  function  corresponding 
to  this  operation  is 

h^°^(fjj  -bC;  C)  » (3-65) 

From  Eq.  3-58,  it  follows  that 

“ ^^^x  ■ |flexp[j2iTCf^]  . (3-66) 

Thus 

h(x;  C)  = expl  j2irCx/Xf  ]exp[ j2T7C^/^f ] • (3-67) 

Therefore,  we  can  write: 

h(x  - C;  C)  = exp[j2iTCx/Xf].  (3-68) 

The  mask  transmittance,  t(x;  C) / required  for  lensless 
processing  is  the  product  of  Eq.  3-68  with  the  transmittance 
of  a spherical  lens.  Choosing  the  lens'  focal  length  to 
be  f,  we  have 

t(x;  C)  = h(x  - C;  5)exp(^  (x^  + 

= exp[=^  (X  - C)^].  (3-69) 

Remarkably,  this  expression  is  recognized  as  the  trans- 
mittance of  a cylindrical  lens  with  focal  length  f^  = f/2 


rotated  45“  in  the  (x;  C)  plane.  The  mask  required  for 
unit  magnification  consists  of  a cylindrical  lens,  or 
equivalently,  a cylindrical  "hololens."  A hololens  is 
simply  a hologram  of  a lens  transmittance.^^ 

In  the  implementation  of  this  processor,  the  cylin- 
drical hololens  was  recorded  as  shown  in  Fig.  3-20 (a). 

The  cylindrical  lens  has  focal  length  f^  = 10cm  and  is 
illuminated  from  the  left  by  a coherent  plane  wave.  The 
photosensitive  medium,  as  shown,  is  placed  a distance  2f^ 
from  the  cylindrical  lens  and  is  illuminated  with  a coherent 
plane  wave  reference  beam  R.  The  reference  beam's  propa- 
gation direction  is  perpendicular  to  the  x axis  and  is 
incident  at  an  angle  of  0 = 30“ . 

The  lensless  1-D  unit  magnifier  is  shown  in  Fig. 

3-20 (b) . For  an  input,  we  use  the  double  pulse  shown  in 
Fig.  3-21 (a).  This  transmittance  is  placed  directly  in 
front  of  the  cylindrical  hololens  and  is  illuminated  from 
the  left  by  a coherent  plane  wave  propagating  in  the  con- 
jugate direction  to  that  of  the  reference  wave  used  in 
the  recording  geometry.  In  Fig.  3-21,  we  show  the  re- 
sulting field  amplitude  square  modulus  perpendicular  to 
the  z axis  at  various  values  of  z.  Figure  3-21 (a)  is  for 
z * 0 and  is  thus  the  double  pulse  input.  Figure  3-21 (b) 
is  the  field  eunplitude  square  modulus  for  a point  between 
z = 0 and  f^.  Note  that  a rotation  of  sorts  is  beginning 


Figure  3.21  Output  of  the  lensless  processor  in 

Figure  3.20b,  for  a)  z * 0 b)  z between 
0 and  f c)  z between  f_  and  2f_  and 

d)  z - 5f_.  *= 


119 


to  take  place.  The  field  amplitude  square  modulus  at  a 
point  between  z » and  2f^  is  shown  in  Fig.  3-21  (c). 
Lastly,  our  desired  processor  output  appears  at  z = 2f 
as  is  shown  in  Fig.  3-21 (d).  The  1-D  output,  taken  along 
the  vertical  axis,  as  predicted,  is  a unit  magnified  ver- 
sion of  our  1-D  input. 

3.6.2.  Lensless  Convolution 

For  this  example,  we  rewrite  the  convolution  integral 
of  Eq.  3-15  as 

- C/Xf)dC.  (3-70) 

Note  that  the  unit  magnifier  discussed  in  the  previous 
example  is  a special  case  of  this  relation  for  h^°^  (f^^  - 
C/Xf)  = ^(f^  - C/Xf) . 

To  determine  the  required  mask  for  this  processor, 
we  write  from  Eq.  3-58 

“x^^x'  " C/Xf)exp(j2TrCf^).  (3-71) 

From  this  relation,  we  can  easily  show  that 

h(x  - C;  5)  * H^°Nx)exp(i|j  xCJ  (3-72) 


where 


H^°Nx)  - /*„  h^®Nfjj)exp(j2Trf^x)df^  . 


(3-73) 


120 


To  find  the  transmittance  of  the  desired  mask  for  our  lens- 
less processor,  it  remains  to  multiply  Eq.  3-72  by  the 
transmittance  of  a spherical  lens  whose  focal  length  we 
shall  choose  to  be  f . The  result  is 


t(x;  C)  » h(x  - C;  C)exp[^(x^  + jj 

= (x)exp[^(x  - . (3-74)  | 

i- 

The  complex  exponential  term  is  again  recognized  as  a 
cylindrical  lens  of  focal  length  f^  = f/2  rotated  45“ 
on  the  (x,  C)  plane.  Thus,  1-D  convolution  can  be  per- 
formed by  placing  the  1-D  transmittances,  H^^^x)  and  f(C) 
back  to  back  adjacent  a rotated  cylindrical  lens.  The  ^ 

desired  processor  output  appears  at  a distance  2f  from 

c I 

this  optical  sandwich.  j 

For  an  experimental  implementation,  two  pairs  of  double 
pulses  [Fig.  3-22 (a)]  were  convolved  [Fig.  3-22 (b) ] . A ' 

hologram  was  made  of  the  1-D  Fourier  transform  of  the 

one-dimensional  transmittance  corresponding  to  h^°^ (C)  i 

using  a single  cylindrical  lens  with  a focal  length  of  20 
cm.  This  hologreun  served  as  the  transmittance  H^°^(f^) 
shown  in  the  processor  in  Fig.  3-23.  The  input  trans- 
mittance f(C)  was  placed  as  shown  adjacent  H^°Mf  ) and 

* fl 

a rotated  cylindrical  lens.  Since  a 20  cm  cylindrical  n 

lens  was  used  to  transform  h^°^  (C) r we  require  that  the 


1 


J 
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rotated  cylindrical  lens  have  half  that  focal  length. 

That  is, 

f^  = f/2  » 10cm.  (3-75) 

The  processor  was  illuminated  from  the  left  with  a co- 
herent plane  wave  with  wavelength  5145S.  The  processor 
output  appearing  a distance  2f  = 20cm  to  the  right  of 
the  rotated  lens  is  shown  in  Fig.  3-24.  (The  processor 
output  was  magnified  by  conventional  means  for  this  pic- 
ture.) Although  the  structure  of  the  vertical  axis  com- 
ponent of  this  picture  is  not  readily  apparent,  the  over- 
lapping portions  compare  quite  favorably  with  the  convolu- 
tion result  in  Fig. • 3-22 (b) . 

3.7.  Design  Cons iderations 

We  have  presented  four  basic  schemes  for  1-D  signal 
processing:  DOD,  OSD,  FDOD,  and  FOSD.  In  this  section, 
we  address  two  basic  questions:  1)  for  a given  space- 

variant  operation,  which  scheme  is  most  easily  implement- 
able  and  2)  if  a linear  operation  cannot  be  implemented 
in  its  obvious  form,  can  we  massage  xt  into  an  alternate 
impleroentable  form?  We  will  illustrate  these  design  chal- 
lenges respectively  with  the  one- dimensional  linear  oper- 
ations of  magnification  and  inverse  Abel  transformation. 
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3.7.1.  Magnifier  Design  Considerations 

The  DOD  processor  for  one-dimensional  magnification 
was  previously  presented  and  is  pictured  in  Fig.  3-3. 

The  line  spread  function 

h(x  - 5;  C)  = 6(x  - M?)  (3-76) 

is  generated  by  focusing  an  incident  plane  wave  to  a line 
source  which  in  turn  is  incident  on  the  1-D  input.  Let 
us  now  consider  three  possible  alternate  one-dimensional 
magnifier  processor  designs. 

Consider,  first,  the  FDOD  processor  for  the  magnifier. 
From  Eq.  2-22,  our  required  frequency  response  mask  must 
have  the  transmittance 

k(x  - cv;  v)  » 6 (x  - MC)  exp ( j2ir5v) d? 

* iexp(j2TTxv/M).  (3-77) 

Implementation  of  this  mask  is  possible  by  holographic 
means.  Compared  with  the  DOD  processor,  however,  this 
scheme  is  obviously  suboptimal  in  terras  of  ease  of  imple- 
mentation . 

Next,  we  inspect  the  OSD  and  FOSD  magnifier  processor 
for  which  we  must  restate  our  input-output  relationship 
as 


(3-78) 
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Thus 

- C;  C)  = 5(f^  - MC).  (3-79) 

From  Eq.  3-58,  it  follows  that 

H^(fx;  C)  = -SCfx  - MC)eifp(j2irCfx) 

= 6(fx  " MC)exp(j27rMC^).  (3-80) 

The  line  spread  function  mask  which  is  required  in  the  OSD 
processor  is  then 

h(xr  5)  = exp(J2itHC^)  - ME) 

= e3^”ME(E+x)  _ ,3.gj, 

Again,  our  mask  is  a complex  sinusoid  and,  again,  the  DOD 
processor  scheme  is  superior  as  far  as  ease  of  implemen- 
tation. 

Lastly,  we  consider  the  FOSD  magnification  processor. 
From  Eqs.  3-63  and  3-79: 

k^°Mx  - cv;  = /Zoo'^^^x  ”MC)exp(j2iTvC)dC 

= ^ exp  (j2irfxV/M) . (3-82) 

From  Eq.  3-61: 

*Here,  b = 1 [see  Eq.  3-58] 


1 


127 

^ S exp[j2uf^v(c  + ^)].  (3-83) 

Inverse  transforming  gives  the  required  mask  transmittance: 

k(x;  V)  = J2„  '')®xp(j2Tif^x)df^ 

= ^[x  + (c  + |)].  (3-84) 

Here,  as  in  the  DOD  processor,  our  "mask"  is  a Dirac  delta 
sheet.  In  the  FOSD  processor,  however,  we  cannot  use  a 
cylindrical  lens  to  generate  the  Dirac  delta  sheet  unless 
we  use  the  input  spectrum  P(v),  rather  than  f(C)»  as  the 
processor  input.  Under  the  assumption  that  we  do,  in 
fact,  favor  f(C)  as  our  processor  input,  the  DOD  scheme 
is  again  superior. 

The  magnifier,  as  we  have  seen,  is  a relatively  easy 
linear  system  to  evaluate.  In  general,  processor  analysis 
in  the  manner  just  presented  becomes  rather  analytically 
intractable.  The  potential  use  of  such  an  approach, 
however,  is  obvious. 

3.7.2.  Inverse  Abel  Transform 

As  has  been  shown,  the  Laplace  and  Fourier  transform 
can  both  be  relatively  straightforwardly  evaluated  by  co- 
herent optical  means.  The  Mellin  transform  can  also  be 
49-50 

thus  evaluated.  All  integral  transforms  are,  in  fact, 

special  cases  of  the  superposition  integral.  The  kernels 


(line  spread  functions)  corresponding  to  many  such  trans- 
forms, however,  lie  outside  the  unit  circle  on  the  complex 
plane  and  thus  cannot  be  directly  represented  by  a hologram. 
An  example  of  such  a transform  is  the  inverse  Abel  trans- 
form defined  by;^^'^® 

-4 /"xi  a! 

(3-85) 

where  y(*)  is  the  unit  step  function.  The  transform  kernel 
is 

h(x  - C;  C)  = Zi  HiizMl  -|^(•).  (3-86) 

The  kernel's  dynamic  range  extends  from  zero  to  infinity. 

This,  and  the  fact  that  the  kernel  contains  a derivative 

operator,  precludes  holographic  representation.  Thus, 

the  DOD  processor  cannot  be  directly  utilized. 

Let  us  then  investigate  the  required  transmittance 

28 

for  the  FOOD  processor.  We  utilize  Parseval's  theorem 
and  rewrite  Eq.  3-86  as 

9(X)  - /Zi4 


X [-j2TrvP  (v)  ]dv 


(3-87) 
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where  we  have  recognized  that 

* -j27rvF(v).  (3-88) 

Our  frequency  response  here  is  thus  given  by 

k(x  - cv,  V)  - 12  r ,3-89) 

1*1  • 

where  the  derivative  operation  is  here  included  in  the 
multiplicative  v term.  From  Gradshteyn  and  Ryzhik^^ 

[Eqs.  3-771-7  and  3-r771-9] , this  integral  reduces  to 

k(x  - cv;  v)  = jirv  [Jq  (2Trvx) -jN^  (2irvx)  ] (3-90) 

where  J (•)  and  N_(*)  are  zeroth  order  Bessel  functions  of 
o o 

the  first  and  second  kind  respectively.  Although  vN^(v) 

52 

is  well  behaved  about  the  origin,  is  not  ),  we 

would  prefer  a mask  that  is  strictly  real.  Equation  3-90 
obviously  does  not  conform  to  this  preference. 

It  is  possible  to  utilize  a variation  of  the  DOD  pro- 
cessor to  compute  the  inverse  Abel  transform  where  the 
required  mask  transmittance  is  purely  real.  This  proces- 
sor, pictured  in  Fig.  3-25,  has  a transmittance  of  J^(2ttvx) 
placed  in  plane  PI.  The  three  cylindrical  lenses,  Ll,  L2 
and  L3,  then  perform  a one-dimensional  Fourier  transform 
with  respect  to  v.  Thus  the  field  amplitude 

iZca  Jq  (2iTxv)exp(-j2'ir5v)  dv 


(3-91) 
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is  incident  on  plane  P2.  Since  Jq(*)  is  an  even  function, 
we  can  rewrite  Eq.  3-91  as 


(2itxv)  cos  (2irCv)  dv 

•w  O 


= 2/~  Jq(2tt|x1  v)cos(2TTCv)dv 


; 5 < |x| 
} ^ > 1x1 


(3-92) 


where  we  have  used  Eq.  6.671.2  of  Gradsteyn  and  Ryzhik.^^ 

The  values  of  the  function  J_(*)  lie  within  the  unit  circle 

o 

on  the  complex  plane  and  thus,  in  principle,  can  be  re- 
presented holographically.  The  form  of  the  relation  in 
Eq.  3-92  looks  similar  to  that  in  Eq.  3-86  except  that 
the  polarities  are  opposite.  That  is,  Eq.  3-84  requires 
that  C > 1*1  while  Eq.  3-92  is  nonzero  only  for  C < lx]. 

We  can,  however,  rectify  this  situation  by  performing  the 

A, 

variable  substitution  5 = -5  in  Eq.  3-85  to  give: 


g(x)  = ? /!»  f (3-93) 

/ x^- 

Here,  our  integration  interval  covers  the  desired  range 
of  C < 1 X I . 

Returning  now  to  our  processor,  we  generate 
electronically  with  the  op-amp  differentiator  shown.  This 
is  fed  into  a one-dimensional  electro-optical  transducer 
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whose  output  transmittance  multiplies  Eq.  3-92.  The  field 
amplitude  immediately  to  the  right  of  plane  P2  is  thus 


f 


1 


1 


IT 


y(-C-lxl) 


d 


(3-94) 


This  product  is  Fourier  transformed  with  respect  to  C by 
lenses  L4,  L5  and  L6  so  that  incident  on  plane  P3  is  the 

[ field  amplitude: 

I lar  f(-5)le-J2"''«d5.  (3-95) 

For  V = 0»  this  relation  reduces  to  that  in  Eq.  3-93.  Our 
desired  inverse  Abel  transform  corresponding  to  an  input, 
f(C)/  thus  appears  in  plane  P3  along  the  x axis. 

The  inverse  Abel  transform  processor  cannot  be  clas- 
sified into  any  of  the  processor  classes  thusfar  consid- 
ered. The  ideas  leading  to  its  design,  however,  parallel 
closely  those  considered  previously.  We  thus  conclude 
that,  given  a space-variant  operation,  attention  should 
not  be  restricted  only  to  the  DOD,  OSD,  FDOD,  and  FOSD 
processors,  since  augmentation  and  perturbations  of  funda- 
mental underlying  concepts  can  be  utilized  to  yield  a de- 

t • 

\ sign  which  is  more  highly  implementable. 

I 

► 

I 3.8.  Conclusion 

I In  this  chapter,  we  have  presented  n unerous  coherent 

i optical  processors  that  are  capable  of  performing  a wide 
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class  of  one-dimensional  operations.  All  of  the  proces- 
sors require  a mask,  possibly  holographic  if  the  kernel 
is  bipolar  or  complex,  and  a one-dimensional  input  trans- 
mittance. The  E>OD  and  OSD  processors  utilize  a line  spread 
function  mask  and  respectively  display  the  output  and  out- 
put spectrvim.  The  FOOD  and  FOSD  processors  do  likewise, 
but  utilize  a frequency  response  mask.  In  all  of  these  pro- 
cessors, the  1-D  system  output  is  viewed  along  a line  in 
the  processor  output  plane.  The  Laplace  transform  and 
aunbiguity  function  operations  were  shown  to  be  two  special 
cases  where  the  entire  output  plane  could  be  utilized. 

The  OSD  and  FOSD  processors  were  shown  to  be  able  to 
be  interpreted  as  displaying  the  output  (instead  of  output 
spectrum)  corresponding  to  a linear  operation  that  could 
be  deduced  from  the  system  line  spread  function.  Along 
with  the  DOD  and  FDOD  processors,  there  are  thus  four 
general  schemes  for  direct  one- dimensional  processing 
each  of  which  require  different  masks.  For  a given  linear 
operation,  we  can  thus,  in  principle,  derive  four  separate 
expressions  for  possible  masks  for  use  in  an  appropriate 
1-D  processor.  This  allows  for  a degree  of  flexibility 
in  processor  design.  Using  these  alternative  schemes  in 
design  considerations  was  exemplified  in  the  case  of  the 
1-D  magnifier.  As  was  shown  for  the  case  of  the  inverse 
Abel  transformer,  however,  attention  should  not  be  re- 
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4.  TWO-DIMENSIONAL  SPACE-VARIANT  PROCESSORS 

In  the  previous  chapter,  the  two-dimensional  nature 
of  coherent  optical  processors  was  utilized  to  perform  one- 
dimensional operations.  This  chapter,  on  the  other  hand, 
describes  various  possible  schemes  for  performing  general 
two-dimensional  operations  with  coherent  optical  processors. 

There  are  a number  of  specific  2-D  linear  space-var- 
iant operations  that  can  be  performed  by  coherent  proces- 
sors. The  thin  lens,  for  example,  performs  a 2-D  Fourier 
transformation  which  is  a linear  space-variant  operation. 

The  astigmatic  processors  encountered  in  the  previous 
chapter,  wherein  orthogonal  Fourier  transformation  and 
imaging  are  performed,  are  also  space-variant . The  oper- 
ations of  Laplace  transform  and  ambiguity  function  display 
can  be  considered  as  special  cases  of  two-dimensional 
space-variant  astigmatic  processing.  The  elementary  act 
of  non-unity  magnification  imaging,  also,  must  rigorously 
be  considered  space-variant.^^  Besides  an  input  trans- 
mittance, the  Fourier  transform,  astigmatic,  and  imaging 
processors  consist  only  of  two  basic  elements:  lenses 
emd  free  space.  Positioning  of  lenses  in  valcious  other 
locations  will  of  course  result  in  other  space-variant 
operations.  The  class  of  possible  operations,  however,  is 
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^ somewhat  limited.  Rather  general  treatment  of  such  multi- 

S 54 

; plane  optical  processors  is  given  by  Carlson. 

f 

[ An  alternate  treatment  of  generalized  2-D  space-var- 

55-57 

I iant  processing  has  been  developed  by  Bryndahl  et.al. 

Inherent  in  his  development  is  the  assumption  of  a "slowly 

f varying"  input.  Both  Bryndahl 's  and  Carlson's  methods, 

I as  well  as  other  space-variant  processing  schemes,  are 

58 

discussed  by  Goodman. 

Our  treatment  of  the  design  of  space-variant  processors, 

I , on  the  other  hand,  will  be  considered  from  the  systems  view- 

|i  point.  In  1-D  processor  design,  the  continuum  orthonormal 

I basis  set  response  characterization  of  linear  systems  was 

Ir  used.  In  2-D  processor  design,  we  will  use  the  character- 

izations developed  in  Chapter  II  where  limiting  assump- 
tions or  approximations  were  made  concerning  the  system 
and/or  the  class  of  allowable  inputs.  These  include  the 
piecewise  isoplanatic  approximation  (PIA) , discrete  ortho- 
\ normal  basis  set  response  (DBR) , and  sampling  theorem  char- 

f 

acterizations . 

In  the  treatment  of  2-D  processors,  two  views  are 
I possible.  The  first  is  that  we  have  an  actual  physical 

I coherent  processor  which,  for  reasons  of  real  space  com- 

pactness and  stability,  we  wish  to  characterize  holo- 

I 

^ graphically.  For  this  case,  we  need  to  physically  probe 
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responding  system  responses  must  then  somehow  be  holographi- 
^ cally  cataloged  in  a straightforwardly  accessable  manner. 

The  volume  hologram  approach  to  be  discussed  is  a scheme 
based  on  this  view. 

The  second  manner  in  which  2-D  space-variant  design 
can  be  viewed  is  similar  to  the  1-D  processor  concept. 

That  is,  instead  of  optically  recording  a required  mask, 

the  mask  transmittance  is  generated  by  alternate  means 

159 

such  as  computer  generated  holograms.  All  schemes  in 

\ 

^ this  chapter,  except  the  volume  hologram  treatment,  can 

I be  implemented,  in  principle,  by  either  the  direct  re- 

f 

[ cording  method  or  the  computer  generated  mask  methods. 

Three  basic  schemes  for  general  two-dimensional 
space-variant  processing  will  be  considered.  The  first 

: scheme  involves  composite  holograims  and  is  a direct  appli- 

i 

cation  of  a space-variant  system  sampling  theorem  developed 
in  Chapter  II.  The  constraints  on  system  input  and  line 

; 

spread  here  impose  severe  limitations  on  allowable  opera- 
tions. The  second  scheme  utilizes  the  volume  hologreum  in 

which  system  responses  are  angle  multiplexed  within  a thick 

19 

medium.  This  scheme  was  first  suggested  by  Burton,  et.al. 

\ Application  of  the  PIA  characterization  to  the  volume 

r 

hologram  was  investigated  by  Marks.®  The  sampling  theorem 

20-21 

' approach  was  similarly  applied  by  Deen  et.al.  Other 

1 results  have  been  reported  by  Walkup  and  Hagler®*^  and  Marks 

8 


2 


et.al 
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The  third  method  of  two-dimensional  space-variant  pro- 

22 

cessing,  proposed  by  Krile  et.al.,  makes  use  of  phase 
coded  reference  beams.  An  effect  inherent  in  this  method 
is  the  smearing  of  unwanted  crosstalk  terms  into  diffuse 
background  noise. 

Since  this  chapter  deals  with  two-dimensional  space- 
variant  processing,  all  previous  linear  system  character- 
izations, which  had  a 1-D  favor,  must  be  generalized  to 
two  dimensions.  For  example,  the  2-D  superposition  inte- 
gral is  written  as 

g(x,y)  = /“^  /"^f (C,n)h(x-c,y-n; C,n)dCdn . 

(4-1) 

The  kernel,  h (*,•;*,•),  for  obvious  reasons,  is  now  re- 
ferred to  as  the  system  point  spread  function.  As  will 
become  clear  in  the  development  to  follow,  mathematical 
and  conceptual  generalization  from  one  to  two-dimensional 
treatment  is  trivally  straightforward. 

We  should  point  out  that  the  line  spread  function 

does  have  meaning  in  the  context  of  2-D  system  character- 
61—62 

ization.  In  this  chapter,  however,  we  will  concern 

ourselves  only  with  the  2-D  system's  point  spread  func- 
tion. 


4.1.  System  Response  Generation 

In  the  system  characterizations  developed  in  Chapter 
II,  it  was  necessary  to  catalog  the  system  responses  to  a 
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class  of  inputs  which  were  specified  by  a limiting  assump- 
tion placed  on  the  allowable  input  class  and/or  system. 

Both  the  piecewise  isoplanatic  approximation  (PIA)  and 
the  seunpling  theorem  characterizations  require  knowledge 
of  the  system  response  to  Dirac  deltas  placed  at  certain 
points  at  the  system  input.  An  alternate  sampling  theorem 
required  scunpling  be  performed  at  the  system  output.  The 
Fourier  duals  of  the  sampling  theorems  require  that  the 
frequency  response  be  sampled  in  either  the  input  or  output 
domain.  For  the  discrete  orthonormal  basis  set  response 
(DBR)  characterization,  we  need  to  catalog  the  system  re- 
sponse to  inputs  consisting  of  each  element  in  the  basis 
set. 

In  this  section,  we  explore  physically  generating 
these  responses.  That  is,  given  an  actual  linear  coherent 
optical  system  with  input  and  output  planes,  how  do  we 
physically  generate  the  required  catalog  of  system  re- 
sponses? Once  the  desired  responses  are  generated,  it  re- 
mains only  to  find  a method  by  which  to  optically  catalog 
the  response  (on  a hologram  for  example)  in  such  a manner 
as  to  be  retrievable  for  system  representation.  These 
methods  include  the  volxime  hologremi  and  phase  coded  ref- 
erence beam  schemes  which  will  be  presented  shortly.  We 
present  the  response  generation  techniques  in  a separate 
section  for  the  simple  reason  that  certain  cataloging  methods 


t 


140 


can,  in  principle,  utilize  more  than  one  of  the  re- 
sponse characterizations  in  their  implementation.  We 
should  also  keep  in  mind  that,  except  for  the  volume  holo- 
gram approach,  the  required  system  response  masks  can  either 
be  computer  generated  or  optically  recorded. 


4.1.1.  Impulse  Response 

In  the  scunpling  theorem  and  PIA  system  characteriza- 
tions, we  require  generation  of  the  system  point  spread 
function.  For  an  arbitrary  coherent  optical  system,  this 
can  be  performed  as  shown  in  Fig.  4-1.  An  incident  co- 
herent plane  wave  is  focused  to  the  point  (C»n)  on  the 
system,  input  plane.  Since  the  spherical  lens  is  circular 

with  radius  A,  this  "point"  is  better  described  by  the 

23 

Fourier  transform  of  a circle  function: 


where 


F [circ  ( 


A 


-)i  =1 


J3^(2TrAp) 


(4-2) 


X f 


(4-3) 


circ  r A 


1 

0 


; |r|<  1 
; I r I > 1 . 


(4-4) 


and  Jj^(»)  is  the  first  order  Bessel  function  of  the  first 

23 

kind . One  possible  expression  for  2-D  Dirac  delta  is 
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«(x^,  y^)  = lim  I Jj^(2TrAp).  (4-5) 

Thus,  we  conclude  that  for  a sufficiently  large  lens  radius, 
A,  our  field  amplitude  on  the  system  input  plane  is  well 
approximated  by  a Dirac  delta.  We  can  move  the  location 
of  the  delta  to  the  coordinates  (x^,y^)  = (C»n)  simply 
through  an  appropriate  lens  translation.  The  correspond- 
ing field  amplitude  distribution  on  the  system  output  plane 
is  then  the  system  point  spread  function  corresponding  to 
the  input  point  (C,ri).  Physically,  this  corresponds  to  a 
sample  point  spread  function.  We  stress  this  by  letting 
the  coordinates  (C,ti)  correspond  to  the  specific  point 
. Then  the  field  amplitude  distribution  on  the 
system  output  plane  is 

^'(^o-^n'yo-^'-^n'V.  <4-6) 

In  this  expression,  x^  and  y^  are  our  spatial  variables 

while  and  can  be  considered  as  indexes.  Thus,  we 
n m 

can  generate  the  sample  point  spread  function  h(x  ,y  ,n  ) 

o o n m 

by  simply  translating  the  origin  of  the  (x^/Yq)  plane  to 

the  point  (C„»n„). 

n m 

Suppose,  now,  we  wish  to  generate  forms  of  the  system 
point  spread  function  sampled  in  x and  y rather  than  in  C 
and  n.  Consider,  again,  Pig.  4-1  with  attention  restricted 
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this  time  to  the  point  (^qiYq)  ~ output  plane. 

As  the  input  point  source  explores  various  values  of  (C,n)/ 
the  corresponding  field  amplitude  at  the  output  point  (x  ,y  ) 
is  h(x  -?,y  -n; C»n)  . 

To  holographically  record  a transmittance  of  the  out- 
put sampled  point  spread  function,  we  must  translate  a 
photosensitive  medium  placed  in  the  (x^,y^)  plane.  Imagine 
first,  a point  source  input  at  (C»n)  = (0,0)  and  the  pho- 
tosensitive medium  in  the  (XQ,yQ)  plane.  Immediately 
in  front  of  the  mediiim  we  place  an  opaque  mask  that  passes 
light  only  at  the  point  (x  ,y  ) » * Then,  as  though 

rigidly  connected,  the  input  point  source  and  photosensi- 
tive medium  translate  through  all  values  of  Dur- 

ing this  translation,  the  opaque  mask  remains  stationary. 

Some  reflection  on  the  reader's  part  reveals  that  the 
transmittance  recorded  utilizing  this  procedure  is  h(Xp, 
yq;C,n) . 


4.1.2.  Frequency  Response  Sampling 

As  with  the  impulse  response,  we  wish  to  formulate 
methods  by  which  the  Scunple  frequency  responses  of  a co- 
; herent  processor  can  be  directly  recorded  both  in  the  in- 

put and  output  domains.  In  two  dimensions,  the  frequency 
response  is  defined  as 


k[Xo-Cx''^»yo-Cy'^riJ 

*S[exp{j2Tr(  v^x^  + 


(4-7) 


Li 


where  the  constants  c and  c retain  dimensional  consis- 

X y 

tency.  In  coherent  optics,  the  complex  exponential  re- 
quired for  our  system  input  is  interpreted  as  a tilted 
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plane  wave.  This  can  probably  be  most  easily  seen  from 
the  geometry  of  Pig.  4-2. 

A Dirac  delta  (point  source)  is  placed  in  the  (a. 6) 
plane  at  the  coordinates  (a,b) . The  thin  lens  performs 
a Fourier  transform  so  that  the  field  eunplitude  incident 
on  the  (x^,y^)  plane  is 

exp[“^£—(Xj^^a  + y^b)  ] . (4-8) 


From  geometrical  optics,  we  see  that  Eq.  4-8  describes  a 

63 

tilted  plane  wave  incident  on  the  (x^,y^)  plane.  Compar- 
ing with  the  argument  of  the  system  operator  in  Eq.  4-7, 
we  have 


(4-9) 


In  an  actual  laboratory  situation,  generating  tilted 

plane  waves  in  the  manner  shown  in  Fig.  4-2  would  probably 

not  be  used.  Rather,  we  would  generate  the  plane  waves  in 

accordance  with  the  direction  cosines,  cos6  and  cos0  . 

X y 

Prom  the  geometry  in  Fig.  4-2,  these  are  given  by 


COS0  » -a/r  ; cos0  = -b/r  (4-10) 
X y 
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where 

(4-11) 

In  general,  a and  b will  be  much  smaller  than  f so  that 
r ::  f.  Using  this  approximation  and  Eq.  4-9,  we  can  then 
rewrite  Eq.  4-10  as 

0 = cos”^  Xv  ; 0 = cos~^  Xv  (4-12) 

X ^ y n • 

These  are  the  desired  plane  wave  propagation  angles  cor- 
responding to  the  input  frequencies  and  v^. 

To  generate  a coherent  optical  system's  frequency  re- 
sponse, we  simply  feed  a tilted  plane  wave  into  the  system 

input.  If  the  input  frequencies  are  and  v then  the 

’n  in 

corresponding  field  cunplitude  distribution  on  the  system 

output  plane  is  the  sample  frequency  response  given  by 

k(x  - c V , y_  - c„v  ; V , V ).  This  corresponds 
° * *^n  ° y 'm  ^n  'm 

to  input  Scunpling. 

For  output  sampling,  we  can  utilize  a scheme  similar 
to  that  used  for  output  sampling  of  the  system  point  spread 
function.  We  first  choose  our  sample  coordinate  in  the 
output  plane,  • Then,  as  we  sweep  the  various  fre- 

quencies at  the  system  input,  our  output  photosensitive 
medium  translates  accordingly.  For  input  frequencies  of 
and  our  photosensitive  mediums  origin  is  located, 

in  the  output  plane  at  coordinates  (c^^v  ,c  ) . in  this 

^n  y m 


w 


147 


manner,  we  generated  the  sample  frequency  response 


4.1.3.  Sine  and  Rect  Response  Generation 

In  the  development  of  the  discrete  orthonormal  element 
basis  set  response  (DBR)  characterization  of  linear  systems 
in  Chapter  II,  we  investigated  specifically  the  sine  and 
rect  orthonormal  basis  sets.  Here,  we  explore  the  gener- 
ation of  these  responses  in  coherent  optical  systems. 

First  consider  the  two-dimensional  (n,m)^^  element 
of  the  rect  basis  set  here  generalized  from  Eq.  2-62: 


2/^ 


exptj2TT(v  X + V y )] 
^n  m 


X.  y , 

X rect  rect  (4-13) 

X y 

where  v » “ ni/2T  , and  where  we  have  restrict- 

ed  our  input  to  be  identically  zero  for  [x^l  > and  jy^l 
> Ty.  The  form  of  Eq.  4-13  is  identical  to  that  of  the 
input  Scunpled  frequency  response  except  for  the  rect  terms. 
Thus,  we  can  generate  the  system  rect  response  in  the  seune 
manner  as  we  generated  the  system  frequency  response:  by 
probing  the  system  input  with  tilted  plane  waves.  The 
rect  terms  can  be  included  by  simply  making  the  input  plane 
opaque  outside  of  the  rectangular  region  |x^|  £ T^  and 
ly^l  ^ Ty.  Prom  Eq.  4-12,  the  angular  propagation  direc- 


148 


tions  for  the  input  plane  wave  are 

6 = cos  ^ Xvp  = cos  ^ 5^ 

*n  ^n 

0„  = cos"^  Xv  = cos"^  5^  . (4-14) 

2Ty 

The  system  sine  response  is  also  straightforwardly 
generated.  The  two-dimensional  (n,m)^^  element  of  the  sine 
basis  set  is  generalized  from  Eq.  2-57  as 

sine  2W^(x^-x^)  sine  2Wy(y^-yj^)  (4-15) 

where  2W  and  2W  are  the  allowable  bandwidths  of  the  in- 
X y 

put  in  X.  and  y.  respectively  and  x = n/2W  and  y_  = m/2W„. 

1 1 ■*  n X m y 

The  sine  basis  set  elements  can  be  generated  as  shown 
in  Fig.  4-3.  A rectangular  pupil  is  placed  on  the  (a, 8) 
plane  and  is  described  by 

X y 

where  f is  the  focal  length  of  the  lens.  On  the  (x^,y^) 
plane  appears  the  Fourier  transform  of  this  pupil  evalu- 
ated at  the  frequency  components  (x^/Xf,  y^/Xf)  and  trans- 
lated to  the  coordinates  (x  ,y  ) . Outside  of  the  propor- 

n m 

tionality  terms / this  field  amplitude  is  identical  to  Eq. 
4-15.  We  can  generate  other  elements  in  the  sine  basis 
set  simply  by  translating  the  rectangular  pupil  and  thin 
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lens  while  leaving  the  (x^,y^)  plane  stationary.  These 
elements  serves  as  the  input  to  our  system  and  the  cor- 
responding field  amplitude  at  the  system  output  is  the 
desired  system  sine  response. 

4.1.4.  Discussion 

The  system  characterizations  considered  in  Chapter  II 
all  utilized  the  principle  of  probing  the  input  with  var- 
ious signals  and  cataloging  the  corresponding  system  fre- 
quency response.  As  has  been  shown,  this  response  gener- 
ation can  be  performed  straightforwardly  in  coherent  opti- 
cal systems.  It  remains  to  develop  some  schemes  for  cata- 
loging. This  is  the  topic  of  the  remainder  of  this  chap- 
ter. 


4.2.  A Composite  Hologram  Approach 

In  Chapter  II,  we  developed  a sampling  theorem  for 
space-variant  systems  whose  line  spread  function  was  band- 
limited  in  X.  From  Eq.  2-73,  the  two-dimensional  gener- 
alization of  this  expansion  is 

g(3C/y)  * I I [f  (x,y)h(x^,yjj^;x,y)  ] 
n m 

*[sinc2Wj^(x-x^)sinc  2Wy(y-yjj^)]  (4-17) 

where  x = n/2W  ,y  = m/2W  , and  2W  and  2W  are  the  (low- 

n X Ml  j X y 

pass)  bandwidths  of  h(x,y;5,n)  respectively  in  the  x and 
y directions.  This  relationship  can  alternately  be  writ- 


ten as 
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L 


g(x,y)  - I Il*(x-Xn.y-y„)h(x„,y„jx-x^,y-y„)l 


n m 


*[sinc  (2W  x) sine (2W  y) ] . 
X y 


(4-18) 


We  can  interpret  Eq.  4-18  as  follows:  The  input, 

f(x,y)  is  shifted  to  the  coordinates  (Xj^»yjj,)  on  the  (x,y) 

plane  whereupon  it  is  multiplied  by  a corresponding  sample 

point  spread  function.  This  process  is  repeated  for  all 

sample  coordinate  pairs,  (Xjj,yjj|)  . All  resulting  products 

are  then  s\immed  and  passed  through  an  appropriate  2-D  low 

pass  filter  which  gives  rise  to  the  required  convolving 

sines  in  Eq.  4-18.  Note  that  no  assumptions  concerning 

the  input  have  yet  been  made. 

Suppose,  now,  we  wish  to  simultaneously  perform  all 

shifted  samplerinput  point  spread  function  multiplication 

on  a single  plane.  In  order  to  assure  that  no  overlapping 

of  terms  occurs,  (i.e.,  no  crosstalk),  we  must  require  for 

each  n and  m that  the  product,  f(x-x„,y-y„)  h (x^,y_;x-x_, 

n Id  n In  n 

y-y  ) , is  identically  zero  outside  the  rectangle  defined 
m 

by 


|x  - Xn>  - \ 

ly  “ yjnl  1 1/4  Wy 


(4-19) 


or,  equivalently 
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2n-l 

4W 


< X < 


2n-H 
4W . 


(4-20) 


2in-l 

4W„ 


2m+l 
4W  . 


This  can  occur  in  two  ways:  1)  Our  input  is  identically 

zero  for  |C|  < 1/4  W and  In  | < 1/4  W or  2)  the  point 

X y 

spread  function  h(x,y;5fn)  is  identically  zero  outside 

the  interval  |C|  < 1/4  W and  |n|  < 1/4  W . Thus,  in 

X y 

order  to  perform  all  required  multiplication  on  a single 
plane,  severe  compact  support  constraints  must  either  be 
placed  on  either  input  or  point  spread  function.  This 
is  in  addition  to  the  constraint  that  h(x,y;C»n)  be  band- 
limited  in  both  X and  y. 

A coherent  processor  capable  of  carrying  out  the 
Seunpling  theorem  in  Eg.  4-3  under  the  discussed  constraints 
is  pictured  in  Fig.  4-4.  In  plane  PI,  we  place  our  input, 
f(C,n).  If  the  compact  support  constraint  falls  on  the 
line  spread  function  rather  than  the  input,  it  is  neces- 
sary to  mask  that  region  of  plane  lying  outside  the 

rectangle  |5l  < 1/4  W , (til  < 1/4  W . This  is  without 

X y 

loss  of  information  since  the  system  we  are  characterizing 
sees  the  input  only  over  these  intervals. 

Spherical  lens,  Ll,  performs  a two-dimensional  Four- 
ier transformation  on  the  input.  The  field  eunplitude 
incident  on  plane  P2  is  thus 
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F(v^,v^)  = If(5,p)l 


(4-21) 


* lZjZa>  f (Crn  )exp[-j2TT(v^C-K)^n  ) ]d5<3n  , 

Placed  in  plane  P2  is  a pin  hole  mask  the  transmittance 
of  which  can  be  modeled  as 


t (v  ^ ,Vn  ) = I (v  ^-2nW^)  6 (v^  -2mWy) 

’ 4wV  coi«b(^)coinb(^) 

X y X y 

where  the  comb  function  is  defined  by 


(4-22) 


comb(v)  4 I<S(v-n) 
” n 


(4-23) 


The  field  amplitude  immediately  to  the  right  of  P2  is  then 
given  by 


IwV  ^ ^ ^ <2W“)  • 

V ” V xr 


(4-24) 


This  product  is  again  Fourier  transformed  by  lens  L2  to 
give  incident  on  plane  P3  the  field  amplitude 


,,  tF(v  )]*F  w )comb(~^j  )] 


5'  ti'J 


(4-25) 


where 


“ /r«on«,<*)®*Pl:j2TT(v^x  + v^y)]dv^dv^  . 


(4-26) 
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Using  the  scaling  theorem  and  the  fact  that 

F^[comb(C)]  = comb(v)  (4-27) 

reduces  Eq.  4-25  to 

f (x,y)* [comb(2W  x)comb(2W  y) ] 

* y 

X y n m 

where  x^  = n/2W„  and  y„  = m/2»<  . Our  input  is  thus  re- 
n X m y 

plicated  over  plane  P3.  Due  to  the  compact  support  of 
the  input  and  our  choice  of  spacing  in  the  pin  hole  mask, 
these  replications  do  not  overlap.  We  stress  this  by  re- 
writing Eg.  4-28  as 

I l£(x-x„,y.y„, 

rect  [2W^(x-x^)  Jrect  [2Wy(y-yjjj)  ] . (4-29) 

Placed  in  plane  P3  is  a mask  of  sample  point  spread  func- 
tions. The  transmittance  of  this  mask  is 

n n 

rect  (2  W^(x-x^)  Irect  [2Wy  (y-y^jj)  ] . (4-30) 

If  the  transmittance  is  not  non-negative  real,  we  require 
a holographic  mask.  The  placement  of  a number  of  holo- 
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grams  side  by  side  on  a plane  results  in  what  is  termed 

64 

a composite  hologrcim.  The  mask  can  either  be  directly 
recorded  or  a computer  generated  transmittance.  The  mask 
is  multiplied  by  Eq.  4-29  to  give  immediately  to  the  right 
of  plane  P3  the  field  amplitude 


(4-31) 


where  we  have  now  suppressed  the  rect  terms.  Comparing 
this  with  Eq.  4-18,  we  see  that  now  only  low  pass  filtering 
is  required.  This  is  performed  by  lenses  L4  and  L5  and 
the  rectangular  pupil  on  plane  P4 . Our  desired  output 
then  appears  on  plane  P5  and  our  space-variant  operation 
has  been  performed. 

Again,  the  class  of  space-variant  operations  that 
can  be  performed  using  this  "composite  hologram  approach" 
is  somewhat  limited.  Our  constraint  is  that  either  the 
point  spread  function  in  C and  n and/or  the  input  have 
compact  support  not  exceeding  the  reciprical  bandwidths 
of  h(x,y;C,n)  in  x and  y respectively. 

To  illustrate  application  of  this  processor  to  a 
specific  operation,  we  will  now  consider  a Fourier  trans- 
form operation  defined  by 


g(x,y)  = ^bf(C/n) 

exp[-j2Trr  (x?  + yti)]dCdri 


(4-32) 
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where  a and  b are  positive  spatial  constants  and  r,  with 
dimension  of  frequency  squared,  retains  the  required  di- 
mensionless nature  of  the  exponent.  The  point  spread 
function  corresponding  to  this  operation  is 

h(x-C,y-ri;C/n)  = exp  [- j2iTr  (x?  + yn)] 

rect(^)rect(5^).  (4-33) 

It  follows  that 

h(x,y;C,ri)  = exp [-j2irr{  (xC  + yn ) + + n^)] 

rect  (|^)  rect  (^) . (4-34) 

The  two-dimensional  transfer  function  for  this  system  is 

then 

Hxy  (fx'^y' ^ j^y[h(x,y;?,n)  ] 

= 6(f  -r5)6(f  -rn) 

X y 

exp[-j2iTr  (C^  + n^)  ] 
rect(^)rect(^)  . (4-35) 

Due  to  the  multiplicative  Dirac  delta  terms,  we  can  malce 

the  substitutions  C * This  gives 

X y 

Hxy (fx'^y'^'^^  * 6 (f^-rC)6 (fy-rn) 

* £/)lrect(^)rect(^). 


exp[- 


(4-36) 


Clearly,  the  point  spread  function  is  bandliinited  in  both 
X and  y with 
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. ra 

Wy  = rb.  (4-37) 

Note  also  from  Eq.  4-34,  that  h{.x,y;^,r\)  is  zero  outside 
of  the  interval  |Cl  £ 1/4  W^,  [nl  < 1/4  if  we  impose 
the  inequalities 


Inl  ib  < jjj 


or,  equivalently,  since  r > 0, 


(4-38) 


a 


(4-39) 


If  these  inequalities  are  assumed,  no  compact  support 
assvimptions  for  the  input  need  be  made. 

To  this  point,  we  have  only  established  that  the 
sampling  theorem  in  Eq.  4-33  can  be  applied.  The  next 
step  is  to  compute  the  required  transmittance  of  the  sample 
point  spread  function  mask.  From  Eq.  4-34,  we  can  write 


h(Xn#yjjj;x,y)  = exp[-j2Trr{x(x  + x^)  + y(y  + yj^)}] 
rect(^)rect(-^) , 


(4-40) 
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It  follows  that 


h(Xn,y^;x-Xn,y-y^) 

= exp[-j2Trr{  (x-x^)x  + (y-yj^)y}] 

x-x  y-y 

^®‘=t(2^)rect(-5^) 


exp[-j2irr  (x^  + y^)  ]rect(-j^)rect  (-jgSi) 


X exp[j2Tir  (xx^  + yyj^)  1 . 


(4-41) 


This  is  our  composite  hologram  mask  transmittance.  Note 

2 2 

that  the  exp[-j2Trr(x  + y )]  term  can  be  interpreted  as 

the  transmittance  of  a thin  convex  lens  with  focal  length 

23 

f determined  by  the  relation 


r » 1/Xf 


(4-42) 


where  X is  the  wavelength  of  the  spatially  coherent  illumi- 
nation. From  the  inequalities  in  Eq.  4-39,  the  focal  length 
must  satisfy 


[max(a,b)l  . 


(4-43) 


For  visible  light,  X <<  max[a,bl. 

The  composite  hologram  approach  just  presented  ob- 
viously has  a Fourier  dual.  The  2-D  system  frequency 
response  is  used  instead  of  the  point  spread  function 
characterization.  The  system  constraints  have,  for  the 


most  part,  different  physical  interpretations.  A con- 
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straint  which  is  identical  in  both  schemes  is  the  band- 
limiting  assvimption  on  the  first  variables  of  the  fre- 
quency response.  That  is,  if  the  point  spread  function 
is  bandlimited  in  x and  y,  then  so  is  the  system  frequency 
response.  The  different  interpretations  appear  in  the 
constraints  on  the  system  input  and  last  two  variables 
of  the  frequency  response.  For  the  point  spread  function 
characterization,  we  required  the  input  to  have  a compact 
support  not  exceeding  the  reciprocal  bandwidths  of  the 
point  spread  function  in  its  first  two  variables.  In  the 
frequency  response  characterization,  we  require  the  input 
spectrum  have  finite  support.  That  is,  the  input  must  be 
bandlimited.  A similar  constraint  applies  to  the  frequency 
response  with  respect  to  its  second  variable. 

The  coherent  processor  capable  of  performing  the  fre- 
quency response  sampling  theorem  is  identical  in  form  so 
that  in  Fig.  4-4  except  that  plane  Pi  and  lens  LI  are 
omitted.  The  bandlimited  input  is  placed  directly  on 
plane  P2  directly  adjacent  to  the  pin  hole  mask.  Plane 
P3  contains  the  appropriate  frequency  response  mask.  Lenses 
L3  and  L4  and  the  rectangular  pupil  in  plane  P4  perform  the 
required  low  pass  filtering.  Our  output  appears  on  plane 
P5. 


i 


r 


r 
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4.3.  Volume  Hologreun  Approach 

The  composite  hologreun  scheme  just  presented  cata- 
logs S2unple  point  spread  functions  side  by  side  on  a holo- 
graphic transmittance.  It  requires,  however,  rather  strin- 
gent compact  support  constraints  on  the  system  input  and 
point  spread  function.  An  alternate  cataloging  scheme, 
not  requiring  these  constraints,  makes  use  of  the  angular 
sensitivity  of  the  volume  hologreun.  Here,  the  system  re- 
sponses are  angle  multiplexed  within  a thick  medium.  As 
we  will  see,  however,  the  scheme  does  not  straightforwardly 
generalize  to  two  dimensions.  It  is  presented  here  with 
the  thought  of  possible  utilization  in  a hybrid  processing 
scheme.  This  topic  is  in  need  of  further  investigation. 

4.3.1.  Diffraction  Efficiency  of  the  Volume  Hologram 

The  diffraction  efficiency  performance  of  the  volume 

hologrcun  has  been  treated  admirably  with  the  coupled  wave 

65 

theory  of  Kogelnik.  Kogelnik's  general  theory  can  be 
found  in  many  optics  texts®^'®®  and  is  presented  via  an 
alternate  derivation  by  Alferness.^^  Kogelnik's  theory 
has  also  been  applied  to  specific  media. 

In  our  treatment  of  volume  hologram  diffraction  effi- 
ciency, however,  we  shall  use  the  one-dimensional  scalar 
optics  model  of  Smith^^  which  follows  from  the  pioneering 
work  of  van  Heerden.^^  We  use  the  scalar  optics  model 
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primarily  due  to  its  relative  analytic  simplicity  in  com- 
parison with  coupled  wave  theory.  The  scalar  optics 
treatment  also  remains  consistant  with  the  other  scalar 
optics  models  used  in  this  thesis. 

An  understanding  o£  the  basic  diffraction  efficiency 
phenomena  of  the  volume  hologram  can  best  be  qualitatively 
understood  by  inspection  of  Fig.  4-5.  In  Fig.  4-5a,  a 
thick  photosensitive  recording  medivim  of  width  t and  re- 
fractive index  n is  illuninated  with  two  coherent  plane 
waves  of  wavelength  X.  The  plane  wave  are  incident  at 
angles  0^  and  6^  where  the  subscripts  r '.pectively  denote 
reference  and  object  waves.  The  thick  media  is  exposed 
and  processed  to  form  a volume  hologram.  We  illuminate 
the  hologram,  as  shown  in  Fig.  4-5b,  with  a plane  wave 
incident  at  an  angle  of  In  general,  three  plane 

waves  will  exit  the  hologram.  The  beam  b^  is  the  zeroth 
order  non-diffracted  component  and  b^  is  the  conjugate  or 
"twin  image"  beam.  The  plane  wave  of  interest  is  b^, 
since,  if  0p  “ then  b^  will  propagate  in  the  same 
direction  as  our  original  object  beam.  The  diffraction 
efficiency  with  respect  to  b^  is  the  ratio  of  the  energy 
of  the  field  amplitude  of  b^  to  the  energy  of  the 
field  amplitude  of  the  incident  playback  beam.  Heur- 
istically  and  in  fact,  the  diffraction  efficiency  is  maxi- 
m\im  when  the  playback  beam  is  in  the  same  direction  as 


Figure  4 . 5 Recording  two  incident  plane  waves  in  a volume 
medium  (a)  The  recording  geometry  (b)  Playback 
geometry  (c)  The  resulting  Bragg  and  extinction 
cones . 


T 
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our  initial  reference  beeun.  That  is,  6, 


e. 


As  6 de- 
P 


P r 

viates  from  0^  say  in  the  clockwise  direction,  the  dif- 
fraction efficiency  is  reduced.  The  angular  deviation 

of  9 from  0 required  to  achieve  zero  diffraction  effic- 
P ^ 

iency  is  appropriately  termed  the  extinction  angle.  Through 

scalar  optics  treatment,  the  extinction  angle,  A0,  is  given 

6 20 

by  the  expression  ' 


2 2 

, n - sxn  0„ 

_ ^ o 

^ / 5 5 7 5 2 

/ n - sin  0^  sin  0^  -/  n -sin  9^  sin  0^ 

(4-44) 

We  assume  that  either  a clockwise  or  counter  clockwise 
deviation  of  A 8 from  0^  will  extinguish  the.  diffracted 
beam,  b^.  While  in  practice,  this  null  actually  occurs, 
rotation  of  the  playback  beam  beyond  the  extinction  angle 
results  in  a finite  diffraction  efficiency.  For  our  pur- 
poses, this  second  order  effect  can  be  safely  ignored  and 
we  can  treat  the  angular  interval 

0 -A0<0  <0  +A0  (4-45) 

r P r 


as  that  in  which  playback  can  be  achieved  with  finite 
diffraction  efficiency.  Note  that  a playback  beam  at 
the  angle  0^  will  also  result  in  maximum  diffraction 
efficiency.  In  future  treatment  however,  we  will  avoid 
this  angular  playback  region  so  that  additional  unwanted 
diffracted  terms  will  not  be  generated. 
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Note  that  the  analysis  to  this  point  has  been  one- 
dimensional.  An  obvious  question  arises:  What  if  the 
playback  bezun  in  Fig.  4-5b  was  incident  on  the  volume 
hologram  in  a direction  not  lying  in  the  plane  defined 
by  two  recorded  beams?  The  diffraction  mechanism  in 
this  more  general  case  can  be  best  understood  with  the 
geometry  of  Fig.  4-5c.  Using  the  recording  beam  direc- 
tions indicated  in  Fig.  4-5a,  we  form  the  "Bragg  cone" 
shown. Next,  visualize  any  plane  passing  through  that 
line  corresponding  to  the  intersection  of  the  volume 
hologram  plane  and  the  plane  defined  by  the  two  recording 
beeims.  For  all  possible  illuminating  plane  waves  with 
propagation  directions  lying  on  this  plane,  the  relative 
majcimum  diffraction  efficiency  occurs  for  those  two  propa- 
gation directions  on  the  Bragg  cone.  Still  on  this  plane, 
we  can  extinguish  the  diffracted  beam  by  rotating  the 
incident  beam  an  angular  interval  corresponding  in  concept 
to  the  extinction  angle.  Thus,  in  a first  order  analysis, 
we  can  speak  of  "extinction  cones"  which,  as  shown  by  the 
dotted  lines  in  Fig.  4-5c,  lie  within  and  without  the 
Bragg  cone.  As  we  shall  see,  this  generalized  view  of  the 
extinction  angle  concept  is  the  very  mechanism  that  pro- 
hibits two-dimensional  generalization  of  the  space-variant 
processing  schemes  to  be  presented. 
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4.3.2.  Space-Variant  System  Representation  i 

We  can  utilize  the  extinction  angle  phenomenon  in- 
herent in  volume  holography  to  implement  space-variant 
systems.  Representation  of  space-variant  systems  with 
volume  holograms  using  the  piecewise  isoplanatic  approxi- 
mation (PIA)  was  earlier  investigated  by  Marks.®  Here, 
we  will  outline  the  implementation  of  a sampling  theorem 
characterization  presented  in  Chapter  II.  Implemen- 

tation of  a discrete  orthonormal  basis  set  response  (DBR) 
characterization  using  the  system  sine  response  is  also 
presented. 

The  sampling  theorem  for  variation  limited  systems 
with  bandlimited.  input  is  repeated  here  from  Eg.  2-84; 

g(x)  = I f (C_)h(x-C_ ;C_)*sinc  2W  x . (4-46) 

n n li  s 

n 

A scheme  for  implementing  this  sampling  theorem  is  pic- 
tured in  Fig.  4-6.  We  excite  our  variation  limited 
system  with  a Dirac  delta  placed  at  input  coordinate 
C = The  corresponding  field  eunplitude  on  the  system 

output  plane  is  thus  h(x-Cj^;  • Also  placed  on  the  system 
output  plane  is  a Dirac  delta  placed  at  x = t a.  The 
desired  value  of  a,  as  we  shall  see,  is  determined  by  the 
spatial  extent  of  the  seunple  line  spread  function. 

The  field  amplitude  on  the  x axis  in  Fig.  4-6  is 


Figure  4.6  Recording  geometry  for  the  volume  hologram  approach. 


given  by 


5(x-Cn-a)  + h(x-Cjj;Cj^).  (4-47) 

This  is  Fourier  transformed  by  the  thin  lens  to  give  a 
field  eunplitude  incident  on  the  thick  photosensitive  re- 
cording medium  of 

exp[-j2TT(Cj^+a)f^l  + H^(f^;Cj^)exp[-j2TTf^5^1  (4-48) 

where  f^  is  the  spatial  dimension  of  the  thick  media's 
plane  divided  by  Xf  and  H (f  ;?)  is  the  system  transfer 

X X 

function  defined  in  Eq.  2-26.  The  photosensitive  medium 
records  the  squared  modulus  of  Eq.  4-48: 

= H^(fx;Cn)exp[j2TTfxa] 

+ [1  + |H^(f^;Cjj)  1^1  (4-49) 

where  denotes  complex  conjugation.  It  is  the  first 
term  in  this  expression  which  is  of  interest. 

Suppose,  now,  we  performed  a number  of  individual 
sequential  exposures  of  this  kind  for  various  values  of 
n.  The  resulting  volume  hologreun  would  then  have  a trans- 


mittance of 
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= I W 

= I H^(f^;?^)exp(j2tTf^a) 
n 

+ I H*  (fj^;5^)exp(-j2Trf^a) 

n 

+ I [1  + 1^].  (4-50) 

n 

Rigorously,  a scale  factor  should  be  included  in  this 
expression  so  that  I(fjj)  is  less  than  or  equal  to  unity 
for  all  f^.  We  omit  it  here  with  no  loss  in  generality. 

With  the  volume  hologram  transmittance  in  Eq.  4-50, 
we  can  implement  the  sampling  theorem  with  the  processor 
pictured  in  Fig.  4-7.  Our  (bandlimited)  input,  f(C),  is 
placed  in  plane  Pi.  Also  placed  in  this  plane  is  a pin 
hole  mask  calibrated  in  accordance  with  the  positioning 
of  the  reference  Dirac  deltas  used  in  the  system  recording 
scheme.  If  we  illuminate  plane  PI  from  the  left  with  a 
normally  incident  plane  wave,  the  field  amplitude  immedi- 
ately to  the  right  of  plane  PI  can  be  modeled  by 

If  (Cn)«(C-Cj^-a)  . (4-51) 

n 

This  is  Fourier  transformed  by  lens  L2  to  give  incident 
on  the  volume  hologram  the  field  amplitude 

I f(5n)exp[-j2Trf^(?^+a)]  . 


m 


(4-52) 


Figure  4.7  Playback  geometry  for  the  volume  hologram  approach. 
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Now,  if  the  hologram  in  plane  P2  were  not  volume,  the  field 
amplitude  immediately  to  the  right  of  plane  P2  would  be 
given  by  the  product  of  Eqs.  4-52  and  4-50.  Since  the 
hologram  is  thick,  however,  angularly  sensitive  diffraction 
efficiency  effects  take  place.  Comparing  the  recording  and 
playback  geometries  in  Figs.  4-6  and  4-7  we  see  that  only  the 
n reference  Dirac  delta  and  the  n^  input  saunple  result 
in  plane  waves  that  are  incident  on  the  volume  medium  at 
identical  angles.  If  we  require  that  the  angular  incidence 
corresponding  to  adjacent  samples  lie  outside  the  equiva- 
lent extinction  angle,  it  then  follows  that  the  only  dif- 
fracted term  corresponding  to  f (5^)  exp  t-j2Trf^  (£;^  + a)  ] is 
Hx(fx»* ^n^ j27rfxa)  . Thus,  we  have  angle  multiplexed  our 
sample  transfer  functions  within  the  emulsion  of  the  vol- 
ume hologram  in  such  a manner  that  they  will  be  accessed 
only  by  appropriately  incident  plane  waves.  It  follows, 
then,  that  the  field  amplitude  exiting  the  volume  holo- 
gram is 

I [f  (Cjj)exp{-j2iTf^(Cjj+a)  }]  tH^(fj^;?^)exp(j27rf^a)  ] 

n 

n 

Also  in  plane  P2,  we  place  a rectangular  pupil,  rect(f  /2W  ) . 

X s 

The  field  amplitude  immediately  to  the  right  of  plane  P2 
is  then  recognized  as  the  Fourier  transform  of  the  sampling 
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theorem  expression  in  Eq.  4-46: 


G(V  = I 


(4-54) 


Lens  L2  performs  a Fourier  transform  on  this  expression. 


Equation  4-46  then  describes  the  field  amplitude  on  plane 


P3  and  our  space-variant  operation  has  been  performed. 


There  are  a number  of  fundamental  limitations  to  this 


scheme  that  should  be  pointed  out.  First  of  all,  with 


reference  again  to  Fig.  4-6,  note  that  there  can  be  no 


overlapping  of  the  Dirac  delta  reference  region  and  the 


seimple  line  spread  function  regions  for  all  n.  This 


must  be  assured  to  avoid  unwanted  crosstalk  on  playback. 


Thus,  the  system  line  spread  function  must  be  "single  sided" 


for  all  C*  The  system  input  must  likewise  be  single  sided. 


But  we  have  assumed  the  input  to  be  bandlimited  and  no 


single  sided  signal  can  be  rigorously  bandlimited.  This 


seeming  inconsistency  can  be  overcome  by  appealing  to  the 


concept  of  the  "essentially"  bandlimited  signal,'  a term 


coined  by  those  applying  sampling  theorem  concepts  to 


causal  (temporal  single  sided)  signals.  A second  funda- 


mental limitation  of  this  scheme  is  its  inherent  one- 


dimensional nature.  That  is,  the  scheme  cannot  be  directly 


extended  into  two  dimensions  since  the  required  diffrac- 


tion  efficiency  mechanism  does  not  generalize  to  two  di- 
mensions. Instead,  we  are  confronted  by  the  previously 
discussed  Bragg  cones. 
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Another  disadvantage,  in  a mathematical  sense,  is 
overscunpling.  In  accordance  with  our  treatment  in  Chapter 
II,  we  are  Scunpling  the  input  at  a rate  2Wg  = 2W^  + 2W^ 
where  2W^  is  the  system  variation  bandwidth  and  2W^  is 
the  bandwidth  of  the  input.  The  actual  minimvim  scimpling 
rate  for  the  input  is  2W^.  We  can,  however,  utilize  the 
minimxun  allowable  (Nyquist)  rate  for  the  input  if  we  angle 
multiplex  the  system's  sine  response  instead  of  its  seunple 
line  spread  function.  Recall,  from  Eq.  2-59,  that  a linear 
system  can  be  characterized  by 


i 

I 


g(x)  * I f(Cj^)S[sinc  2Wj(5-5j^)l  (4-55) 

n 

where,  now,  5^  = n/2W^.  This  relation  is  computationally 
similar  to  the  sampling  theorem  expression  in  Eq.  4-46 
except  for  the  fact  that  the  interpolation  functions  here 
are  sine  responses  instead  of  low-passed  seunple  line  spread 
fxinctions.  In  fact,  the  recording  and  playback  schemes 
in  Figs.  4-6  and  4-7  can  be  used  to  implement  the  sine  re- 
sponse characterization  with  only  two  minor  modifications. 
First,  we  need  to  excite  the  system  with  sinc2W^ in- 
stead of  The  sines  can  easily  be  generated  in  the 
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manner  previously  discussed.  The  second  modification  is 
removal  of  the  rectangular  pupil  on  the  hologram  plane 
in  the  playback  geometry. 

Implementation  of  the  sine  response  characterization 
has  the  scune  fundamental  limitations  as  the  sampling 
theorem  characterization  implementation.  There  are,  how- 
ever, certain  distinct  advantages.  First  of  all,  the 
class  of  representable  systems  is  increased.  That  is, 
the  sampling  theorem  requires  a variation  limited  system 
while  the  DBR  characterization,  of  which  the  sine  response 
characterization  is  a special  case,  is  essentially  inde- 
pendent of  the  system's  pareuneters.  Secondly,  since 
seunpling  is  performed  at  a lower  rate  in  the  sine  response 
characterization,  we  would  expect  to  require  fewer  holo- 
greun  exposures  for  a given  degree  of  output  accuracy. 

The  primary  disadvantage  of  the  volume  hologram  im- 
plementation of  both  the  sine  response  and  sampling 
theorem  characterizations  is  the  resulting  one-dimensional 
resolution  capability.  It  is  possible,  however,  that  they 
could  be  used  as  an  element  in  either  a hybrid  or  more 
elaborate  processing  scheme. 

4.4.  Phase  Coded  Reference  Eeam  Approach 

Phase-coded  reference  beams  have  been  used  for  color 

"TO 

holography  ' and  for  multiplexing  point  source  objects 
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for  digital  storage.  Krile  et.al.  ' have  applied  the 
phase  coded  reference  beeun  idea  to  two-dimensional  space- 
variant  system  characterization.  The  result  is  the  smear- 
ing out  of  unwanted  crosstalk  teirms  into  diffuse  background 
noise.  This  is  in  contrast  to  the  volume  hologram  approach 
where  unwanted  crosstalk  terms  were  filtered  by  the  volume 
hologr2un's  angular  bandpass.  An  overwhelming  advantage  of 
the  phase-coded  reference  beam  approach  over  the  volume 
hologreun  approach,  however,  is  its  capability  of  two-dimen- 
sional system  representation. 

We  will  here  illustrate  utilization  of  phase  coded 
reference  beams  in  implementing  the  s^unpling  theorem  of 
Eq.  2-84.  As  we  shall  see,  mathematical  relations  describ- 
ing its  operation  became  quite  lengthy.  For  this  reason, 
analysis  will  be  performed  in  one  dimension  although  the 
results  are  straightforwardly  extended  into  two  dimensions. 

The  recording  geometry  for  this  scheme  is  pictured  in 
Fig.  4-8  and  closely  resembles  the  structure  of  the  vol- 
ume holograun  recording  scheme  of  Fig.  4-6.  The  system 
is  excited  with  a displaced  Dirac  delta.  The  sample  line 
spread  function,  appears  at  the  system  output. 

Placed  on  this  same  plane  is  a random  diffuser  with  trans- 
mittance d(x)  which  is  displaced  clear  of  the  system  out- 
put plane.  We  illuminate  the  diffuser  from  the  left  with 
a coherent  plane  wave  over  the  interval  5^+a-T<x<Cj^+a+T 


Figure  4.8  Recording  geometry  for  the  phase  coded  reference  beam  approach 
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where  a denotes  the  displacement  of  the  diffuser  and  2T  is 
the  interval  width.  The  field  amplitude  in  the  front  focal 
plane  of  lens  Ll  is  thus' 


h(x-Cn;Cn)  + d^(x) 


(4-56) 


where 


dn(x) 


d(x)rect[- 


x-Cn-a 


-]  . 


(4-57) 


This  is  Fourier  transformed  by  the  thin  lens  to  give  in- 
cident on  the  photosensitive  media  the  field  cunplitude 


(4-58) 


where  Djj(fjj)  = transmittance  of  the  re- 
sulting hologram  is  the  squared  modulus  of  this  expres- 
sion: 


* l''x«x'5n'l"  * lo„<Vl^ 


It  is  the  first  term  in  this  expression,  of  course,  which 
^ ' I is  of  interest. 

Suppose,  then,  that  a number  of  such  exposures  were 
recorded  on  the  same  media  for  various  values  of  n.  The 
resulting  hologram  transmittance  is  then 

* 

■ 
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Note  that,  in  order  for  the  same  portion  of  the  diffuser 
not  to  be  used  in  more  than  one  of  these  exposures,  that 
we  must  require  Let  us  assume  adjoined  pat- 
ches so  that  C ® 2T  for  all  n. 

n+i  n 

Consider  next  the  playback  geometry  in  Pig.  4-9. 

The  diffuser,  lens  LI,  and  the  hologram  are  left  in  the 
same  position  as  in  Fig.  4-8.  To  the  left  of  the  diffuser 
an  input  transmittance,  f(C),  is  sandwiched  between  two 
flys  eye  lens  arrays.  Each  lenslet  is  assumed  to  have 
width  2T.  The  left  most  lenslet  is  illuminated  from  the 
left  by  a coherent  plane  wave.  Thus,  incident  on  our 
input  transmittance  is  the  field  amplitude. 

I 6(C-Cjn-a).  (4-61) 

m 

That  is,  each  of  the  lenslets  in  the  left  most  array  is 
assiamed  to  focus  the  plane  wave  incident  on  its  pupil  to 
a point  source.  Equation  4-61  then  follows  as  a result 
of  the  discussion  on  Dirac  delta  generation  discussed 
previously  in  this  chapter. 

Equation  4-61  now  multiplies  the  input  transparency 
to  give 

ICl 


(4-62) 


Playback  geometry  for  the  phase  coded  reference  beam  approach 
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Each  weighted  Dirac  delta  in  this  sum  is  then  laterally 
transformed  by  the  right  most  flys  eye  array  to  yield  a 
series  of  parallel  weighted  plane  waves.  Thus,  the  field 
aunplitude  incident  on  the  diffuser  is 

I f(£)rect (1-^22-^).  (4-63) 

m 2T 

Thus,  in  the  front  focal  plane  of  lens  Ll  is  the  distribu- 
tion 

I f(C^p)rect(''  d(x) 

m 

m 

This  is  Fourier  transformed  by  lens  Ll  to  give  the  field 
amplitude  incident  on  the  hologram  plane; 


I f(CjDj^(f^)rect(2Tf^)  (4-65) 

m 

where  we  have  placed  a rectangular  pupil  in  the  holograun 
plane  for  low  pass  filtering.  This  multiplies  the  holo- 
gr2un  transmittance  in  Eq.  4-60  to  give  immediately  to  the 
right  of  the  hologram  the  distribution 


D(£^) 


'5  I 

m r. 


= T,(f^) 


^2(V 


(4-66) 
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where 


T3(V=  I I 

m n 


TjCV-  I I *(5„)°m<V«J'*x'«n' 

m n 

T>^  ( f jj)  exp  [ j 2irf  rect  ( 2Tf (4-67) 


Ti<fx).  I I *(V'>m<V«x<^x'5„> 


m n 


D*  ( f^)  exp  [- j 2iTf ^5^]  rect  ( 2Tf  ^ 


Lets  interpret  each  of  these  stuns  individually.  T2(*) 
is  the  zero  order  "through"  term  which  results  from  the 
incident  beam  propagating  unaltered  in  direction  through 
the  hologram  transmittance.  Term  T2(*)  is  the  amplitude 
of  the  so  called  conjugate  or  "twin  image"  diffracted 
term.  As  is  shown  in  Fig.  4-9,  this  term  will  be  gen- 
eral exit  the  hologram  at  a rather  steep  angle  for  the 
chosen  recording  geometry.  It  is  that  we  are  in- 

terested in.  As  shown,  the  general  propagation  direction 
of  the  wave  with  amplitude  Tj^(*)  is  down  the  system  axis. 
Let  us  then,  suppose  (more  or  less  realistically)  that  we 
can  either  suppress  the  effects  of  T2 ( • ) and  ( * ) by 
filtering,  or  that  their  appearance  on  the  processor  out- 
put plane  does  not  overlap  with  the  Tj^(«)  term.  With 
this  assumption,  we  can  restrict  attention  to  T^(*). 
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Consider,  then,  the  effect  of  lens  L2  which  performs 
a Fourier  transform  on  to  give  the  amplitude  tj^(x) 

on  the  processor  output  plane. 

ti(x)  = y.  I f (5j„)h(x-Cj^;Cj^)*sinc  ^ 
n m 

*[d„(x)*d„(x)].  (4-68) 

n m 

Here  we  have  made  use  of  the  convolution  and  autocorrela- 

2 8 

tion  theorems  of  Fourier  transform  theory. 

To  this  point,  we  have  made  no  comment  concerning  the 
character  of  the  diffuser.  Comparison  of  Eg.  4-68  with 
the  seunpling  theorem  expansion  in  Eg.  2-84,  however,  makes 
clear  the  desired  property  of  diffuser.  Specifically,  with 
1/2T  = 2W  , we  would  like 

S 

d„(x)ird„(x)  » 6(x)5^^  (4-69) 

n m nm 

where  is  the  Kronecker  delta.  That  is,  for  n = m, 
nm 

the  autocorrelation  of  a single  patch  of  the  diffuser 
should  be  peaked  like  a Dirac  delta.  Otherwise,  for  n ^ m, 
the  crosscorrelation  of  two  diffuser  patches  should  ideally 
be  zero.  If  this  condition  could  be  met,  then  Eg.  4-68 
would  reduce  to 

t3^(x)  = y f(Cjj)h(x-Cn;Cn)*sinc  . (4-70) 

n 

Outside  of  a proportionality  constant,  this  is  the  ex- 
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pression  for  the  sampling  theorem  in  Eq.  2-84  for  2W^»1/2T. 
Thus,  our  space-variant  operation  has  been  performed  and 
tj^(x)  * g(x)  . 

Unfortunately,  implementing  a diffuser  which  has  the 
property  of  Eq.  4-69  is  obviously  not  rigorously  possible. 
Thus,  we  must  settle  for  a diffuser  which  exhibits  this 
property  to  a good  approximation.  Such  a diffuser  would 
exhibit  peaked  autocorrelation  for  all  patches  and  a broad 
low  level  crosscorrelation  for  each  unequal  patch  pair.  A 
direct  consequence  of  a latter  characteristic  is  the  occur- 
ence of  cross  talk  terms  in  Eq.  4-68.  These  terms,  however, 
will  be  smeared  out  into  diffuse  background  noise.  An 
investigation  into  various  diffusers  for  use  in  the  phase 
coded  reference  beeun  scheme  is  presently  being  undertaken 
by  Krile  et.al.^® 

As  a final  note  in  this  section,  we  again  stress  that, 
unlike  the  volume  hologram  scheme,  the  phase  coded  ref- 
erence beam  approach  can  be  implemented  in  two  dimensions 
and  thus,  in  this  regard,  is  much  more  powerful. 

4.5.  Discussion 

In  this  chapter,  we  have  discussed  system  response 
generation  and  two-dimensional  space-variant  system  repre- 
sentation. Certain  types  of  system  responses  were  shown 
to  be  straightforwardly  performed.  These  include  the 
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system  sine  and  rect  responses.  The  system  point  spread 
function  was  also  shown  to  be  easily  generated  for  both 
input  and  output  sampling  cases. 

Three  methods  of  2-D  system  representation  were  pre- 
sented. Let  us  briefly  contrast  the  advantages  and  dis- 
advantages of  each. 

The  composite  hologram  implementation  for  systems 
whose  point  spread  functions  are  bandlimited  in  their 
output  variables  is  probably  the  most  straightforwardly 
implementable  of  the  three  schemes.  The  required  mask  can, 
in  principle,  either  be  obtained  through  direct  optical 
recording  or  from  computer  generated  hologram.  The  pri- 
mary drawback  of  this  scheme,  however,  is  the  extremely 
narrow  allowable  class  of  systems  and  inputs  to  which 
the  scheme  can  be  applied. 

The  volume  hologram  approach,  although  containing  a 
two-dimensional  flavor,  permits  straightforward  represen- 
tation only  of  one-dimensional  systems.  In  this  regard, 
the  one-dimensional  processors  in  Chapter  III  would  seem 

to  be  superior.  Another  fundamental  drawback  is  the 

✓ 

necessity  of  direct  optical  recording.  That  is,  no  scheme 
for  computer  generated  volume  holograms  has  been  formu- 
lated to  date.  The  beauty  of  the  volume  hologram  approach, 
however,  is  the  clean  manner  in  which  it  filters  out  un- 


wanted crosstalk  terms.  Utilization  of  the  volume  holo- 
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gram  in  future  hybrid  two-dimensional  space-variant  pro- 
cessors is  a distinct  possibility.  Its  use  in  supressing 

part  of  the  crosstalk  in  the  phase  coded  reference  beam 

22 

approach  has  also  been  suggested. 

In  terms  of  application  to  the  widest  class  of  2-D 
systems,  the  phase  coded  reference  beam  approach  seems 
to  be  the  best  of  the  three  schemes.  The  required  holo- 
graphic mask  can  be  either  directly  recorded  or  obtained 
from  a computer  generated  hologram.  The  primary  disad- 
vantage of  the  scheme  is  the  appearance  of  crosstalk  terms 
in  the  form  of  diffuse  background  noise.  This  problem, 
however,  can  be  minimized  by  construction  of  more  optimal 
diffusers. 


CHAPTER  V 


5.  CONCLUSIONS 

In  this  report  we  have  introduced  new  views  on  linear 
system  characterization  theory  and  have  proposed  various 
schemes  for  coherent  linear  space-variant  processing.  A 
recap  of  these  results  and  reflections  on  possible  future 
work  are  in  order. 

The  continuum  orthonormal  response  system  characteri- 
zation introduced  in  Chapter  II  offers  a generalized  view 
of  linear  system  description.  The  familiar  line  spread 
function  and  frequency  response  characterizations  were 
shown  to  be  special  cases.  Of  more  practical  interest 
are  the  three  system  characterizations  wherein  constraints 
are  placed  on  the  system  and/or  input.  The  piecewise 
isoplanatic  approximation  (PIA)  is  such  a case.  The  con- 
straint here  is  that  the  system  line  spread  function 
change  "slowly"  in  shape  over  an  isoplanatically  modelevi 
interval . 

Additional  investigation  of  the  PIA  could  be  per- 
formed in  several  areas.  One  as  of  yet  vinsolved  problems 
is  the  convergence  of  the  PIA  to  the  true  system  output  as 
each  patch  grows  arbitrarily  small  and  the  patch  density 
becomes  arbitrarily  large.  Interestingly,  the  energy 
(integral  of  the  modulus  squared)  of  the  PIA  and  true 
outputs,  for  certain  systems  even  in  the  limit  is  not 
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equal. ^ More  drastically,  the  PIA  of  a finite  energy  output 

7 

might  have  an  infinite  energy.  It  is  this  author's  conjec- 
ture that  only  very  loose  constraints  can  be  placed  on 
the  convergence  of  the  PIA  for  the  most  general  case. 

Possibly  tighter  bounds  can  be  placed  on  a specific  system 
or  system  class.  A similar  area  of  investigation  would 
be  the  study  of  filtering  the  PIA  output  to  improve  con- 
vergence . 

A second  interesting  yet  unsolved  area  of  investigation 

concerning  the  PIA  is  definition  of  the  isoplanatic  patch. 

26 

Lehman  and  Paris  have  offered  such  a definition  which, 

unfortunately,  seems  limited  in  application  to  imaging 

77 

systems.  Arsenault  and  Brousseau  have  also  noted  defic- 
iencies in  the  Lohman-Paris  definition.  We  can  conclude 
from  Lohman  and  Paris,  however,  that  a definition  of  the 
isoplanatic  patch  need  be  only  system  dependent  (and  thus 
input  independent) . Possibly  the  concept  of  the  variation 
spectrum  introduced  in  this  thesis  offers  the  basis  for  a 
new  more  universal  definition  of  the  isoplanatic  patch. 

A second  constrained  system  characterization  pre- 
sented in  Chapter  II  was  the  discrete  orthonormal  basis 
element  response  (DBR)  characterization.  Here,  the  system 
is  completely  defined  for  an  input  class  by  cataloging 
the  system  response  to  each  element  in  an  orthonormal 
basis  set.  The  inputs  are  limited  to  that  class  spanned  by 
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the  basis  set.  The  sine  and  rect  responses  were  shown  to 

be  particularly  straightforwardly  implementable  in  coherent 

processor  applications  respectively  for  band*  and  space- 

limited  inputs.  Due  to  vignetting  and  finite  pupil  effects, 

only  bandlimited  and  space-limited  signals  are  present 

78 

within  any  space-variant  processor  (within  the  treatment 

79 

of  scalar  optics  ) . 

As  with  the  PIA,  an  interesting  study  of  the  DBR 
characterization  is  convergence.  A given  input  conver- 
gence does  not  necessitate  a given  output  convergence. 

For  excunple,  as  our  system,  consider  a Fourier  transformer 
with  bandlimited  input.  Using  the  sine  basis  set,  our 
input  is  reconstructed  via  the  conventional  Whittaker- 
Shcuinon  sampling  theorem  which  displays  uniform  converg- 
ence.^^ The  system  output,  which  will  be  reconstructed 

24 

via  a Fourier  series,  displays  Gibbs  phenomena  and 
thus  does  not  uniformly  converge  (convergence  is  in  the 
mean) . Thus  the  input  and  output  display  different  con- 
vergence properties.  The  converse  holds  for  a Fourier 
transformer  with  space-limited  input  using  the  rect  basis 
elements.  It  would  seem  that  for  "physically  realizable 
systems”  with  finite  energy  inputs  that,  if  our  input 
space  were  spanned,  then  so  would  our  output  space.  We 
could  then  expect,  it  seems,  at  least  mean  square  converg- 
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Another  possible  area  of  investigation  of  the  DBR  is 
the  best  choice  of  a basis  set.  A criteria  on  which  the 
basis  set  could  be  chosen  is  output  convergence.  That  is, 
for  a given  allowable  output  error,  we  would  want  to  choose 
that  basis  set  whose  expansion  converged  quickest  on  the 
output  with  fewest  terms.  Fabrication  and  use  of  ortho- 
normal  element  masks  for  characterization  of  coherent  pro- 
cessors would  also  be  an  interesting  undertaking. 

The  final  constrained  system  characterizations  pre- 
sented in  Chapter  II  were  the  sampling  theorems.  The 
S2unpling  theorems,  in  principle,  seem  to  be  the  most  di- 
rectly applicable  to  coherent  processor  implementation. 

The  system  sine  response,  herein  denoted  as  a DBR  char- 
acterization, could  also  be  argued  to  belong  to  the 
Seunpling  theorem  class. 

Besides  direct  application  to  coherent  processing, 
the  Seunpling  theorems  are  also  useful  in  other  areas. 

In  circuit  theory,  they  give  rise  to  methods  of  time- 
variant  filter  synthesis. Also,  in  treatment  of  linear 
systems  on  the  digital  computer,  the  sampling  theorem 
results  give  an  idea  of  the  sampling  rate  required  for 
good  approximation.  Although  such  sampling  rates  have 

g 

been  discussed  elsewhere  in  a somewhat  more  limited  scope, 
no  analytic  expressions  for  the  resulting  system  charac- 


terization have  been  previously  presented. 
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A curious  question  arises  concerning  utilization  of 
the  proposed  sampling  rates.  As  is  discussed  in  the  main 
text,  sampling  is  performed  on  all  members  at  a rate  above 
the  minimum  allowable  sampling  rate  for  each  member  indi- 
vidually. Seemingly,  we  gain  by  this  overseunpling  a 
simple  computational  form  that  is  a direct  discrete  case 
of  the  superposition  integral.  If  sampling  of  each  member 
is  performed  at  its  corresponding  minimum  sampling  rate, 
the  resulting  expressions  are  vastly  more  complex.  How- 
ever, in  actual  digital  implementation,  this  manner  of 
S2unpling  would  seem  to  require  less  memory  for  a given 
degree  of  allowable  error.  Thus,  a tradeoff  between  the 
two  sampling  schemes  seems  to  be  simplicity  vs.  required 
computer  memory.  A more  definitive  study  of  this  tradeoff 
in  term  of  space-bandlimited  products  would  be  an  inter- 
esting undertaking. 

There  are  also  other  possible  avenues  of  investigation. 
One  is  effects  of  trunctation.  It  would  seem  that  this 
would  be  a straightforward  extension  of  work  already  per- 
formed on  trxuication  effects  in  the  conventional  Whittker- 
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Shannon  sampling  theorem.  Another  interesting  topic 

83 

would  be  a study  of  the  effects  of  jitter  (irregularly 
spaced  saunple  points)  . 

All  of  the  constrained  system  characterizations  in 
Chapter  II  are  deterministic  and,  in  principle,  should 


hold  relatively  true  in  the  presence  of  low  level  noise. 

One  can,  however,  envision  cases  where  an  intermediate  noise 
level  would  give  rise  to  non-negligible  error.  The  allow- 
ance of  this  occurance  uncovers  an  entire  new  area  of  pos- 
sible study.  For  space-variant  systems  that  are  essentially 
deterministic,  however,  the  system  characterizations  as 
presented  in  Chapter  II  offer  a new  wealth  of  both  analytic 
and  practical  treatment. 

Let  us  now  move  on  to  Chapter  III  in  which  numerous 
schemes  for  one-dimensional  space-variant  processing  were 
presented.  All  processors  considered  make  use  of  either 
an  astigmatic  or  Fourier  transform  type  processor.  The 
class  of  space-variant  operations  that  can  be  performed 
by  these  processors  is  extremely  broad  being  fundamentally 
limited  only  by  mask  implement ability.  Even  in  this  re- 
gard, we  have  a degree  of  flexibility.  That  is,  if  a mask 
for  a given  operation  is  not  implementable  for  a given 
processor,  it  might  be  implementable  for  an  alternate 
processor.  As  was  shown,  the  one-dimensional  space-var- 
iant processors  have  a variety  of  applications  includ- 
ing coordinate  distortion,  integral  transform  display, 
convolution,  correlation,  and  cross-spectral  density  dis- 
play. There  were  also  two  serendipities  of  the  processor 
schemes  wherein  the  entire  processor  output  plane  was 
utilized.  These  were  the  Laplace  transform  and  ambiguity 
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function  display  schemes.  The  simplicity  of  the  ambiguity 
function  processor  has  allowed  for  its  inclusion  in  a col- 
lection of  projects  designed  to  introduce  the  undergrad- 
uate Electrical  Engineer  to  coherent  optical  processing. 

A major  result  of  1-D  processor  study  is  the  lensless 
processor.  Here,  the  space-variant  operation  is  performed 
by  a single  mask  directly  adjacent  to  which  we  place  our 
1-D  input.  It  is  hard  to  imagine  a coherent  processor  in 
a less  compact  form. 

Despite  the  diversity  already  demonstrated  for  the 
1-D  processors  group,  a number  of  fundaimental  investiga- 
tions remain.  At  best,  our  treatment  of  the  processors 
was  only  first  order.  We  completely  neglected  such  fun- 
damental parameters  such  as  finite  pupil  size  and  vignet- 
ting which  in  certain  instances  can  cause  significant 
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degradation  in  processor  performance.  Mueller  and  Carlson 
for  example  have  shown  that  bandlimiting  effects  of  the 
astigmatic  processor  impose  fundamental  limitations  on 
the  performance  of  the  Laplace  transform  processor.  As 
such,  it  would  be  instructive  to  perform  a general  inves- 
tigation on  the  first  order  degrative  effects  of  the  var- 
ious 1-D  processors. 

Mask  fcHsrication  is  another  fundamental  area  in  need 
of  further  research.  Three  possibilities  are  1)  optically 
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holographically  recorded  masks,  2)  masks  generated  by  CRT 
display  and  3)  computer  generated  hologreuns.  A combi- 
nation of  these  methods  might  also  be  employed. 

One  of  the  possible  potential  uses  of  the  1-D  pro- 
cessors is  in  filtering  of  temporal  signals.  For  example, 
1-D  processors  have  distinct  advantages  over  the  conven- 
tional 2-D  convolution/correlttMon  coherent  processors 
where  movement  or  Fourier  encoding  is  required.  Utiliz- 
ing one-dimensional  electro-optical  transducers  to  provide 
the  1-D  processor  input,  such  operations  could  be  performed 
in  near  real  time. 


Of  more  fundeuaental  importance  in  temporal  signal 
processing  is  the  concept  of  "time-variant  filtering." 
Through  discussions  with  colleques  and  inspection  of  the 
literature,  it  is  this  author's  opinion  that  this  topic 
has  not  yet  even  begun  to  be  investigated.  There  are 
even  recent  publications  in  which  rather  elementary 


time-variant  filters  are  being  analysed  let  alone 
applied  in  signal  processing.  As  such,  it  would  seem 
that  application  of  the  1-D  processor  to  temporal  signal 
processing  can  be  applied  to  time-variant  filtering  only 
when  the  advantages  and  applications  of  such  filters  are 
known.  We  should  here  stress  that  these  comments  apply 
only  to  use  of  the  term  "time-variant  filtering"  as  a 
generalization  of  the  concept  of  time-invariant  filtering 
by  electric  circuits. 
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Lastly,  let  us  reflect  on  Chapter  IV  wherein  various 
possible  schemes  for  general  two-dimensional  space-variant 
processing  were  presented.  The  first  scheme  (the  composite 
hologram  approach)  was  a straightforward  implementation  of 
a sampling  theorem  applicable  when  the  system  point  spread 
fxinction  is  bandlimited  in  its  output  variables.  The  method, 
however,  places  severe  constraints  on  the  support  of  the 
input  and  is  thus  applicable  in  only  a small  number  of  cases. 

The  volume  hologram  scheme  was  discussed  next  and  was 
shown  to  be  fundamentally  one-dimensional  in  its  application 
to  space-variant  processing.  The  ability  of  the  volume 
hologram  to  store  a number  of  wavefronts  with  angular 
accessibility,  however,  still  leaves  the  possibility  for 
its  use  in  a hybrid  two-dimensional  processor.  Its  use 
in  partial  supression  of  diffuse  background  noise  in  the 
phase-coded  reference  beam  scheme  is  also  a distinct  pos- 
sibility. 

Utilization  of  phase-coded  reference  beams  is,  to  date, 
the  most  promising  method  for  two-dimensional  coherent  > 

space-variant  processing.  The  fundamental  drawback  to 
this  scheme  is  the  smearing  of  unwanted  crosstalks  terms 
into  diffuse  background  noise.  Although  this  effect  can-  ' 

not  be  totally  eliminated,  it  can  be  minimized  by  the  use 
of  codes  which  display  peaked  Dirac  delta  type  autocorre- 
lation and  nearly  zero  crosscorrelation.  An  investiga-  ’ 
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APPENDIX  A 


SPATIAL  AND  FREQUENCY  INVARIANCE 

23  . . 24 

The  concept  of  space-invariant  and  time-invariant 
systems  is  well  established.  For  linear  space- invariant 
systems,  the  system  line  spread  function  takes  on  the 
form 

h(x-C;C)-^h(x-C).  (A-1) 

As  discussed  in  Chapter  1,  2-D  linear  space-invariant 
systems  can  be  easily  implemented  with  coherent  optical 
processors.  Furthermore,  linear  space-invariant  systems 
can  be  straightforwardly  analyzed  through  application  of 
conventional  Fourier  transform  theory. 

There  are  some  linear  systems  whose  classification  as 

77 

variant  or  invariant  is  not  clear.  Arsenault  and  Brousseau 
have  discussed  a class  of  lineai^  systems  which  are  invar- 
iemt  for  a certain  input  class.  An  example  of  a line 
spread  function  for  such  a system  is 

h(x-C;C)  = h^(x-C)rect(|^).  (A-2) 

Due  to  the  rect  term,  this  line  spread  function  does  not 
rigorously  conform  to  the  invariance  criterion  in  Eq.  A-1. 

If,  however,  we  limit  our  input  signal  class  to  those 
functions  identically  zero  outside  the  interval  |5|  T, 
then  the  rect  term  is  superfluous.  Arsenault  and  Brous- 
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seau  have  termed  such  systems  "quasl-linear"  although 
"quasi-invariant"  might  be  more  appropriate. 

A second  system  class  which  has  a vague  classifica- 
tion contains  those  linear  systems  with  line  spread  func- 
tions of  the  form 


h(x-5;5)-^h(x-MC) 


(A- 3) 


where  M is  a fixed  constant.  Although  similar  in  form, 

this  line  spread  function  does  not  rigorously  conform  to 

the  invariance  criterion  in  Eq.  A-1.  A system  considered 

elsewhere  in  this  thesis  which  is  contained  in  this  system 

class  is  the  non-unity  magnification  imaging  system  [See 

23 

Eq.  3-8] . Goodman  classifies  the  magnifier  as  space-in- 

53 

variant  while  Sawchuk  calls  it  "rigorously"  space-variant. 

Alternate  linear  system  classification  categories 
are  frequency-variant  and  frequency-invariant  systems.  A 
linear  system  is  said  to  be  frequency-invariant  if 


k ( x-cv ; V ) -►k  ( x-c  v ) 


(A-4) 


where  c is  a consteuit  and  k(*,*)  is  the  system  frequency 
response  defined  in  Eq.  2-19.  If  the  frequency  response 
does  not  conform  to  Eq.  A-4,  the  system  is  said  to  be 
frequency- variant.  The  term  "frequency"  can  here  apply 
to  either  spatial  or  temporal  frequency. 

An  example  of  a temporal  frequency- invariant  system 
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is  a delay  line  whose  delay  time,  x,  is  proportional  to 
the  applied  frequency  (v) . An  example  of  a spatial  fre- 
quency-invariant system  is  the  familiar  thin  lens  Fourier 
transformer.  Using  Eq.  2-22,  the  frequency  response  of 
the  Fourier  transformer  is 

k(x-cv;v)  = /"^exp x?] 

= Xf6(x-xfv].  (A-5) 

Note  that  the  constant  c is  here  determined  by  system 
pareuneters 


c * Xf  . (A-6) 

We  can  easily  show  that  a space-invariant  system  can 
never  be  frequency- invariant.  Using  Eq.  2-22  as  applied 
to  a space-invariant  system,  we  have 

k(x-cv;v)  = /“^h(x-S)e^^^^'’dC 

where  * Fjjlh(x)].  Clearly,  this  relation  can  never 

confrom  to  the  frequency  invariance  criterion  by  Eq.  A-4. 
(Note  that  the  constant  c is  here  not  specified  by  any 
system  parameters.)  We  can  similarly  show  that  a fre- 
quency-invariant system  can  never  be  space- invariant. 
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In  two-dimensional  coherent  processors,  a space-invar- 
iant system  can  be  straightforwardly  be  made  into  a £re- 
[ quency- invariant  system  (or  visa-versa)  simply  by  placing 

i 

F'  a thin  spherical  lens  a focal  length's  distance  f in  front 

t 

of  the  system's  input  plane.  The  input  plane  is  rede- 
fined  to  be  the  lens'  front  focal  plane.  Then,  for  an 
[ input  exp(j2TrCv),  the  system  output  is 

I k(x-cv;v)  = Xfh(x-Xfv)  (A-8) 

j ^ where  h(»)  is  the  original  space-invariant  system's  im- 

pulse response.  Clearly,  the  augmented  system  is  fre- 
quency-invariant with  c = Xf. 

In  the  definition  of  frequency-invariance  in  Eq.  A-4, 
the  constant  c is  assumed  to  have  dimension  of  inverse 
I frequency  squared.  This  retains  dimensional  consistancy 

between  the  spatial  frequency  v and  the  spatial  variable 
X.  Inclusion  of  such  a constant  in  the  spatial  invariance 
definition  in  Eq.  A-1  was  not  required  since  x and  C are 
dimensionaly  consistant.  Note,  however,  that  if  we  had 
included  an  arbitrary  constant  in  the  spatial  invariance 
. criterion,  then  we  would  arrive  at  Fq.  A-3.  Thus,  we 

I might  logically  adopt  Eq.  A-3  as  a generalized  definition 

j of  spatial  invariance. 

I We  have  here  considered  space- (time-)  and  frequency- 


invariant  systems.  A term  used  in  the  literature  which 


APPENDIX  B 


UTILIZING  MINIMUM  SAMPLING  RATES  IN  A SPACE-VARIANT  SAMPLING 
THEOREM 

The  S2unplin9  theorem  expressions  in  Eqs.  2-83  and  2-84 
are  not  optimum  in  the  sense  of  utilizing  maximum  allowable 
Seunpling  intervals.  That  is,  we  are  seunpling  both  the  in- 
put and  line  spread  function  at  a rate  of  2W,  while  the 
minimum  required  sampling  rates  are  2W^  and  2W^,  respect- 
ively. As  will  be  shown,  however,  the  resulting  expres- 
sion employing  these  minimum  sampling  rates  is  rather  un- 
attractive for  computation  and  implementation  purposes. 

Consider,  then,  the  following  seuapling  theorem  ex- 
pansion of  a space-variant  system's  line  spread  function: 

h(x;C)  = I h(x;C  )sinc2W^(C-C  ),  (B-1) 

where  2W^  is  the  variation  bandwidth  and  5g=q/2W^.  One 
may  similarly  apply  the  sampling  theorem  to  the  system  in- 
put to  give 

f(C)  » If(C  )sinc2Wj(5-C  ),  (B-2) 

p P ^ 

where  Cp  = p/2W^  and  2W^  is  the  input's  bandwidth.  Sub- 
stituting Eqs.  B-1  and  B-2  into  Eq.  2-25  gives 
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r V p q ^ V 


(B-3) 


X exp(-j2TTf^5^)  ]*  [rect(2^)exp(-j2Trf^5^)  ] }. 


x-q 


x-p' 


where  "*"  denotes  convolution.  Equation  B-3  is  identical 
to  Eq.  2-84  yet  employs  larger  sampling  intervals.  The 
above  relationship,  however,  has  the  disadvatage  of  not 
assigning  each  sample  input  value  to  a single  correspond- 
ing sampled  line  spread  function.  Note  that  the  two 
convolving  rects  in  Eq.  B-3  give  an  upper  bound  of  2W  =2W_ 

+ 2W^  on  the  output  spectrum's  bandwidth.  This  is  the  same 
constraint  contained  in  Eq.  2-83. 

Inverse  transforming  Eq.  B-3  gives 

g(x)  = I If(£  )h(x;E  ) 

P ^ ^ {B-4) 

*[sinc  2W^(x-Cp)sinc  2W^(x-Cq)]. 

This  is  the  spatial  domain  sampling  theorem  expansion  for 
minimum  sampling  rates. 

We  now  consider  the  version  of  Eq.  B-4  which  arises 
when  the  line  spread  function  is  bandlimited  in  x with 
low  pass  bandwidth  of  2W^.  In  this  case,  we  can  write 

h(x;Cq)  = 5|.h(xj.;?q)  sine  2W^(x-x^) 


(B-5) 
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where  x = r/2W  . Siibstituting  into  Eq.  B-4  gives 

JT  X 

g(x)  = I I I 

where  our  interpolation  function  is 


Ipqr(3c)  “ sine  2W^(x-x^) 


(B-7) 


*[sinc  2W^(x-Cp)sinc  2W^(x-Cq)]. 


This  relationship  is  obviously  less  computationally  attrac- 
tive than  that  in  Eq.  2-83.  Note  that  evaluation  of  Eqs. 
B-6  and  B-7  at  x^^^  = m/2W  gives  Eqs.  2-101  and  2-102. 


